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SUMMARY 

Jhe  prime  purpose  of  this  research  project  is  to  apply 
the  systems  approach  and  the  technique  of  computer  programming 
in  an  attempt  to  improve  the  process  of  tunneling  by  rapid 
excavation  methods.  One  specific  objective  is  the  optimiza¬ 
tion  of  the  materials  handling  function  for  tunneling  systems. 
This  report  contains  information  on  the  methods  of  simulation 
on  the  digital  computer,  the  development  of  the  computer  model 
from  the  basic  concepts,  the  guide  to  the  use  of  the  program, 
and  the  computer  program  itself. 

The  review  of  simulation  methods  contains  a  basic  defini¬ 
tion  of  simulation  and  a  description  of  the  most  useful  types 
of  models  for  digital  computers.  Primary  attention  was  given 

r 

to  the  stochastic  and  deterministic  conceptual  models  which 
are  used  in  the  simulation  of  tunneling  systems  and  to  methods 
of  updating  the  computer  simulation  time  variable.  In  addi¬ 
tion,  a  discussion  of  computer  languages  for  possible  use  in 
a  simulation  model  is  presented.  The  FORTRAN  language  was 
chosen  for  use  in  the  model  based  primarily  upon  its  wide 
acceptance  and  its  familiarity  to  potential  users. 

The  description  of  the  computer  model  contains  a  synopsis 
of  the  objectives  of  the  model  as  well  as  an  outline  of  the 
concepts  used  in  the  simulation  program.  The  discussion  of 
the  specific  concepts  applied  is  divided  into  sections  dealing 
with  the  individual  unit  operations:  muck  generation,  mater¬ 
ials  handling,  roof  support,  and  environmental  control.  An 
additional  section  deals  with  the  general  concepts  used 
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throughout  the  program.  The  most  detailed  discussions  deal 
with  the  muck  generation  and  the  materials  handling  sub¬ 
systems  which  are  the  most  complex  unit  operations  to  be 
modeled  in  the  program.  The  materials  handling  methods 
receive  the  greatest  analysis  since  they  are  the  most  complex 
from  a  systems  standpoint  and  are  the  most  difficult  to 
simulate . 

The  attempts  at  testing  the  computer  program  are  des¬ 
cribed  in  a  separate  section  of  the  report.  At  presents,  the 
program  has  been  debugged  and  the  logic  and  behavior  of  the 
model  during  simulation  have  been  studied  using  data  ob¬ 
tained  primarily  in  the  field.  No  attempts  have  yet  been 
made  to  check  the  accuracy  of  the  computer  model  as  this 
phase  of  the  testing  program  is  scheduled  in  the  near  future. 

The  final  section  of  this  report  is  a  users'  guide  to 
the  computer  program,  a  list  and  description  of  the  most 
important  variables  contained  in  the  program,  and  a  listing 
of  the  computer  program. 
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SIMULATION  BACKGROUND 

The  term  "simulation"  is  used  quite  frequently  in 
modern  technical  literature  as  the  methods  of  computer  model¬ 
ing  have  become  more  widely  applied  and  accepted.  This 
section  of  the  report  provides  a  brief  outline  of  the  simu¬ 
lation  methods  and  languages  and  defines  the  simulation  terms 
used  throughout  the  report. 

Definition  ol  Simulation 

Simulation  has  been  defined  in  numerous  ways,  but  a 
definition  that  appears  in  a  book  by  Pritsker  and  Balintfy  (9)* 
appears  to  be  most  applicable  here.  They  have  said  that 
"Simulation  is  the  use  of  a  model  to  study  a  problem."  This 
simple  yet  concise  definition  describes  very  well  the  ap¬ 
proach  used  in  this  project.  Our  problem  is  to  improve  or 
optimize  the  unit  operations  of  a  rapid  excavation  system, 

especially  the  materials  handling  subsystem.  The  model  used 

•»* 

will  be  a  mathematical  one,  constructed  for  use  on  a  digital 
computer . 

The  process  of  modeling  is  normally  carried  out  in  a 
series  of  five  steps  generally  referred  to  as  the  scientific 
method.  The  scientific  method  of  making  decisions  is  often 
referred  to  as  the  systems  approach  and  generally  consists 
of  the  following  steps: 

1)  Definition  and  breakdown  of  the  system 

2)  Construction  of  a  model  of  the  system 


*The  numbers  in  parentheses  refer  to  the  numbered 
publications  in  the  REFERENCES  section. 
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3)  Testing  of  the  model 

4)  Solution  of  the  problem 

5)  Implementation  of  the  solution 

This  year’s  work  is  involved  primarily  with  the  first  two 
steps  above,  beginning  with  the  definition  of  the  problem 
and  continuing  through  the  construction  of  the  computer 
model  of  the  system. 

Types  of  Simulation  Models 

In  order  to  simulate  any  particular  process  or  system, 
some  type  of  a  model  is  required.  Several  types  of  models 
exist,  but  only  three  general  types  are  extensively  used. 
These  are  the  physical  models,  the  analog  models,  and  the 
conceptual  models. 

A  physical  model  is  a  physical  model  or  replica  of  a 
system,  generally  scaled  down  to  a  size  which  is  more  easily 
handled  than  the  full-size  system.  The  usual  reason  for 
using  a  physical  model  is  economy  of  operation.  The  model 
can  be  used  to  simulate  the  operation  of  the  actual  system 
without  incurring  the  cost  of  the  full-scale  system.  Phys¬ 
ical  models  are  seldom  used  in  systems  analysis  but  can 
often  be  used  in  other  fields  of  engineering  such  as  in 
aeronautical  evaluation  of  aircraft  design.  A  physical 
model  is  easy  to  "understand"  since  looks  like  the  object 
that  it  represents  or  models. 

The  second  class  of  simulation  models  are  the  analog 
variety.  An  analog  model  is  a  system,  such  as  an  electrical 
or  hydraulic  circuit,  which  can  be  constructed  to  relate 
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to  another  system  in  such  a  manner  that  the  behavior  of 
the  model  can  solve  problems  in  the  analogous  system  we 
are  interested  in.  A  typical  example  of  this  type  of  model 
is  the  electrical  network  analyzers  used  to  solve  prob¬ 
lems  related  to  mine  ventilation  circuits.  Analog  models 
are  useful  only  in  certain  types  of  problems,  but  pro¬ 
vide  rapid,  convenient  answers  in  situations  where  they 

apply. 

The  conceptual  models,  often  called  logical  or  mathe¬ 
matical  models,  are  the  prevalent  model  type  and  are  put 
to  use  on  a  wide  variety  of  problems.  For  this  type  of 
model,  the  components  of  the  system  are  represented  by 
mathematical  formulas,  probability  distributions,  or 
numerical  data  which  is  used  to  model  the  system.  A 
mathematical  model  is  normally  written  in  a  computer 
language  so  that  the  massive  chore  of  performing  the 
simulation  may  be  done  by  computer.  Most  of  the  mathemati¬ 
cal  simulation  models  fall  into  the  class  known  as  the 
Monte  Carlo  methods.  In  these  methods,  the  general  ap¬ 
proach  is  to  run  and  rerun  the  simulation  process  as  a 
statistical  experiment,  measuring  the  results  in  order  to 
learn  something  about  the  process  simulated.  The  Monte 
Carlo  methods  are  subdivided  into  the  stochastic  and  deter¬ 
ministic  models. 

Stochastic  Simulation.  Stochastic  or  probabilistic 
simulation  models  are  used  in  situations  where  the  elements 
of  the  model  are  probabilistic  or  random  in  nature,  i.e., 
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the  elements  of  the  model  cannot  be  predicted  with  cer¬ 
tainty.  A  stochastic  model  operates  with  the  probability 
distributions  of  each  element  in  the  model  and  empirically 
determines  just  what  will  happen  in  a  particular  system  by 
modeling  the  system  under  specific  sets  of  conditions.  By 
studying  the  responses  which  occur  due  to  changing  the  con¬ 
trollable  variables,  the  system  can  be  optimized.  The 
principal  advantage  of  this  class  of  model  is  that  it  may 
be  used  to  solve  many  problems  which  cannot  even  be  ap¬ 
proached  using  conventional  theoretical  methods. 


Deterministic  Simulation.  Deterministic  simulation 
has  been  described  by  Hammersley  and  Handscomb  (1)  as  an  at¬ 
tempt  to  "exploit  the  strength  of  theoretical  mathematics 
while  avoiding  its  associated  weakness  by  replacing  theory 
by  experiment  whenever  the  former  fails."  Deterministic 
simulation  is  used  to  model  processes  which  are  governed  at 
least  in  part  by  specific  laws  or  rules  and  which  will  yield 
predictable  results.  For  this  reason,  deterministic  simula¬ 
tion  has  been  used  to  simulate  such  activities  as  truck 
haulage  (8),  rail  haulage  (7),  and  the  operation  of  bucket 
wheel  excavators  (14).  In  these  applications,  physical  laws 
were  used  to  determine  accelerations,  speeds,  distances, 
power  consumptions,  etc.,  as  a  function  of  the  operating 
characteristic  curves  for  the  equipment  used.  Normal 
practice  in  a  model  of  this  sort  is  to  calculate  the 
required  variables  at  equal  intervals  of  time  in  an 
interative  fashion.  At  each  iteration,  the  theoretical 
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laws  can  be  used  to  calculate  the  desired  variables,  thus 
using  the  power  of  the  digital  computer  to  eliminate  the 
necessity  for  extensive  mathematical  development.  The 
model  can  be  used  to  study  an  activity  based  on  a  theoret¬ 
ical  basis  and  possibly  optimize  the  activity  by  inter- 
pretting  the  outcome  of  the  simulation  experiments. 

Deterministic  simulators  are  sometimes  further  sub¬ 
divided  into  event-oriented  and  time-oriented  models.  The 
time-oriented  model  is  perhaps  more  widely  used  than  the 
other  and  often  is  the  easiest  to  program.  In  this  type 
of  model,  a  specific  increment  of  time  is  chosen  previous 
to  each  computer  run.  The  program  updates  the  simulation 
by  that  time  increment  and  calculates  all  the  variables 
of  record  at  the  new  time.  The  calculations,  are  repeated 
at  each  incrementation  in  the  time  variable.  By  using  the 
proper  logic,  any  variable  can  be  accurately  determined  in 
the  simulation  if  the  concepts  for  simulating  that  variable 
are  valid. 

An  event-oriented  deterministic  simulator  is  a  simu^ 
lator  which  does  not  update  its  time  variable  by  a  constant 
value  but  instead,  updates  the  time  variable  only  when, 
specific  predetermined  events  occur  in  the  simulation.  The 
events  chosen  to  result  in  updating  are  generally  the  com¬ 
pletion  of  activities  after  which  decisions  must  be  made. 
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The  principal  advantage  of  this  method  is  that  all  the 
variables  of  record  may  not  need  updating  during  a  parti¬ 
cular  time  span.  By  using  the  knowledge  about  specific 
events  in  the  process  to  be  modeled,  only  those  variables 
of  record  which  require  updating  are  calculated  by  the 
computer.  A  disadvantage  of  the  method  is  that  it  may 
require  more  programming  work  than  the  time-oriented  model 
for  the  same  system.  The  choice  between  the  emphasis  on 
time  increments  or  events  will  depend  upon  the  system  to 
be  modeled.  It  may  not  be  obvious  which  is  the  most  advan¬ 
tageous  before  the  program  is  initiated. 

Choice  of  a  Simulation  Language 

One  of  the  first  important  tasks  involved  in  construc¬ 
ting  a  computer  model  is  the  choice  of  a  medium,  i.e.,  a 
computer  language,  in  which  to  write  the  model.  There  are 
numerous  computer  languages  to  choose  from,  including  gen¬ 
eral  languages  and  those  specifically  designed  for  appli¬ 
cation  to  simulation. 

Several  general  simulation  languages  are  available  for 
use  such  as  GPSS  '  (General  Purpose  System  Simulator)  and 
SIMSCRIPT.  These  languages  are  designed  to  handle  variations 
of  standard  simulation  problems  which  are  often  encountered. 
GPSS,  for  example,  is  best  suited  to  problems  related  to 
scheduling  or  to  systems  involving  queueing  while  SIMSCRIPT 
is  most  applicable  to  inventory  and  similar  problems.  Sev¬ 
eral  other  languages  are  available  which  are  designed  to 
study  situations  of  a  more  specific  nature.  DYNAMO  and 
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SIMULATE  are  languages  which  are  used  to  simulate  economic 
systems.  More  complete  descriptions  of  these  programs 
can  be  obtained  in  the  computer  language  manuals  and  in 
books  on  computer  simulation  (6). 

One  language  which  merits  special  attention  here  is 
GASP  II.  This  is  a  FORTRAN-based  language  which  is  widely 
applicable  and  which  has  numerous  advantages.  The  origin¬ 
ators  of  the  language  outlined  these  advantages  in  their 
manual  on  GASP  II  (9).  The  most  important  advantages  are 
related  to  GASP's  base  in  a  common  computer  language.  As 
a  result,  the  user  does  not  have  to  learn  a  new  language 
or  obtain  a  new  compiler  for  his  present  machine.  Thus, 
two  of  the  major  problems  related  to  using  a  simulation 
language  are  eliminated.  In  addition  to  these  points,  GASP 
is  a  versatile  tool  which  will  have  appeal  in  many  simula¬ 
tion  analyses. 

Another  possible  language  for  use  in  simulation  is  a 
general  purpose  language  such  as  FORTRAN.  While  this 
language  was  not  designed  for  specific  use  as  a  simulation 
language,  it  is  widely  used  as  such  and  has  several  advan¬ 
tages  as  a  simulation  language.  The  advantages  that  GASP 
II  offers  to  simulation  can  also  be  obtained  from  FORTRAN. 
Thus,  FORTRAN  is  advantageous  since  it  is  widely  understood 
and  does  not  require  a  special  compiler.  FORTRAN  does 
present  some  problems  for  simulation.  These  include  the 
lengthy  input-output  formatting  and  the  lack  of  inherent 
debugging  aids.  However,  these  disadvantages  will  not  be 
serious  ones  if  the  programmer  is  quite  familiar  with  the 
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language  and  will  not  effect  the  program  users. 

With  these  facts  in  mind,  the  choice  of  FORTRAN  was 
made  for  the  simulation  model  being  constructed.  The  main 
factors  affecting  this  decision  are  its  wide  use  and  its 
ease  of  transfer  from  one  machine  to  another.  Most  of  the 
important  simulation  work  done  in  the  mining  and  construc¬ 
tion  industry  has  been  performed  by  FORTRAN  programs  to  date. 
In  addition,  nearly  every  digital  computer  has  FORTRAN 
capability  and  this  will  enable  the  model  to  be  used  on  the 
maximum  number  of  computers.  To  further  minimize  transfer 
problems,  the  authors  of  the  model  have  attempted  to  follow 
USA  Standard  FORTRAN  IV  as  published  by  the  United  States 
of  America  Standards  Institute  (13).  This  will  minimize 
the  machine-dependent  statements  which  will  require  changing 
when  the  program  is  used  on  other  machines. 


description  of  the  model 
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The  computer  model  presented  here  Is  a  Monte  Carlo  type 
model  written  in  the  FORTRAN  language  using  both  determin¬ 
istic  and  stochastic  simulation  methods  to  model  the  over¬ 
all  tunneling  system.  The  program  is  written  in  an  event- 
oriented  manner  with  program  updating  being  accomplished 
after  specific  jobs  or  events  are  completed.  Most  of  the 
unit  operation  submodels  are  written  in  stochastic  form 
although  the  materials  handling  subsystem  contains  much  in 
the  way  of  deterministic  calculations.  Emphasis  has  been 
placed  upon  supplying  a  number  of  options  within  the  program 
to  make  the  program  applicable  to  various  types  or  forms  of 
rapid  excavation  systems.  This  portion  of  the  report  deals 
with  the  model  objectives,  the  description  of  the  simulation 
concepts,  the  logic  used  and  the  outline  of  program  organi- 
zat ion. 

Model  Objectives 

The  primary  goal  of  this  model  is  to  simulate  the  common 
methods  of  driving  a  tunnel  with  a  boring  machine.  To  accom¬ 
plish  this  goal,  it  is  necessary  to  think  in  terms  of  a 
general  computer  program  which  contains  a  number  of  options 
which  allow  a  user  to  vary  the  simulation  of  the  unit  oper¬ 
ations  and  the  way  that  they  interact  during  the  tunneling 
process.  Primary  attention  is  paid  in  this  model  to  the 
materials  handling  process  as  this  is  one  unit  operation 
which  promises  to  yield  results  from  a  systems  evaluation. 

This  conclusion  is  based  upon  observations  about  the 
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materials  handling  function  creating  a  bottleneck  in  the 
operation  (2,3,4,10)  and  due  to  the  fact  that  more  control 
may  be  exercised  over  the  design  and  operation  of  the  mater¬ 
ials  handling  process  than  over  the  other  unit  operations. 

For  this  reason,  the  most  significant  programming  time  and 
attention  was  devoted  to  the  modeling  of  the  materials 
handling  function. 

To  meet  the  basic  objective  of  studying  primarily  the 
materials  handling  process,  models  for  both  cyclic  and 
continuous  handling  methods  have  been  provided  so  that  either 
type  may  be  studied.  The  cyclic  systems  have  been  pro¬ 
grammed  in  a  fashion  which  will  allow  either  a  track  or  a 
rubber-tired  haulage  system  to  be  modeled  providing  that 
the  characteristic  curves  of  the  driving  mechanism  are 
available.  For  continuous  systems,  similar  accommodations 
have  been  provided  so  that  either  belt  or  hydraulic  con¬ 
veyors  may  be  simulated. 

Outline  of  Simulation  Concepts  Used 

The  outline  of  the  logic  and  concepts  used  in  the 
simulation  model  will  deal  first  with  the  general  principles 
or  concepts  used  throughout  the  program.  Afterwards,  those 
concepts  which  apply  primarily  to  the  individual  operations 
will  be  discussed.  For  purposes  of  outlining  these  specific 
concepts,  the  tunneling  process  will  be  divided  into  the 
following  unit  operations; 

1)  muck  generation 

2)  materials  handling 


3)  roof  support 

4)  environmental  control 

Each  of  these  unit  operations  will  be  discussed  separately 
even  though  they  may  not  be  programmed  in  separate  units 
in  the  program  itself. 

General  Concepts.  The  first  of  the  discussions  on  gen¬ 
eral  concepts  should  perhaps  be  centered  around  the  method 
of  introducing  the  necessary  probability  functions  into  the 
program.  For  versatility  and  ease  of  input,  all  the 
probability  functions  which  are  used  in  the  program  are 
introduced  as  piecewise  linear  cumulative  probability 
functions  which  are  sometimes  also  referred  to  as  cumulative 
frequency  polygons  or  ogives  (11).  Figure  1  illustrates 
the  method  for  reading  the  cumulative  probability  functions 
into  the  program.  Several  things  should  be  mentioned  here 
regarding  these  functions: 

1)  Neither  the  abscissa  nor  the  ordinate  values 
must  be  evenly  spaced. 

2)  The  first  ordinate  value,  shown  in  Figure  1  as 
CP  (l)s  must  equal  zero. 

3)  The  final  ordinate  value,  shown  in  Figure  1  as 
CP(NPOINT),  must  equal  one. 

4)  The  ordinate  and  abscissa  values  are  read  into  the 
program  as  pairs  and  must  be  arranged  in  terms 

of  increasing  ordinate  or  cumulative  probability 


values  . 
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5)  The  numoer  of  abscissa  and  ordinate  values  read  in 
may  be  up  to  13.  If  more  are  necessary,,  the  dimen¬ 
sions  of  the  necessary  variables  may  be  easily 
changed  to  provide  the  additional  storage  space. 
Should  the  user  decide  that  a  constant  value  is  to  be  read 
into  the  program  for  a  particular  variable  instead  of  a 
distribution  of  values,  he  may  do  so  under  the  framework 
of  the  above  method.  The  procedure  that  should  be  used  is 
to  read  in  two  ordinate  and  abscissa  values;  the  first 
ordinate  value  should  be  zero  and  the  second  should  be  one 
while  both  abscissa  values  should  be  equal  to  the  constant 
desired  for  that  variable.  For  example,  if  the  user  wished 
to  read  in  a  constant  value  of  10.5  for  a  specific  variable, 
he  would  read  in  the  following  values  for  the  cumulative 
frequency  polygons 

CP(1)  =  0.0  T ( 1 )  =  10.5 

CP  ( 2 )  =  1.0  T ( 2 )  =  10.5 

The  computer  would  then  automatically  assign  a  value  of  10.5 

to  the  variable  in  question  every  time  it  is  called  in  the 
program . 

Since  the  computer  program  described  here  is  classi¬ 
fied  as  an  event-o r lent ed  model  rather  than  a  time-oriented 
model,  a  simplified  explanation  of  an  event-oriented  model 
is  presented  to  provide  a  basic  background  for  users.  The 
computer  will  store  the  clock  time  for  all  pertinent  events 
in  storage.  In  searching  for  the  activity  which  should  be 
updated  next,  the  computer  will  go  to  the  activity  with  the 
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shortest  clock  time.  An  example  of  how  this  would  work  is 
illustrated  in  Figure  2.  In  Figure  2(a),  the  status  of  the 
five  activities  assumed  to  exist  in  the  problem  are  shown 
in  Gantt  chart  fashion.  The  tick  marks  shown  are  indicative 
of  specific  events  such  as  the  completion  of  certain  jobs 
or  tasks.  Since  Activity  2  exhibits  the  shortest  clock 
times  the  computer  must  deal  with  or  act  upon  Activity  2 
before  it  proceeds  to  the  activity  with  the  next  shortest 
clock  time.  If  it  is  possible  to  update  Activity  2  beyond 
its  present  clock  time,  then  this  is  done  as  shown  in 
Figure  2(b)  and  the  computer  then  focuses  attention  on  the 
new  activity  which  has  the  shortest  clock  time,  Activity  3. 

If  the  situation  requires  that  the  updating  of  Activity 

2  is  restrained  by  another  activity,  then  the  start  of 
another  cycle  of  Activity  2  may  not  begin  immediately.  A 
very  simple  example  is  presented  in  Figure  2(c)  where  the 
start  of  Activity  2  is  assumed  to  be  restrained  by  Activity 

3  (and  only  Activity  3).  This  type  of  situation  may  arise 
because  of  manpower,  space,  sequencing,  or  other  limitations 
In  any  case,  in  this  situation  Activity  2  must  wait  until 
Activity  3  is  completed  before  it  can  be  reinitiated.  Thus, 
the  wait  is  indicated  by  a  dotted  rather  than  a  solid  line. 
After  Activity  3  is  completed,  Activity  2  is  simulated  to 
completion  and  Activity  3  then  has  the  shortest  clock  time 
and  is  considered  for  updating  next.  In  reality,  the 
simulation  of  an  activity  may  be  restrained  by  numerous 
other  activities  of  different  types.  However,  the  general 
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(b)  Normal  Updating  of  Activity  2 


(c)  Updating  of  Activity  2  After 
Completion  of  Activity  3 


Figure  2 


Simple  Examples  of  Updating  in 
Oriented  Simulation  Model. 


CLOCK  TIME 


CLOCK  TIME 


CLOCK  TIME 


an  Event- 


principle  of  focusing  on  the  activity  with  the  shortest 
clock  time  will  apply  no  matter  how  complex  the  logic,  pro¬ 
viding  that  the  model  has  been  properly  programmed. 

The  assignment  and  utilization  of  manpower  is  another 
consideration  which  applies  throughout  the  computer  model. 
The  method  of  allocation  of  manpower  was  aimed  at  maximum 
versatility  in  the  number  of  men  assigned  to  a  particular 
job.  For  each  task  in  the  tunneling  process,  an  upper  and 
lower  limit  on  the  number  of  workmen  assigned  is  read  into 
the  program.  The  lower  limit  will  reflect  the  minimum  num¬ 
ber  of  men  required  to  safely  carry  out  a  task.  The  upper 
limit  will  generally  be  determined  by  space,  productivity, 
safety,  or  other  practical  limitations  of  the  activity.  The 
computer  program  will  always  assign  at  least  the  minimum 
number  of  men  to  a  job  before  it  is  initiated  and  will  assign 
as  many  men  as  it  can  subject  to  availability  and  upper 
limit  restrictions.  As  more  men  are  assigned  to  a  job,  the 
time  to  accomplish  the  job  is  reduced  proportionally.  This 
policy  is  based  upon  the  assumption  that  the  upper  and  lower 
limits  of  manpower  are  reasonable  and  that  all  men  are  gain¬ 
fully  occupied  on  any  particular  job.  As  each  job  is  com¬ 
pleted,  the  men  assigned  to  that  job  are  reassigned  to  other 
jobs  if  it  is  possible.  When  several  jobs  require  manpower 
simultaneously,  the  largest  job  in  terms  of  manhours  required 
is  assigned  men  first. 

On  additional  general  topic  of  discussion  here  is  the 
options  available  for  outputting  information  from  the  com¬ 
puter.  At  the  termination  of  each  simulation  run,  a  listing 
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of  summary  statistics  is  printed  routinely.  In  order  to 
allow  users  to  determine  jus^  what  is  taking  place  in  the 
computer  program,  a  log  of  operations  which  outputs  informa¬ 
tion  on  each  significant  event  in  the  simulation  as  it  occurs 
can  also  be  optionally  implemented.  If  the  log  of  opera¬ 
tions  is  not  desirable  or  necessary,  the  user  may  suppress 
this  series  of  output  statements  and  the  computer  will  print 
only  the  simulation  summary  statistics. 

Muck  Generation  Subsystem.  The  muck  generation  subsys¬ 
tem  includes  all  the  activities  taking  place  at  the  face  of 
the  tunnel  concerned  with  the  operation  of  the  tunneling 
device.  Thus,  the  muck  generation  subsystem  is  concerned 
primarily  with  the  rate  of  advance,  the  inspection  and  re¬ 
pair  and  replacement  of  bits,  and  the  repair  and  maintenance 
of  the  tunneling  device.  The  bits  are  one  of  the  most 
important  of  the  considerations  in  the  generation  of  muck, 
particularly  in  large  tunnels  driven  in  hard  rock.  Each 
bit  on  the  face  of  the  mole  must  be  numbered  for  the  pur¬ 
poses  of  the  computer  program.  This  can  be  done  as  shown  in 
Figure  2  of  a  previous  report  (5)  or  in  any  other  suitable 
manner.  After  numbering  each  bit,  a  time-to-f ailure  pro¬ 
bability  function  is  assigned  to  each  bit  with  the  proba¬ 
bility  being  expressed  in  terms  of  the  feet  of  advance.  In 
addition,  another  distribution  for  the  replacement  manhours 
required  is  assigned  to  each  bit  location  in  order  to  dif- 
f grenti-t ts  between  bits  in  terms  of  the  replacement  time 
required.  A  separate  time-to-f ailure  and  repair  time 
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distribution  is  provided  to  similate  repair  work  on  the  bits 
which  does  not  require  replacement,  e.g.9  welding  or  other 
repair  work  on  the  bit  housings.  When  a  bit  reaches  the 
point  of  failure,  it  is  not  replaced  immediately  but  is 
replaced  at  the  first  inspection  after  the  failure  has 
occurred  9  i.e.9  at  the  first  opportunity  for  the  failure 
to  be  discovered. 

The  inspection  of  the  bits  are  assumed  in  the  program 
to  be  completed  in  conjunction  with  the  resetting  of  the 
jacks  after  completion  of  a  normal  stroke  or  on  any  occasion 
in  which  the  machine  is  down  for  other  purposes.  It  is 
assumed  to  be  made  normally  after  any  integer  number  of 
cycles 9  i.e.9  after  the  jacks  have  been  reset  a  predetermined 
number  of  times.  If  the  bits  are  found  to  be  in  condition 
for  more  boring9  the  boring  is  reinitiated.  If  failed  or 
vjorn  bits  are  detected,  the  replacement,  operation  is  simu¬ 
lated  before  the  boring  is  continued.  Some  tolerance,  input¬ 
ted  in  terms  of  feet  of  advance,  is  allowed  in  the  program 
so  that  worn  bits  do  not  have  to  be  replaced  the  instant 
their  generated  lifetime  is  assumed  to  end. 

The  muck  generation  subsystem  also  includes  provision 
for  repairs  and  maintenance  which  must  be  performed  on  the 
tunneling  machine.  Those  repairs  which  result  in  the  shut¬ 
down  of  the  system  are  compiled  into  a  time-to-f ailure  dis¬ 
tribution.  A  distribution  of  manhours  required  for  these 
repairs  is  also  provided  to  complete  the  simulation  of  this 
part  of  the  process.  In  all  cases  of  simulating  repairs 
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associated  with  the  mole,  the  tunneling  machine  is  assumed 
to  be  down  in  the  model  and  as  many  crewmen  as  possible 
under  the  circumstances  are  assigned  to  the  repair  action 
in  order  to  expedite  the  boring  operation.  All  of  the  above 
processes  are  simulated  in  a  relatively  straightforward 
stochastic  manner.  This  is  accomplished  by  placing  each 
event  (bit  failure,  mole  failure,  etc.)  in  an  event  matrix 
and  testing  at  each  update  time  to  see  if  any  action  is 
required.  In  this  manner,  all  events  in  the  muck  genera¬ 
tion  subsystem  are  handled  in  the  same  matrix  and  are  scanned 
at  the  same  time  in  the  program. 

The  final  important  element  In  this  subsystem  is  the 
rate  of  generation  of  muck  during  the  operation  of  the  mole. 
This  process  is  accomplished  in  the  program  through  the 
advance  rate  distribution  and  the  geometry  of  the  face.  The 
advance  rate  potential  of  the  tunneling  device  in  feet  per 
hour  can  be  formed  into  a  probability  distribution.  A  ran¬ 
dom  sample  from  this  distribution  is  chosen  to  obtain  an 
advance  rate  which  applies  for  the  advance  of  one  stroke 
length  of  the  machine.  This  advance  rate  is  then  combined 
with  the  tunnel  cross-section  to  determine  the  muck  flow 
rate.  An  instantaneous  advance  rate  would  have  been  more 
precise  but  the  result  in  terms  of  the  simulation  would  have 
been  negligible,  i.e.,  the  long-term  production  of  the  machine 
does  not  appear  to  be  sensitive  to  this  variable.  In  the 
computer  program,  the  simulation  of  the  muck  generation  sub¬ 
system  is  carried  out  in  the  main  program  and  in  SUBROUTINE 
MUCK. 
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Materials  Handling  Subsystem.  The  materials  handling 
subsystem  was  the  most  complex  portion  of  the  overall  model 
to  program.  This  situation  existed  as  a  result  of  the  em¬ 
phasis  placed  upon  the  materials  handling  process  in  the 
model  and  the  physical  complexity  of  some  of  the  muck  hand¬ 
ling  systems.  Simulation  of  both  cyclic  and  continuous  sys¬ 
tems  have  been  provided  for  in  the  model.  The  computer  model 
is  designed  in  such  a  manner  that  the  haulage  distance  is 
increased  as  the  tunnel  is  advanced.  This  is  accomplished 
by  keeping  track  of  the  advance  and  increasing  the  haulage 
length  each  time  a  predetermined  advance,  DELTH,  is  attained. 
This  also  results  in  changes  in  the  inby  end  of  the  haulage 
system  which  must  be  reflected  within  the  model. 

(1)  Cyclic  Systems.  Materials  handling  using  cy¬ 
clic  systems  are  the  most  complex  methods  from  a  systems 
standpoint.  The  cyclic  materials  handling  systems  were 
modeled  primarily  with  single-track  haulage  systems  in  mind 
but  a  haulage  system  using  rubber-tired  vehicles  can  be 
accommodated  using  the  same  model  since  the  simulation  pro¬ 
gram  is  designed  with  this  in  mind.  The  initial  concern  of 
the  cyclic  materials  handling  model  to  be  discussed  here  is 
the  method  of  introducing  the  tunnel  grade  characteristics 
into  the  program.  This  is  accomplished  by  dividing  the 
tunnel  into  sections  with  each  section  having  a  constant 
grade.  In  case  of  a  tunnel  with  continuously  varying  grade, 
the  tunnel  profile  may  have  to  be  approximated  by  the  assumed 
linear  grade  segments.  The  segments  are  read  into  the  program 
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in  order  proceeding  from  the  dumping  point  and  continuing 
to  the  face  of  the  tunnel  as  shown  in  Figure  3  where  a  tunnel 
profile  with  five  sections  is  illustrated.  In  all  cases, 
the  distances  are  measured  along  the  center  line  of  the 
tunnel  and  changes  in  azimuth  are  ignored  and  assumed  to  be 
of  little  or  no  consequence  in  the  movement  of  the  haulage 
devices  as  they  traverse  the  tunnel.  For  programming  reasons, 
the  sections  outside  the  portal  are  counted  separately  from 
the  sections  within  the  tunnel.  The  program  will  accommodate 
a  tunnel  profile  with  100  sections  without  alteration. 

The  switches,  or  switchpoints  in  the  case  of  rubber- 
tired  vehicles,  are  assumed  to  be  evenly  spaced  along  the 
tunnel  route.  For  rubber-tired  vehicles,  a  bored  tunnel  is 
not  an  ideal  roadbed  and  thus  it  is  not  usually  possible  for 
the  vehicles  to  pass  anywhere  except  vhere  special  passing 
points  have  been  blasted  out  of  the  tunnel.  For  this  reason, 
the  simulation  model  is  assumed  to  be  able  to  model  this 
type  of  haulage  system  with  passing  points  at  equal  inter¬ 
vals  along  the  tunnel.  The  cyclic  materials  handling  sub¬ 
model  simulates  the  movement  of  the  vehicles  on  a  switch-to- 
switch  basis  in  SUBROUTINE  TRANS.  For  example,  assume  that 
a  train  is  waiting  on  the  inbound  side  of  Switch  B  of 
Figure  4  on  one  of  its  empty  trips  to  the  face  of  the  tunnel. 
When  the  track  is  cleared,  SUBROUTINE  TRANS  controls  the 
movement  of  the  empty  train  by  calling  SUBROUTINE  MOTION 
which  simulates  the  motion  of  the  train  from  Switch  B  to 

In  order  to  obtain  clearance  to  use  the  section 


Switch  A. 
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Figure  4  -  Diagram  of  Two  Adjacent  Switches 
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of  track  between  the  two  switches,  the  track  section  must 
be  clear  and  the  empty  train  must  have  priority  as  deter¬ 
mined  by  the  decision  or  control  function  in  SUBROUTINE  TRANS. 
The  decision  as  to  which  train  has  priority  to  a  particular 
section  of  track  is  made  on  a  f ir st-in-f ir s t- out  basis. 
Adjacent  switches,  such  as  Switch  A  and  Switch  B„  are  al¬ 
ways  considered  together  in  determining  this  priority.  For 
example,,  if  an  inbound  train  reaches  Switch  B:  before  an 
outbound  train  reaches  Switch  A,  then  the  inbound  train  has 
the  priority  for  the  use  of  the  connecting  tiack  and  it  com¬ 
pletes  its  movement  to  Switch  A  before  the  olitbound  train 
can  initiate  its  move  from  Switch  A  to  Switch  B.  By  con¬ 
sidering  all  the  switches  simultaneously,  SUBROUTINE  TRANS 
can  control  the  operation  of  all  the  trains  in  an  event- 
oriented  fashion  while  SUBROUTINE  MOTION  simulates  the  ac¬ 
tual  swi tch- to~switch  movements. 

SUBROUTINE  MOTION  handles  the  motion  of  the  train  in 
an  event-oriented  deterministic  fashion  based  upon  the  physi¬ 
cal  laws  of  motion.  One  of  the  first  publications  dealing 
with  this  basic  simulation  method  for  haulage  systems  was 
introduced  by  Nelson  (7).  For  this  application,  his  basic 
deterministic  approach  has  been  changed  to  one  which  does 
not  make  use  of  equal  time  increments  but  instead  concen¬ 
trates  upon  specific  events  in  the  movement  of  the  train  as 
its  travel  is  simulated.  The  basic  physical  law  used  is 
Newton's  second  law  of  motion  which  for  the  case  of  a  rolling 
vehicle  (12)  can  be  written  ass 
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(T  -  Ff  -  Fg)G 

a  =  w„  +  w  +  w 

il  c  m 

where;  T  =  tractive  effort  of  the  driving  wheels  in  pounds 

F^  =  force  required  to  overcome  friction  in  pounds 

F  =  force  required  to  overcome  the  gravity  com- 

6 

ponent  in  pounds 

Wj  =  weight  of  the  locomotive  in  pounds 
''  =  weight  of  the  cars  in  pounds 

W  =  weight  of  the  muck  in  pounds 

a  =  acceleration  in  feet  per  second  per  second 
g  =  acceleration  of  gravity,  32.2  feet  per  second 
per  second 

Since  the  tractive  effort  does  not  remain  constant  for 
changes  in  the  speed  of  the  tractive  unit,  some  method  of 
applying  the  formula  above  must  be  used  so  that  the  changes 
in  the  speed  and  the  tractive  effort  are  reflected  in  the 
program.  To  accomplish  this,  the  characteristic  curve 
of  the  tractive  unit  which  relates  its  speed  and  tractive 
effort  must  be  made  available  for  use  in  the  computer  model. 
A  number  of  selected  points  along  this  characteristic 
curve  are  read  into  the  computer  program  as  shown  in  Figure 
5  for  a  hypothetical  two-speed  locomotive  unit.  The  program 
then  assumes  that  the  characteristic  curve  is  linear  between 
succeeding  points  so  that  the  effect  is  an  approximation 
of  the  actual  curve  by  a  piecewise  linear  function  defined 
by  the  points  selected  for  input.  The  degree  of  simulation 
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accuracy  required  in  the  runs  will  dictate  the  number  and 
spacing  of  the  points  selected.  At  present,  the  proper 
variables  in  the  computer  program  are  dimensioned  to  allow 
reading  in  up  to  30  points  along  this  characteristic  curve. 

The  use  of  the  tractive  effort-speed  curve  in  SUB¬ 
ROUTINE  MOTION  is  carried  out  on  an  iterative  basis  using 
certain  specified  events  to  indicate  the  need  for  recal¬ 
culation  of  the  variables  of  motion.  Normally  this  is  done 
based  upon  the  assumed  linear  segments  of  the  characteris¬ 
tic  curve  as  follows.  A  train  (or  other  vehicle)  which  is 
starting  from  rest  is  assumed  to  do  so  at  the  average 
tractive  effort  value  for  the  first  assumed  linear  segment 
along  the  curve  in  Figure  5,  i.e.,  at  a  tractive  effort 
value  of  [TE  (1)+TE  (2) ] /2 .  An  acceleration  is  calculated 
based  upon  this  tractive  effort  and  the  train  moves  until 
the  acceleration  results  in  the  train  achieving  the  speed 
at  the  end  of  the  first  linear  segment,  S(2).  When  this 
occurs,  a  new  average  tractive  effort  value,  [ TE (2)+TE (3) ] / 2  , 
is  applied  for  the  period  of  time  required  for  the  train’s 
speed  to  reach  S(3),  and  so  on.  This  iterative  method  con¬ 
tinues  until  the  train  reaches  its  maximum  allowable  speed 
or  until  it  reaches  a  new  grade  section  in  the  tunnel.  At 
the  maximum  speed,  the  train's  speed  is  not  permitted  to 
accelerate  any  further  and  it  continues  with  a  constant 
velocity.  When  a  change  in  grade  occurs,  this  changes  the 
gravity  force  component  and  thus  the  acceleration  is  auto¬ 
matically  recalculated  within  the  program  even  though  the 
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Figure  5  -  Method  of  Inputting  Data  From  the  Characteristic 
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train  has  not  speeded  up  to  the  next  input  point  on  the 
tractive  effort-speed  curve.  To  perform  this  calculation, 
the  computer  will  interpolate  to  determine  the  current 
value  of  the  tractive  effort  and  average  this  value  with 
the  next  higher  tractive  effort  value  read  in  along  the 
curve.  This  average  will  then  be  used  to  calculate  the 
initial  acceleration  on  the  new  grade.  It  is  assumed  in 
this  method  that  the  mass  of  the  train  is  a  point  mass 
located  at  the  locomotive  unit.  This  assumption  will  not 
effect  the  simulation  significantly  unless  the  tunnel  pro¬ 
file  is  changing  rapidly  and  considerably  in  grade,  a 
situation  which  does  not  occur  in  rapid  excavation  tunnel¬ 
ing  jobs. 

The  dumping,  loading,  and  switchout  times  for  the 
cyclic  materials  handling  systems  are  handled  separately. 
The  loading  times  are  determined  by  the  interaction  of 
the  muck  generation  and  the  materials  handling  systems. 

The  cars  of  a  train  are  loaded  by  the  action  of  the  mole 
as  it  advances  into  the  face.  Thus,  the  loading  time 
for  each  car  is  a  stochastic  function  which  is  dependent 
on  the  rate  of  advance  which  is  generated  in  the  program 
for  the  tunneling  device.  The  dumping  time  of  each  train 
is  also  determined  stochastically  to  allow  for  the  varia¬ 
tions  which  will  certainly  occur  in  the  process.  Thus, 
a  dumping  time  cumulative  probability  function  must  be 
read  into  the  computer  as  illustrated  in  Figure  1.  The 
switchout  time  mentioned  above  is  the  name  given  here  to 


the  time  required  for  an  empty  and  a  loaded  train  to 
switch  out  under  the  gantry  conveyor  using  the  switch 
normally  located  directly  behind  the  conveyor.  This 
process  may  be  deterministically  simulated  under  ideal 
circumstances.  However,  operators  often  use  incoming 
trips  to  haul  the  tunnel  supplies  and  these  must  be  un¬ 
loaded  when  the  train  reaches  the  face  area.  Thus,  it  is 
necessary  to  use  a  probabilistic  approach  on  the  switch- 
out  time  in  order  to  reflect  the  variations  in  time  due 
to  the  necessity  of  unloading  the  supplies  at  the  face. 
This  can  be  done  by  utilizing  a  bimodal  distribution,  the 
first  or  shortest  mode  reflecting  switchout  times  where  no 
supplies  are  unloaded  and  the  second  mode  related  to  times 
necessary  to  complete  the  switchout  operation  when  the  un¬ 
loading  time  is  included  in  the  switchout  time.  When  un¬ 
loading  of  supplies  is  not  a  problem,  a  unimodal  distri¬ 
bution  may  be  suitable  for  this  variable. 

At  the  start  of  a  simulation  run,  the  trains  are  po¬ 
sitioned  behind  the  tunneling  device  in  such  a  manner  that 
they  are  spaced  one  switch  apart.  This  setup  places  the 
trains  in  as  favorable  a  state  of  readiness  as  can  be 
achieved  in  the  tunnel.  This  initial  setup  scheme  was 
chosen  since  it  was  felt  that  the  trains  would  be  in  a 
ready  state  during  a  normal  startup  of  a  tunneling  opera¬ 
tion,  e.g.,  at  the  beginning  of  the  first  shift  of  the 


simulation. 


(2)  gP.ntlnuous  Svalna.  The  simulation  of 
a  continuous  materials  handling  system  is  simple  in  com 
parison  with  the  cyclic  systems.  To  model  the  actual 


transport  of  the  muck,  the  concept  of  effective  cross- 
aection  is  defined  as  the  area  occupied  by  the  broken 
muck  in  a  cross-section  of  the  material  flow  when  the 
materials  handling  method  is  operating  at  its  maximum 
capacity.  The  value  would  be  a  constant  for  any  system 
and  would  be  independent  of  the  flow  velocity  and  the 
material  density.  The  effective  cross-section  for  a 
belt  conveyor  and  those  for  a  hydraulic  or  pneumatic 
system  would  differ  as  shown  in  Figure  6.  In  all  cases, 
however,  when  the  effective  cross-section  is  multiplied 
by  the  velocity  of  transport  and  the  proper  density  value, 
the  result  should  be  the  maximum  mass  flow  rate  of  the 
muck  for  the  specific  materials  handling  system  used. 

Care  should  be  taken  in  expressing  the  value  of  the  density 


as  the  effective 


cross-section  of  the  belt  is  based  upon 


the  profile  of  broken  rock  while  those  for  the  systems 

using  pipe  are  based  upon  solid  material.  Once  the  muck 

has  entered  the  flowstream,  the  actual  transport  can  be 

easily  simulated.  This  can  be  modeled  deterministically 

based  upon  the  flow  velocity  and  the  length  of  the  haulage 
sys  tern . 

One  of  the  most  important  considerations  in  the 


materials  handling  subsystem  for 


continuous  systems  is 


Figure 
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THE  EFFECTIVE  CROSS-SECTION  OF  A  CONVEYOR 
IS  THE  CROSS-SECTIONAL  AREA  OF  THE  MUCK 
WHEN  THE  CONVEYOR  IS  OPERATING  AT  ITS 
MAXIMUM  CAPACITY  (INDICATED  BY  THE  CROSS- 
HATCHED  AREA  ABOVE) 


THE  EFFECTIVE  CROSS-SECTION  OF  A  HYDRAULIC 
OR  PNEUMATIC  CONVEYOR  IS  THE  CROSS- 
SECTIONAL  AREA  OCCUPIED  BY  THE  MUCK  WHEN 
THE  CONVEYOR  IS  OPERATING  AT  ITS  MAXIMUM 
CAPACITY  (INDICATED  BY  THE  SHADED  AREA 
ABOVE) 


Illustration  of  the  Effective  Cross-Section 
for  Continuous  Materials  Handling  Systems 
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the  interaction  with  the  muck  generation  subsystem.  The 
primary  function  is  to  regulate  the  flow  of  muck  into  the 
materials  handling  subsystem.  If  the  muck  generation 
rate  is  greater  than  that  which  can  be  handled  by  the 
materials  handling  subsystem,  then  the  rate  of  muck  genera¬ 
tion  is  slowed  to  permit  the  materials  handling  subsystem 
to  accommodate  the  muck.  This  would,  of  course,  slow  the 
advance  of  the  overall  system.  When  the  materials  handling 
device  can  handle  the  flow  of  muck,  the  muck  generation 
subsystem  can  then  be  allowed  to  operate  in  an  unconstrained 
manner . 

At  the  other  end  of  the  materials  handling  subsystem 
where  the  muck  is  dumped,  another  possibility  for  inter¬ 
ruptions  in  the  flow  of  -ck  occurs.  This  can  arise  be¬ 
cause  of  an  interaction  with  another  transport  system, 
beca»ue  of  the  condition  of  a  holding  device,  or  due  to 
numerous  other  factors  which  can  effect  the  flow  of  material 
from  the  tunnel.  Because  of  the  varied  nature  of  the 
possibilities  which  may  be  encountered  on  this  end,  no 
specific  delay  has  been  programmed.  However,  if  a  specific 
type  of  delay  is  expected  to  occur  at  the  discharge  point 
of  the  continuous  materials  handling  system,  this  can  be 
added  to  the  model  in  the  manner  which  will  correctly 
affect  the  simulation  of  this  characteristic  of  the  system. 

Tunnel  Support  Subsystem.  The  support  function  for 


tunneling  is  quite  variable  because  of  the  nature  of  the 
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geologic  materials  through  which  tunnels  are  driven.  Many 
excavations  are  provided  with  support  in  a  fashion  which  may 
be  considered  to  be  cyclic,  i.e.,  a  cycle  of  jobs  is  carried 
out  to  advance  the  support  by  one  "set."  In  the  computer 
program,  the  simulation  of  such  a  cyclic  method  is  carried 
out  by  assigning  a  probability  distribution  to  the  number 
o f  manhour s  required  to  advance  the  support  through  a  single 
cycle  of  support  work.  This  makes  the  interrelationship  be¬ 
tween  the  muck  generation  and  the  support  subsystems  an  easy 
one  to  handle  in  the  model.  The  time  to  advance  the  tunnel¬ 
ing  device  the  length  of  one  set  can  be  compared  to  the  time 
required  to  complete  one  cycle  of  support  and  tunnel  advance 
can  be  limited  to  the  speed  of  the  slower  process.  This 
procedure  will  permit  the  support  subsystem  to  keep  up  and 
provide  the  support  which  is  required  to  safely  advance  the 
tunnel . 

Other  methods  of  providing  support  in  a  tunnel  are  much 
less  cyclic  in  nature  and  vary  significantly  from  the  methods 
suggested  above.  Examples  of  this  type  of  support  methods 
include  roofbolting  and  guniting.  For  methods  which  are  not 
cyclic  in  nature,  the  simulation  must  be  handled  differently. 
This  can  be  done,  however,  within  the  framework  of  the  cyclic 
support  methods  outlined  above  by  shortening  the  length  of 
a  "set"  to  a  value  which  is  short  compared  to  the  stroke  of 
the  tunneling  machine.  In  this  manner  the  simulation  will 
approach  the  installation  of  support  which  occurs  contin¬ 
uously  rather  than  one  which  causes  the  support  to  be  advanced 
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in  spurts.  As  an  illustration,  the  action  of  installing 
roof  bolts  may  be  modeled  by  inputting  the  probabilistic 
number  of  manhours  required  to  advance  the  support  a 
relatively  short  distance  along  the  tunnel  e.g.,  one  foot. 
As  the  support  is  advanced,  the  advance  of  the  mole  can 
be  checked  to  insure  that  it  does  not  exceed  the  advance  of 
the  available  support  exactly  as  was  done  for  the  cyclic 
systems.  Since  the  support  is  not  advanced  in  long  incre¬ 
ments,  however,  the  model  is  realistic  in  relation  to  the 
actual  system. 


Environmental  Control  Subsystem.  The  primary  tasks  in 
providing  an  adequate  environment  throughout  the  tunnel 
normally  involve  extension  of  the  ventilation  system  and 
maintaining  a  water  supply  if  used  on  the  cutting  head  to 
aid  in  dust  abatement.  The  process  of  supplying  these 
auxiliary  needs  will  normally  be  performed  at  specific  inter¬ 
vals  of  tunnel  advance.  The  installation  of  the  ventilation 
tubing  is  normally  undertaken  at  intervals  of  advance  equal 
to  the  length  of  the  tubing  sections.  The  simulation  of 
the  installation  is  performed  stochastically  by  providing 
a  probability  function  for  the  number  of  manhours  required 
to  install  one  length  of  the  ventilation  tubing.  Provision 
has  been  made  for  allowing  the  tunnel  to  advance  by  more  than 
one  length  of  the  tubing  before  the  installation  of  the  tub¬ 
ing  must  be  undertaken.  A  similar  method  is  applied  to  the 
process  of  maintaining  the  supply  of  water  at  the  face.  A 
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separate  probability  distribution  for  the  demands  of  this 
system  is  read  into  the  computer  for  each  run. 

Additional  auxiliary  services  may  be  necessary  at  the 
face  which  may  or  may  not  be  directly  related  to  the  environ¬ 
mental  control  function.  These  may  include  such  functions 
as  the  advancing  of  the  track,  the  extension  of  the  sump 
lines,  or  other  jobs  which  must  be  carried  out  on  a  periodic 
basis.  These  processes  may  be  simulated  within  the  environ¬ 
mental  control  subsystem  just  as  those  functions  directly 
connected  to  the  environment  in  the  tunnel.  A  third  periodic 
process  of  this  type  can  be  simulated  by  using  the  probabil¬ 
ity  distribution  already  provided  within  this  subsystem. 

Other  functions  of  ;•  similar  nature  can  be  handled  if  neces¬ 
sary  by  providing  additional  distributions  and  using  the 
framework  of  logic  inherent  in  the  environmental  control 


subsystem . 
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TESTING  OF  THE  MODEL 

The  testing  of  the  computer  model  was  only  partially 
completed  at  the  end  of  the  project  year.  The  initial 
testing  phase  concerned  with  checking  the  logic  of  the  pro¬ 
gram  and  its  macro  behavior  was  accomplished  using  data 
obtained  mainly  in  the  field.  Howevers  more  exhaustive  eval¬ 
uation  and  development  was  scheduled  for  the  second  year  of 
the  project  and  is  yet  to  be  undertaken. 

Data  Collection 

In  the  testing  of  the  computer  program,  as  much  data 
as  possible  from  the  field  was  used  to  supply  the  computer 
program.  In  the  muck  generation  subsystem,  data  obtained 
through  the  courtesy  of  the  White  Pine  Copper  Company  was 
used  in  the  simulation.  The  bit  life  distributions  were 
compiled  from  actual  bit  records  kept  by  the  mine  personnel 
during  the  period  of  experience  with  their  Robbins  machine. 
The  bit  lives  available  were  formed  into  a  histogram  for 
each  bit  on  the  head  of  the  machine.  The  histograms  were 
formed  from  the  raw  data  and  then  converted  into  cumulative 
frequency  diagrams  by  a  computer  program  written  for  that 
purpose.  The  repair  times  for  each  of  the  bits  were  not 
determined  from  actual  data  but  were  instead  estimated  by 
company  officials.  The  repair  times  for  each  of  the  bits 
on  the  machine  were  individually  assumed  to  be  constant 
values  but  higher  constant  repair  times  were  assigned  to 
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bits  near  the  periphery  of  the  cutting  head  where  working 
conditions  were  more  difficult  due  to  space  restrictions. 

The  time-to -failure  and  repair  distributions  for  the  tunnel¬ 
ing  device  were  determined  by  reconstructing  the  operating 
record  from  shift  reports  and  obtaining  the  individual  times 
between  failure  and  the  number  of  manhours  required  to  com¬ 
plete  each  repair.  These  were  then  formed  by  computer  into 
the  necessary  distributions  for  use  in  the  computer  program. 

The  data  for  testing  of  the  cyclic  materials  handling 
subsystem  was  not  hard  to  gather,  although  actual  field  data 
was  not  available  for  some  of  the  variables.  A  tunnel  pro¬ 
file  with  many  grade  changes  was  hypothesized  for  use  in  the 
test.  Trains  corresponding  to  present  practice  were  assembled 
for  the  simulation.  Three  two-speed  diesel  locomotives  with 
a  weight  of  fifteen  tons  were  selected.  Eight  fifteen-ton 
cars  with  an  empty  weight  of  three  tons  were  chosen  for  each 
train.  The  distribution  of  the  weight  loaded  in  each  of  the 
trains  was  assumed  to  be  normal  with  a  standard  deviation 
equal  to  5%  of  the  mean  value.  A  bimodal  switchout  time  dis¬ 
tribution  was  hypothesized  to  ind icate  a  practice  of  unload¬ 
ing  supplies  from  the  incoming  trains.  The  distribution  of 
dumping  time  was  estimated  from  one  contractor’s  experience 
on  a  previous  tunneling  project. 

Data  for  the  tunnel  support  and  the  environmental  con¬ 
trol  subsystems  was  obtained  from  available  records  on  the 
White  Pine  system.  The  individual  samples  were  collected 
by  studying  the  shift  reports  and  extrapolating  as  best  as 
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possible  the  number  of  manhours  spent  during  specific  acti¬ 
vities  involving  each  of  the  subsystems.  By  collecting  in~ 
formation  on  a  large  number  of  occurrences,  distributions  of 
the  manhours  required  for  specific  advances  of  these  two 
subsystems  were  formed. 

Testing  Procedure 

The  initial  test  of  the  program  was  made  with  the  idea 
of  eliminating  the  programming  problems  in  the  model,  i.e., 
eliminating  the  bugs  and  errors  in  logic  in  the  model.  This 
was  accomplished  simply  by  attempting  to  run  the  program  and 
check  the  validity  of  the  results.  The  most  complex  portion 
of  the  program  was  the  materials  handling  subsystem  and  this 
subsystem  was  the  most  difficult  to  debug.  When  the  obvious 
debugging  problems  were  out  of  the  way9  the  program  was  then 
checked  to  be  certain  it  was  operating  logically  and  out- 
putting  data  in  the  log  of  operations  which  agreed  with  cal¬ 
culations  made  by  hand.  This  procedure  probably  did  not  re¬ 
sult  in  testing  all  the  possible  branches  of  the  program 
even  though  an  attempt  was  made  to  cover  as  much  of  the  logic 
as  possible.  After  several  problems  were  eliminated,  the 
program  seemed  to  be  at  least  superficially  correct  and  free 
of  obvious  bugs. 

No  attempt  was  made  to  test  the  accuracy  of  the  simula¬ 
tion  model  in  terms  of  the  overall  results  as  this  step  in 
the  testing  procedure  was  planned  for  the  second  year  of  the 
project.  The  testing  of  the  accuracy  of  the  model  was  to  be 
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undertaken  using  the  data  obtained  at  White  Pine  as  input 
to  a  simulation  run  which  would  model  the  tunneling  opera¬ 
tion  for  about  one  month's  time.  The  results  of  the  simu¬ 
lation  in  terms  of  the  tunnel  advance  and  the  times  spent 
in  the  various  unit  operations  would  then  be  compared  with 
the  actual  values  of  these  variables  obtained  from  tunneling 
records  for  the  time  period  in  question.  Attempts  could  then 
be  made  to  adjust  or  improve  the  computer  program  in  areas 
where  its  performance  was  concluded  to  be  unsuitable. 

Present  Status  of  the  Program 

Since  the  development  of  the  program  is  not  complete 
at  the  present  time,  users  should  recognize  that  parts  of 
the  model  may  still  be  in  rather  unfinished  form  in  the  pro¬ 
gram.  In  particular,  the  program  may  still  contain  bugs  which 
have  not  been  detected.  In  addition,  options  which  would 
make  the  program  more  versatile  and  useful  may  not  be  included 
due  to  the  limited  period  of  use  of  the  model.  As  an  example, 
it  was  hoped  to  expand  the  program  to  include  the  logic  for 
systems  using  both  cyclic  and  continuous  materials  handling 
systems,  the  cyclic  system  being  applied  to  the  handling  of 
supplies  while  the  continuous  system  was  applied  to  the 
handling  of  muck.  Such  logic  does  not  presently  exist  in 
the  model.  These  inadequacies  are  to  be  taken  care  of  during 
the  latter  stages  of  development  and  use  of  the  program. 

,  the  program  is  still  in  a  state  of 


At  present,  however 
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development  and  testing  and  should  not  be  considered  a  fin¬ 
ished  product. 

One  of  the  most  important  aspects  of  the  testing  of  the 
model  which  has  not  been  completed  is  the  testing  of  the 
accuracy  of  the  program  in  modeling  actual  tunneling  situa¬ 
tions.  For  this  reason.,  the  fact  that  the  model  will  com¬ 
plete  a  run  and  output  data  is  not  sr.ficient  reason  to  have 
complete  confidence  in  the  results.  Inaccuracies  may  be 
caused  by  bugs  in  the  program  or  by  the  assumptions  of  the 
model  not  being  valid  for  all  or  some  of  the  conditions  under 
vhich  the  model  is  to  be  applied.  Users  should  note  these 
warnings  before  making  use  of  the  program. 
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Definition  of  Important  Variables  in  the  Program 

This  section  contains  the  definition  of  the  most  important 
variables  and  the  units  in  which  they  are  expressed  in  their 
usage  in  the  program.  The  variables  which  must  be  input 
into  the  program  are  also  defined  in  this  list.  To  prepare 
an  input  data  deck,  a  user  must  refer  to  the  main  program 
and  SUBROUTINE  TRANS  where  all  the  data  is  input.  The  data 
prepared  for  the  main  program  should  appear  first  in  the 
data  deck  while  the  data  for  SUBROUTINE  TRANS  follows. 

All  the  input  variables  are  defined  in  this  list,  which  is 
alphabetized  for  convenience  in  locating  specific  variable 
names.  Users  may  refer  to  the  program  for  the  order  and 
format  information  on  variables  and  then  to  this  list  for 
the  definition. 

ACCFC  —  available  accelerating  force,  tons 

ACCMAX  —  maximum  acceleration  rate  allowed  in  the  tunnel, 
feet  per  second  per  second 

ACCR  —  acceleration  rate  of  a  train,  feet  per  second  per 
second 

ACT(I) _ the  reduced  time  to  complete  activity  i_  after 

redistributing  the  manpower. 

ACTIM(I)  —  the  time  required  to  complete  the  i_th  activity 

ADRT  --  the  tunnel  advance  rate,  feet  per  hour 

AFT  _  feet  of  advance  required  to  load  one  train 

AVAMH  —  manhours  available  for  the  support  function 

AVATF  —  available  tractive  effort  of  a  locomotive  at  its 
current  speed,  pounds 
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BINSMH  —  manhours  required  to  inspect  the  bits  and  regrip  tho 
mole 

CAPMH  --  the  capacity  of  the  continuous  materials  handling 
system,  tons  per  minute 

CCTM  —  the  time  in  minutes  required  to  generate  one  trainload 
of  muck 

CF(I,J)  —  the  ^th  ordinate  value  read  in  from  the  ith  cumulative 
probability  curve  in  the  main  program 

I  =  1  to  NBITS  correspond  to  the  probability 

distributions  for  the  time-to-f ailure  of  the 
bits  in  feet  or  operating  hours. 

1  =  NBITS  +  1  to  2 aNBITS  correspond  to  the  probability 
distributions  for  the  manhours  required  to  replace 
the  bits 

I  =  2*NB ITS  +  1  corresponds  to  the  probability  distrib 
ution  for  the  time  between  bit  repairs,  hours 
I  =  2*NBITS  +  2  corresponds  to  the  probability  distrib 
ution  for  the  time  between  mole  repairs,  hours 
I  =  2*NBITS  +  3  corresponds  to  the  probability  dis¬ 
tribution  for  the  time  between  repairs  of  the 
third  (optional)  equipment,  hours 
I  =  2*NBITS  +  4  corresponds  to  the  probability  dis¬ 
tribution  for  the  advance  rate,  feet  per  hour 
I  =  2*NBITS  +  5  corresponds  to  the  probability 

distribution  for  the  manhours  for  repair  of  the 
bits 

I  =  2*NBITS  +  6  corresponds  to  the  probability  dis¬ 
tribution  for  the  manhours  required  for  repair 
of  the  mole 

I  =  2*NBITS  +  7  corresponds  to  the  probability  dis¬ 
tribution  for  the  manhours  required  to  repair  of 
the  third  (optional)  equipment 
Altogether,  2*NBITS  +  7  probability  distributions 
are  read  into  the  main  program. 


CFD(I) 

CFL(I) 

CFR(I) 

CFS(I) 

CROSEC 


the  ith  ordinate  value  read  in  from  the  cumulative 
probability  curve  for  the  dumping  time 

the  ith  ordinate  value  read  in  from  the  cumulative 
probability  curve  for  weight  of  the  muck  in  a  car 

the  ith  value  of  the  cumulative  probability  read 
from  the  support  requirement  function 

the  ith  ordinate  value  read  in  from  the  cumulative 
probability  curve  for  switching  time 

cross-sectional  area  of  the  tunnel,  square  feet 


CT(I)  --  the  ith  abscissa  value  read  in  from  the  cumulative 
probability  curve  for  dumping  time,  minutes 
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CTIME  --  clock  time  from  the  start  of  the  simulation,  minutes 

CTLOC(I)  --  total  or  clock  time  of  the  ith  locomotive  in 
minutes 

CTM(I)  --  the  i^th  value  of  the  ACT(I)  array  if  arranged  in 
ascending  order 

D1  --  the  distance  in  feet  from  one  stop  to  the  next  stop  of 
a  train  excluding  the  distance  required  to  stop 

D2  —  distance  in  fee*;  to  the  end  of  the  present  grade  section 

D ( I )  --  horizontal  length  of  section  i_  of  the  tunnel  profile  in 
feet 

DECEL  --  maximum  deceleration  rate  allowed  in  the  tunnel,  feet 
per  second  per  second 

DELTH  —  increment  added  to  the  tunnel  length  as  the  face 
advances,  feet 

DISTR(I)  —  distance  traveled  by  the  i^th  locomotive  in  feet 

DISW  —  distance  between  two  switching  points  in  feet 

DMS  —  current  distance  between  the  switch  closest  the  face 
and  the  next  switching  point,  feet 

DS(I)  —  distance  from  the  dumping  station  to  the  i^th  switch 

DSTOP(I)  —  distance  required  for  the  i^th  locomotive  to  stop, 
feet 

FCAR  —  the  friction  coefficient  of  each  mine  car  in  pounds 
per  ton 

FLOCO(I)  —  friction  coefficient  of  locomotive  i^  in  pounds 
per  ton 

FRFC  —  force  required  to  overcome  the  frictional  resistance, 
pounds 

FTA(I)  —  the  ith  abscissa  value  read  from  the  support 

requirement  curve,  manhours  per  foot  of  advance 

G(I)  —  present  grade  of  section  i^  of  the  tunnel  profile 

GAMMA  --  specific  weight  of  the  muck  in  the  solid,  pounds 
per  cubic  foot 

GFC  —  force  required  to  overcome  the  grade  resistance, 
pounds 

GLEFT  • —  distance  in  feet  remaining  to  be  traveled  in  the 
track  section 
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HAUL  —  the  current  haulage  length  in  feet 

HRPSH  —  working  hours  per  shift,  i.e.,  the  total  shift  time 
minus  travel  and  other  idle  time 

ICYCLE  —  the  variable  which  indicates  the  type  of  material 
handling  system 

ICYCLE  =  0  indicates  a  continuous  system 

IDE ( I )  —  queueing  number  of  the  i^th  locomotive  while  waiting 
empty  at.  the  dumping  station  to  enter  the  tunnel; 
IDE(I)  =  0  means  the  i^th  locomotive  is  not  in  the 
queue 

IDEOS  --  the  variable  which  indicates  that  the  simulation  is 
to  terminate;  IDEOS  =  1  indicates  the  termination 

IDL(I)  —  queuing  number  of  the  i^th  locomotive  as  it  waits  to 

dump  its  muck  at  the  dumping  station;  IDL(I)  =  0  means 
the  ith  locomotive  is  not  in  the  queue 

IDLOAD  —  indicates  whether  any  trains  were  loaded  or  not; 

IDLOAD  =  1  indicates  trains  have  been  loaded 

IL  —  the  number  of  the  locomotive  which  has  the  shortest  clock 
time  but  which  is  awaiting  the  movement  of  another  loco¬ 
motive 

ILC  --  the  number  of  the  locomotive  which  has  the  same  clock 
time  as  that  of  the  main  program 

ILS  —  controls  the  input  statements  in  SUBROUTINE  TRANS; 

ILS  =  0  means  no  simulation  is  performed 

ILWTID  —  the  variable  which  indicates  the  beginning  of  the 
simulation;  ILWTID  =  1  indicates  the  beginning 

IMAN  --  number  of  men  currently  available 

INLC  —  the  number  of  the  loaded  locomotive  at  the  loading 
point 

INSPM  --  the  number  of  men  required  to  i  sped  the  bits 

IR  —  the  subscript  used  to  obtain  the  repair  manhours  for 
ITEM 

ITEM  —  the  number  of  the  unit  which  has  the  shortest  life 

KK  —  the  next  lower  speed  point  on  the  characteristic  curve 

KMAX  —  number  of  points  on  the  characteristic  curves  of  the 
locomotive  at  which  input  data  will  be  read 
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KOUNT  —  the  number  of  bits  which  need  to  he  replaced 

LC(I)  —  the  queuing  number  of  the  trains  in  the  ILC  list 

LCLAS  --  number  of  points  read  in  from  the  cumulation  fre¬ 
quency  function  for  the  weight  of  muck  in  one  muck 
car 

LIL  --  the  variable  which  retains  the  numbers  of  the  loco¬ 
motives  which  were  in  the  previous  IL  list 

LL(I)  —  the  switch  on  which  the  jLth  locomotive  is  located 

LLW(I)  --  the  number  of  the  locomotive  in  the  jLth  spot  in 
the  LIL  queue 

LOAD(I)  —  indicates  the  status  of  the  i^th  train 

I.0AD(I)  =  0  indicates  the  train  is  empty 
LOAD(I)  =  1  indicates  the  train  is  loaded 


LOGPRT  --  print  option  variable 

LOGPRT  =  0  indicates  that  the  complete  log  of 
operations  is  printed 

LOGPRT  f  0  indicates  that  only  the  summary  of 
the  simulation  is  printed 


LS  ( I) 


the  variable  which  indicates  the  status 
switch 

LS(I)  -  0  indicates  the  switch  is  empty 
LS(I)  =  1  indicates  the  switch  contains 
train 

LS(I)  =  2  indicates  the  switch  contains 
train 

LS(I)  =  3  indicates  the  switch  contains 
empty  and  a  loaded  train 


of  the  ith 

an  empty 
a  loaded 
both  an 


LW(I) 


the  number  of  the  i^th  locomotive  in  the  clock  time 
queue 


LWTID  --  indicates  whether  or  not  there  is  an  empty  train  at 

the  loading  point;  LWTID  =  0  indicates  no  empty  train 


MAD  —  number  of  men  available  to  be  reassigned  when  a  repair 
activity  is  completed 


MAN(I,J)  --  the  variable  which  stores  the  upper  and  lower  limits 
on  the  number  of  men  assigned  to  each  activity 
1=1  corresponds  to  the  lower  limit 
1=2  corresponds  to  the  upper  limit 
J  =  1  to  NBITS  corresponds  to  the  limits  of 
manpower  for  the  replacement  of  the  bits 
J  =  NBITS  +  1  corresponds  to  the  limits  of 
manpower  for  the  repair  of  the  bits 
J  =  NBITS  +  2  corresponds  to  the  limits  of  manpower 
for  the  repair  of  the  mole 
J  =  NBITS  +  3  corresponds  to  the  limits  of  manpower 
for  the  repair  of  the  third  (optional)  equipment 
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MANAW(I)  —  the  number  of  men  assigned  to  the  i^th  job 

MAXSHT  --  maximum  number  of  shifts  that  the  simulation  is  to 
be  run 

MH  --  variable  which  indicates  which  option  was  employed  in 

reading  in  the  muck  generation  cumulative  frequency  curves 
MH  =  0  indicates  the  abscissa  values  are  in  terms  of  hours 
MH  £  0  indicates  the  values  are  in  terms  of  the  feet  of 
advance 

MM  —  grade  section  number  which  train  NL  is  presently  traversing 

ML  --  number  of  the  locomotive  currently  being  moved 

MNBITL  —  lower  limit  on  the  number  of  men  required  to  repair 
bits 

MNBITU  —  upper  limit  on  the  number  of  men  required  to  repair 
bits 

MOTM  —  the  time  in  minutes  required  for  the  hauling  of  the  muck 
generated  by  TEMSTR 

MREST  —  cumulative  number  of  men  who  spent  idle  time  during 
the  computer  run 

MSS(I)  —  number  of  the  locomotive  occupying  the  i^th  switch 

MTB  —  number  of  men  who  are  reassigned  when  a  repair  activity 
is  completed 

NACF  —  the  number  of  events  to  be  simulated  in  the  muck  genera¬ 
tion  subsystem  in  addition  to  the  events  related  to  bit 
replacement 

NBITS  --  the  number  of  bits 

NCARS  —  number  of  muck  cars  assigned  to  each  train 

NCF  —  total  number  of  cumulative  frequency  diagrams  read  into 
the  muck  generation  subsystem 

NCLAS(I)  --  the  number  of  points  read  in  for  the  jLth  cumulative 

probability  function  of  the  muck  generation  subsystem 

THIRD1  —  cumulative  time  spent  in  doing  the  third  event,  minutes 

NCREW  —  the  number  of  men  in  the  crew 

NDC  --  number  of  points  read  in  from  the  cumulative  frequency 
function  for  the  dumping  time 

NEVENT  —  the  number  of  separate  repair  activities  currently 
being  performed 
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NL  —  locomotive  number  presently  being  simulated 

NLDL  --  number  of  loaded  trains  waiting  at  the  dumping  station 
to  dump 

NLDE  —  number  of  empty  trains  at  the  dumping  station 

NLDL  _  the  number  of  loaded  trains  at  the  dumping  point 

NLOCO  —  number  of  locomotives 

NRBG  —  the  number  of  points  read  in  from  the  cumulative 
probability  curve  for  the  support  function 

NS  --  the  switch  from  which  locomotive  NL  is  moved 

NCS  —  number  of  points  read  in  from  the  cumulative  frequency 
function  for  the  time  to  switch  trains  behind  the  mole 

NSCF  —  the  number  of  time-between-r epair  cumulative  probability 
functions  read  into  the  muck  generation  subsystem 

NSDP  —  number  of  sections  of  the  haulage  profile  between  the 
dumping  point  and  the  tunnel  mouth  read  into  the 
program 

NSECS  —  number  of  sections  of  the  haulage  profile  within  the 
tunnel  read  into  the  program  (after  input,  NSECS  is 
the  number  of  sections  in  the  tunnel  profile  at  the 
time  of  simulation) 

NSW  —  number  of  switching  points  currently  in  the  haulage 
system 

NSHIFT  _  the  number  of  shifts  simulated  so  far  in  the  current 

run 

OTRD  —  distance  in  feet  that  the  train  overtravels 

pyx  _  the  time  the  continuous  materials  handling  system  can 

operate  before  a  breakdown,  minutes 

RADIUS  —  radius  of  the  tunnel,  feet 

REQMH  —  required  manhours  of  support  work  for  one  foot  of 
advance 

REQTF  —  required  tractive  effort,  pounds 
RESTMH  —  cumulative  number  of  idle  manhours 

S  (1 ,  J)  —  speed  of  the  i  th  locomotive  at  the  j_th  point  on  its 
characteristic  curve 

_  cumulative  length  of  advance  since  the  last  value  of 

DELTH  was  added  to  HAUL 


SAFT 
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SCCTM  —  the  cumulative  time  in  minutes  to  advance  by  TEMSTR 

SGL(I)  —  distance  in  feet  from  the  i_th  switch  to  the  inby 

end  of  the  track  section  on  which  the  switch  exists 

SLEFT  —  distance  in  feet  to  the  next  switch  point 

SP  --  former  3peed  of  the  train,  feet  per  second 

SPEED  (I)  --  velocity  of  the  i_th  locomotive,  feet  per  minute 

SSCC  --  incremental  time  in  minutes  that  a  train  waits  for 
the  completion  of  another  event 

ST(I)  --  the  i^th  abscissa  value  read  in  from  the  cumulative 
probability  curve  for  switching  time,  minutes 

STROKE  --  stroke  of  the  mole,  feet 

SWTTIM  --  the  cumulative  delay  time  in  minutes  due  to  the 
support  subsystem 

T('I,J)  --  tractive  effort  of  the  i^th  locomotive  at  the  j_th 
point  on  its  characteristic  curve 

T1  --  the  time  in  seconds  required  to  travel  the  distance  D1 

T 2  —  time  in  seconds  to  reach  the  end  of  the  present  grade 
section 

TBELT1  --  operating  time  of  the  continuous  materials  handling 
system,  minutes 

TBELT2  --  delay  time  due  to  the  continuous  materials  handling 
system,  minutes 

TBELT3  —  downtime  of  the  continuous  materials  handling 
system,  minutes 

TBIT1  --  cumulative  working  time  of  the  bits,  minutes 

TBIT2  --  cumulative  idle  time  of  the  bits,  minutes 

TBIT 3  --  cumulative  time  the  bits  are  under  repair,  minutes 

TBIT4  --  cumulative  time  the  bits  are  under  replacement, 
minutes 

TBIT5  --  cumulative  time  the  bits  are  under  inspection, 
minutes 

TDUMP(I)  --  dumping  time  in  minutes  of  the  i^th  locomotive 
during  the  last  dumping  cycle 

TEMPWT  —  the  weight  in  tons  of  the  portion  of  the  material 
remaining  to  be  loaded  in  the  current  train 
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TEMSTR  --  portion  of  the  stroke  which  remains  to  be  completed 

TFT  --  the  number  of  feet  the  mole  can  advance  before  being 
stopped 

TFTA  --  the  incremental  number  of  feet  the  mole  is  to  be 
adv  anc  ed 

THIRD!  --  time  the  third  (extra)  subsystem  spends  working, 
minutes 

THIRD2  —  time  the  third  (extra)  subsystem  spends  in  waiting, 
minutes 

THIRD3  —  time  the  third  (extra)  sybsystem  undergoes  repair, 
minutes 

TIMAX  --  maximum  clock  time  in  minutes  that  the  simulation  is 
to  be  run 

TIME(I)  --  time  required  in  minutes  for  the  i_th  locomotive  to 
get  from  one  switch  to  the  next  minus  the  value  of 
TP AS S ( I)  or  TSTOP(l) 

TLOAD(I)  —  loading  time  in  minutes  of  the  i^th  locomotive  when 
it  was  last  loaded,  minutes 

TLOCl(I)  —  cumulative  time  the  i_th  locomotive  spends  in  the 
loading  process,  minutes 

TL0C2(I)  --  cumulative  time  the  f^th  locomotive  spends  in  the 
dumping  process,  minutes 

TL0C3(I)  —  cumulative  time  the  jLth  locomotive  spends  in 
motion,  minutes 

TL0C4(I)  --  cumulative  time  the  _ith  locomotive  spends  waiting, 
minutes 

TMH  --  manhours  requirsd  to  advance  by  TFTA 

TM0LE1  --  cumulative  working  time  of  the  mole,  minutes 

TM0LE2  --  cumulative  idle  time  of  the  mole,  minutes 

TM0LE3  --  cumulative  time  the  mole  is  under  repair,  minutes 

TNL  —  maximum  length  of  advance  of  the  tunnel  in  feet  for 
the  simulation  run 

--  the  tolerance  placed  upon  the  repair  starting  times 
in  minutes;  i.e.,  when  one  repair  action  is  initiated, 
the  potential  repairs  are  checked  and  are  also  initiated 
if  they  are  within  the  tolerance  time  of  requiring 
repair 


T0L1T 
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TPASS(I)  --  time  required  for  the  ith  locomotive  to  travel 
through  a  switch  without  stopping,  minutes 

TPM  —  the  muck  generation  rate  in  tons  per  minute 

TSEC  —  the  current  length  of  the  ith  section  which  has  been 
driven  and  added  to  the  variable  HAUL 

TSTOP(l)  —  time  required  for  the  ith  locomotive  to  decelerate 
and  stop  on  a  switch,  minutes 


TSUPPT 


cumulative  time  expended  for  support  activities 
minutes  5 


the  time  in  minutes  required  to  switch  out  the  loaded 
train  at  the  loading  point 

the  time  in  minutes  that  the  mole  can  advance  before 
being  stopped 


TUNNEL 


the  length  of  tunnel  bored  to  the  present,  in  feet 
from  the  portal 


TV  ( I ,  J)  --  the  j_th  time  or  other  abscissa  value  read  in  from 
the  j^th  cumulative  probability  curve  in  the  main 
program  in  units  of  feet  or  operating  hours  (for 
a  definition  of  the  meanings  of  each  of  the  values 
of  1,  see  the  variable  CF(1,J)) 

VELMAX  —  maximum  velocity  allowed  in  the  tunnel,  feet  per 
second 


WAITIM 

WTCAR 
WTD  — 
WTG  — 
WTIM  — 


cumulative  idle  time  of  the  muck  generation  sub¬ 
system  in  minutes 

-  weight  in  tons  of  each  muck  car  while  empty 

cumulative  weight  of  muck  dumped,  tons 

cumulative  weight  of  muck  generated,  tons 

the  time  in  minutes  to  move  an  empty  train  to  the 
loading  point 


WTL(I)  — 

WTLDG  — 
WTLOAD  ( I ) 
WTLOC(I) 
WTMUCK  — 


the  i^th  abscissa  value  read  in  from  the  cumulative 
probability  curve  for  weight  of  muck  in  a  car,  tons 

the  weight  in  tons  of  the  load  to  be  generated  by  TFTA 

—  weight  of  the  muck  in  the  _ith  train  in  tons 

--  weight  of  locomotive  i.  in  tons 

the  weight  of  muck  in  tons  to  be  loaded  in  one  train 

weight  in  tons  of  the  _i_th  locomotive  and  its 
cars 


WTTRN(I) 


empty 


WWTM 


-incremental  time  in  minutes  that  a  train  waits  for 
the  completion  of  another  event 


XX(I) 


the  abscissa  value  as  determined  from  SUBROUTINE 
CALCUM 
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Computer  Logic  Diagrams 

The  Computer  logic  diagrams  of  the  main  program  and 
two  of  the  seven  subrotines,  SUBROUTINE  MOTION  and  SUBROUTINE 
TRANS,  appear  on  the  following  pages.  The  remainder  of 
the  subroutines  are  not  represented  in  this  section  since 
they  perform  relatively  simple  functions  for  which  the 
logic  diagrams  were  considered  unnecessary.  The  diagrams 
presented  are  not  intended  to  be  a  detailed  flowchart  of 
all  the  calculations  and  manipulations  that  take  place 
in  the  computer  program.  Instead,  they  are  meant  to  convey 
the  macro  logic  of  the  simulation  and  way  that  it  fits 
together  in  the  model.  Most  of  the  variables  which  appear 
in  the  logic  diagram  are  identified  in  the  previous  section 
of  this  Appendix.  In  the  logic  diagrams,  two  types  of 
offpage  connectors  are  used.  The  connectors  appear  as 
small  circles  with  numbers  or  letters  enclosed.  Connectors 


containing  numbers  indicate  the  actual  program  statement 
at  which  the  connection  is  Co  be  made.  This  gives  the  reader 
one  extra  bit  of  help  in  following  the  program  using  the 
logic  diagram.  The  connectors  containing  letters  are  those 


for  which  no  exact  statement  number  to  which  the  program 
proceeds  could  be  named. 


- ~ - - - -  - - - — - -  - rw—XW*  - 
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MAIN  PROGRAM 


r 

I 


i. 


>Sf 


‘Kj^VffTKiVfrwr- 


MAIN  PROGRAM  (CONT'D) 


MWWWWVWS 


MAIN  PROGRAM  (CONT'D) 


CALCULATE  WTLDG j 
SCCTM  =0.0 


NO  /  ICYCLE  XYES 
V  EQ  0  / 


,siwsm 
EO-O 
OR  1 


fK 

LWTID  =  0  I 

r 

CALL  TRANS 


DETERMINE  CTIME, 
THIRD2,  TM0LE2 
TBIT2 


WDLDG=  WTLDG-TEMPWT 
TEMPWT  =  0.0 
ILWTID=  ILWTID  + 1 
I DLOAD  =  0 


(CALL  TRANS 


LWTID  =  t 


_i/WTLDG- 

\tempwt 


TEMPWT  =  0.0 
ILWTID=  ILWTID +  1 
IDLOAD= 1 


-  FTAD  =  TFT 
|  CALL  TRANS 


TEMPWT  =  TEMPWT- WTLDG 
FIND  LOCO  JJ  AT  NSW 
FIND  TLOCKJJ),  CTIME (JJ) 


ideos\X§s/'~\ 
•  EQ.  1 


CALL  TRANS 


CALL  WIUCK 


TEMPWT  = 

WTMUCK 

_ J 

X  loeos^orES 

\-E0  1  y? 

JNO 

CALL  MUCK 

— 

TEMPWT  =  WTMUCK  I 


VW]  trWNVR  « 


l! 


i  i 


T»» 


MAIN  PROGRAM  (CONT'D) 


MAIN  PROGI 


DISTRIBUTE  MANPOWER 


YES 


TEMSTP^TEMSTR-TPT 


FIND  CTIME,  TBIT3,  TM0LE3 

THIRD3,  TBIT2,  TMOLE2,  THIRD2 
XX(I  )=XX(I)-XX(ITEM) 
XX(ITEM) 


FIND  SHFTMH, 
WORK(I),  ACTIM(I) 
ACT(I), 


H 
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RAM  (CONT'D) 


PRINT  REPAIR  DATA 

j  INSPECT  BITS 
,TEMSTR= STROKE 


6 


*rr*jua«r.”T«*.r*TV.‘:w 


SUBROUTINE  MOTION 


SUBROUTINE  MOTION  (CONT’D) 


S) 
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SUBROUTINE  MOTION  (CONT'D) 


ACCR-\+ 
ACC MAX  - 


FIND  D1,T1 
FIND  02, T2 


ACCR=ACCMAX 


oi-  X+ 

D2  s>~Q 


SPEED(NL)  =  S(NL,KK+1 ) 


D1  =  D2 
T1  =  T2 

SP  =  SPEED(NL) 
SPEED(NL)  =  SP  +  T1*ACCR 


CALCULATE  DSTOP(NL) 
DISTR(NL),  TIME(NL) 
GLEFT,  SLEFT 


SLEPT 

-Q.S 


GLEFT  V 
-0.3  /C 


CALCULATE 
GFC,  REQTF 


SUBROUTINE  MOTION  (CONT'D) 


70 


SUBROUTINE  TRANS 


SUBROUTINE  TRANS  (CONT'D) 
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CALCULATE  CAPMH 
I  TPM,  MOTM 


I  CALCULATE  OIFTM, 
MOTM,  CTIME, 
SHFTMH,  NSHIFT 


MOTM  =  MOTM-SPWT 
FIND  T8ELT1 
SPWT=0.0 

CALCULATE  DOWNMH 
FIND  TBELT3 
FIND  NEW  SPWT 


PRINT  SUMMARY 
OF  CONTINUOUS 
[HANDLING  SYSTEM 


READ  INPUT 
DATA 

HCZ 

INITIALIZE 

VARIABLES 

rr~~ 

CALCULATE 
PWT,  SPWT 


SPWT  =  SPWT- MOTM 
FIND  TBELT1 


SUBROUTINE  TRANS  (CONT’D) 

fSi 


(3»1 

Jnq 

1s{ne1 

■EO.  1 
.OP  3. 


t  LC 
■EQ.O 


FIND  INLC 


JtJLC 

EG 

LCC11 


1oad(nlNyes_ 

EO  1  / 


NS-EO.1 


VES  X  NS  EO 

^S.  NSW 

YES^^  ,£q  | 


.SINS*! 
EO  Z 
on  3^ 


COUNT  CTIME 
DISTR(NL).  TL0C3(NL) 


1 

' —  ILC=ILC+1 

- — n 

I  l 

LC(ILC)SNL 

CALL  MOTION  j 

-  IfiUfJ— 

I 

PRINT  MOTION 


COUNT  CTLOC(NL) 
TL0C3(NL) 
CHANGE  LL(NL) 
LS(NS _ _ 


LOGPRT 


MS  EO 
NSW 


t  LC 
EO.O 


FIND  LOCO  II 
AT  LS(NS-I) 


lc{h  x 

VEG  11 1/ 
yesT" 
!lc*ilc+i 

LC(ILC)  =  NL 
I DLC  =  1 
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SUBROUTINE  TRANS  (CONT'D) 


SUBROUTINE  TRANS  (CONT'D) 


ftWCVtStr- 
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SUBROUTINE  TRANS  (CONT'D) 
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SUBROUTINE  TRANS  (CONT'D) 


■k. 

'umr.n'  — 


NLOCO 


|CTTM  =  CTIME 
+SCCTM 
FIND  ILC 


FIND  NL  FOR 
SMALLEST  CTLOC(I), 
IFILC _ 


».*w±vw 
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SUBROUTINE  TRANS  (CONT’D) 


SUBROUTINE  TRANS  (CONT'D) 


78 


«*wm**. OTK***'  «•  ffnw  S55T«” 
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SUBROUTINE  TRANS .  (CONT'D) 


DEtHLi\ HO 
■EQ  T 


LS(NS+I ) 
dEQOORl 


FIND  CTLQC  NL), 
TLOC  1  (NL) 
CHANGE  LL(NL), 
LS  (NS) 


PRINT  MOTION 
DATA 


FIND  CTLOC(JJ), 
TL0C4(  JJ ) 


FIND  SSCC 
TL0C4(NL), 
CTLOC(NL) 
ILC  =  ILC+NL 


CHANGE  IL 
TO  ILC 
DETERMINE 
TLOC  4(J  Jh 
CTLOC  JJ 


FIND  LOCO  II 
ON  LS  (NS+1 ) 


CHANGE 

NLDE,IDE(I) 


DETERMINE 

ILC.CTLOC(I) 


CALL  MOTION 


^•.■WWVHSMW wnMWWfl* 
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The  Computer  Program 

The  computer  model  which  is  presented  on  the  following 
pages  consists  of  a  main  program  and  seven  subroutines. 

The  jobs  performed  by  the  individual  subroutines  are  explained 
by  the  comment  cards  located  at  the  beginning  of  each  subroutine 
deck.  The  program  was  written  for  the  Univac  1108  in  standard 
FORTRAN  IV  and  should  be  rather  easy  to  transfer  to  other 
machines  since  few  machine  dependent  statements  were  used. 

One  aspect  of  the  program  which  may  need  attention  is  the 
random  number  generator.  The  Univac  1108  used  at  the  Univer¬ 
sity  of  Utah  uses  the  function  RAND(N)  to  assign  a  random 
number  uniformly  distributed  between  zero  and  one  to  any 
variable.  For  example,  the  statement  Y=RAND (N)  will  result 
in  a  random  number  between  zero  and  one  being  assigned 
to  the  variable  Y.  Users  wishing  to  use  the  program  will 
have  to  check  the  random  number  function  for  their  machines 
and,  if  necessary,  replace  all  the  statements  calling  random 
numbers  with  statements  specific  to  their  own  machines. 


r>  n  r 
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DIMENSION  I XRf<  ( 20 ,»  » IXa (20 ) 
l) I mens i on  man ( i oa ,  '■> i  ju/..., 

mS:  £i!t;,u,*cy' 

»^LMAA^iIEctU™Sc!L^M^T?'lAUU,rtECSXNsS*2TCARLF^MJC'ACCMAX, 

*CR0S£C»  AMMA,cTI‘-tE:#Lw  I  JO#  TT^'inl^VimTx*  aIRHcTra^* 

*FL0C0f liTTfef iSo;^DClOO>?DS«Vo{ J WTLOAD?So > isPEelf 1 f  J°> ' ^TLOr  C 101 » 

fa  CF(l2H»  13)  ,T\i  M?n.Mi  T  in,  »4FS(13)  #TLOAlH10)  #TDUMH(10)  , 

*  *  {0l-0AU # SCCTM,  TB I T2 ,  THOLE 2  #  THJLKD2 1  dELTh !  IOFOS » TfJLT  WTID 

*  ^rAD'7MOLEl,7R1T1'TH1KDl'fJCR^W»SHFTMH,HRFSHrWSHIF-T,lMAN,ICYCLE 

head  IN  CONTROL  CARDS 


Rt  AO  94  #  MRUG 

94  FORMAT  (  15) 

i  cfRin-  ft*(,).  isi.nrbs  ) 

is  FORMAT  (  2F10.5) 

READ  95#  RAD I OS# GAMMA 

95  FORMAT  (  2F10.3) 
CU0SEC=3,14i59*RADlUS'FRADlUS 

READ  9b*  MBITS#  NACF,  NCREW,lcYCLE 
9o  FORMAT  (  4llu) 

READ  97,  T1MAX,  TNL#  hAXSHT 

97  FORMAT  (  2F15.4,  HO  ) 

RLaD  98#  LOGPRT#  MM 

98  FORMAT  (  215  ) 

Rt.AU  99#  CTIMc#  TUNNEl#  HRPSH#  TO<  tT 

99  FORMAT  (  4F10.3)  UT 

RLAD  if)0#  BINSMH#STR0XE»INSPM 


10''  FORMAT  (  2F10,3,I5) 
NSCF=NDITS*2  +  NACF 
NCF=NtiTTS*2  +  NACF  *2-l 

Initialise  the  variables 


DO  101  I=1»NCF 
DO  102  JJ=1#2 

i02  MaN( I » jj)=0 

DO  101  0=1,13 
CF ( I  # j) =0,0 

Id  TV(I#U)=o,o 
MKtSTsO 
RLSTMH=0,0 
WAI TM=o , 0 
NSHIF  r=o 


r.  r  r  t~  c.  c  nor 


•tvrztt'.ivmwzp. 
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TtiITlrO.0 
TbIT2=n , « 

TtilTizC  •  0 
TBIT4zO.O 
TdlTb=0.n 
TMOLElzu.O 
Ti'10LE2=U»G 
T.10LE3-0.0 
TiURDIsO.O 
TnIRU2=U  «  0 
THIRUjzO.C 
10=10 
IlSzO 
I  L-v  T 1  Ozu 

IulOS=0 
S  ia/  T  T  IMzO  •  0 
tsupptzo.o 
c 

c  READ  1M  CUMULATIVE  PRuBADILIT  I  FUNCTIONS 

00  103  1  =  1 »  NCK 
READ  104 r  NCLASCi) 

104  FORMAT  (  HO) 

NJ=NCi.AS(i) 

00  110  JzlrNj 

110  READ  1  Ob r  CF(l.J)*  TV(I»J) 
lCb  FORMAT  (  2F10.5) 

103  continue 

READ  IN  THE  UPPER  AND  LOWER  LIMITS  OF  NUMBER  OF  MtN 
REQUIRED  FOR  EACH  ACTIVITY 
MAM(I,l)zLOWER  LIMIT 
MAN ( I »?) sUPPER  LIMIT 

READ  10b»  MNBITL*MNBI TU 
NlsNBITS 
DO  107  IslrNl 
MAN  ( I  r  1 )  =MN|3ITL 
1C7  MAN ( I » 2 ) sMNBITU 
10b  FORMAT  (  2110  ) 

N2=N1*?+1 
N3=NSCF-1 

READ  109  ((MAw(I»J)t  vj=l»?)»  1=N2,n3) 

109  FORMAT  (  1215) 

CALL  SUBROUTINE  CALCUM  TO  OBTAIN  ABSCISSA  VALUES  CORRESPONDING 
TO  CF ( I )  =  RAND 


'  r,  r  v 


8  3 


Du  ISO  i=l»NCF 
ICF-I 

CALL  CALCUM ( ICF  » X ) 

ISC  XX(I)=X 

COMPUTE  MANHOURS  AND  uUMBER  Oh  MEN  AVAILABLE  FOR  I  HE  SHIFT 

NSMIFTrNSHlFT-f  1 
SI  lFTMH=MRPSH*NCREtf 
IMAN=NCREW 

C  II  WANTED*  PRINT  LOG  u F  OPFRAIlONS 

C 

IF  (LOGPRT)  995*lb0#995 
it"  PRINT  lbl 

lol  FORMAT  (  IH1 ,  * *****LOb  OF  OPERATIONS***** *  //) 

iuc  format  t  iho»*  clock  i ime» *5x» ’completed  event  *»3x»*  man  hour' 

*  JX,* TUNNEL  LENGTH*  ) 

c 

SEARCH  FOR  THl  SHORTEST  LIFE  OF  TilE  UNIT  IN  THE  SYSTEM 


995  IF  (  ILS  .NE.  0  )  GO  TO  99 
L 

CALL  TRANS 
C 

09  CONTINUE 

ILS=ILS+1 
ALRT=XX (NSCF) 

XX(NSCF)=XX(NSCF)/bO,U 
IF  (  ICYCLE  ,EQ.  OJ  GO  TO  829 
C 

CALL  i-il.'CK 
C 

TEMPW irrfTMUCK 
o29  TEMSTRr STROKE 
999  In=NSCF“1 

IL1=NdITS*2+2 

M=NBITS*2+1 

KK=M 

DO  170  1=;L1,IN 
Hi  IF  (  XX(I)  .LI.  XX  (M)  )  M=I 
ITEMzM 

c 

c  count  number  of  events  having  the  same  life 

c 

ICOUNTiO 
LLL=ITFM-i 
MMzITem+1 
SX=XX ( ITEM) 


P>  C'.  000 
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If  (  ITEM  .EQ.  KK  >  Go  TO  3m 
DO  311  1=KK,LuL 
XI=XX( I } 

DIF=ABF(5X-Xl) 

IF  <  DIF  .GT.  TOLIT  )  GO  TO  3li 
IC0UNT=IC0UNT+1 
I I-ICOUNT+1 
IX  ( 1 1 ) rl 

31 1  CONTINUE 

IF  (  ITEM  .EG.  IN  )  GO  TO  315 
>314  00  312  JsMM.jn 
X1=XX{JJ 
DXFsAUS(SX-Xl) 

IF  (  DIF  .GT.  TOLIT)  GO  TO  31* 

ICOUNTsICOUNT  +1 
II=IC0UNT+1 
Ia(II)=J 
ol,2  CONTINUE 
315  CONTINUE 

NEVENT.-1C0UNT+1 

IX(1)=ITEM 

COMPUTE  THE  TIME  FOR  ADVANCING  AND  THE  DISTANCE  TO  BE  ADVANCED 

IF  (MH)  355.3b6.355 

356  TTM=Xa{ITEM) 

TFT=XX{,\ISCF)*XX(ITlM) 

Go  TO  357 

355  TFT=Xx { ITEM ) 

T  T  ^=XX( ITEM)/XX{NSCF) 

357  TMHzT  TM*NCR£W  /60.G 

COMPARE  TFT  WITH  STROKE 
^67  IF  {  TFr-TEMSVR)  350,352,352 

ADVANCE  PY  TEMSTR  (  CASE  OF  TEMSTR  .LT.TFT  ) 

052  TIM=TENSTR/XX(NSCF) 

TMHsT TM*NCREW/60 . 0 
WTLDGoCKOSEC# TEMSTR  *oAMMA/20oO , 0 

tfta=temstr 

CALL  SUPPRT (TFTA  fREQMH) 

SCCTMon.O 

IF  ( ICYCLE  ,NE.  0)  GO  TO  760 
FTAD=TEMSTR 

CALL  TRANS 
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GO  TO  413 

tOO  IF  (  LS (NSW )  ,e:Q.  U  »OR.  LS(.njSW)  .HQ,!)  GO  to  751 
LwTIO=l 
GO  TO  752 

751  LWTIO=0 

CALL  TRANS 

CTIME=CTIME+WTIM 
TM0LE2=TM0LE2fWTIM 
TuIT2=Tl3lT2+WT  IM 
THIRU2=THIRD2+WTIM 

752  IF  (  WTLOG-TEMPWT)  755#754>755 

753  TEMPW1 s TEMPWT-WTLDG 
DO  756  1  =  1 f NLOCO 

756  IF  (  LL ( I )  *EU.  NSW)  oO  TO  75/ 

PRINT  758 

75«  FORMAT  (  1H0, ’LOADING  WAS  ATTEMPTED  WITHOUT  EMPTY  TRAIN  AT  LOADING 
i  POINT’ ) 

STOP 
/ 57  JJ=I 

C  TLOC ( J J ) =CTLoC ( J J ) +TMUCK 
TLCCK  JJ)=TLOCl  ( JJ)+TMUCK 
GO  TO  413 

754  TEMPWTsu.O 
ILWTIUsILWTID+1 
IDL0AD=U 


c 

C 

c 

6 


C 

C 

C 


CALL  TPmNS 

IF  <  lPLOS  ,EU.  1)  GO  TO  520 

CALL  MUCK 

TLMPWTrwTMUCK 
GO  TO  413 

755  WTLOG=wTLOG-TtMPWT 
TEMPWTrO.C 
ILWTID; ILWTId+1 
ULOAUrU 

CALL  TRANS 

IF  (  IDEOS  ,E0,  1)  GO  TO  520 
CALL  MUCK 


r 


TEMPWT=WTMUCK 
GO  TO  760 


o  n 


wMsassasjraftvpwji  =r  aww. 
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0  CHECK  TO  SEE  IF  THE  AuVANCE  CAN  BE  COMPLETED  IN  THE  SHIFT 

Hi.?  If-  (  SHFTMH-TMH)  411, <412,412 
*♦1 1  SHFTMH=SHFTMh+NCREW*HKPSH 
NSHIFTzNSHlFT+l 
GO  TO  413 

412  SfiFTMi(=bHFTMH-TMH 

TUNNELrTUNNEL  +TEMSTR 
TBIT1=TUIT1+T1M 
TMOLE1=TMOLE1+TTM 
THIRUItTHIRDI+TTM 

ctime=ctime+ttm 
AVAMH=TTM*2.0/60eC- 
IF  (  KFOMH-AVAMH)  835*835,836 
836  WTTIM=1 (REQMh-AVAMH)/u.0)*60. u 
SwTTlMrbWTTlM+WTTlM 
CTIME=C  r TME+WTTIM 
TBIT2=TBIT2+WITIM 
TMOLESs TM0LE2+WTTIM 
THIR02=THIRD2+WTTIM 
MRESTsVREST  +  (NCREW-,;) 

RESTMHrRESTMH+WTTlM*(NCREW-2)/oO.O 

waitm=waitm+wttim 

WTTMH=WTTIM*nCREW/60.0 

890  IF  (  SHFTMH-WTTMH)  830,832*832 
o30  SHFTMHrbHFTMH+NCREw^HKPSH 

NSfHFT=NSHlFT  +  l 

GO  TO  890 

«52  shftmh=shftmh-wttmh 

TbUPPT=TSUPPT+TTM+WTTIM 

GO  TO  891 

o35  TSUPP1=TSUPPT+TTM 

891  CONTINUE 

tft=tft-temstk 

TEMSTR=STROK£ 

IF  (LOGPRT)  481,480,481 
460  PRINT  ]o2 

PRINT  362,  CTIME,  TMh,  TUNNFL 
INSPECTION 


481  IF  ( IMAN-INSPM >450*451*451 


45°  mkest=iman+mrest 
RESTMHrbHFTMH+RESTMH 
IF  (  IMAM  .EQ,  0)  GO  TO  892 

WAITM=WAITM+(SHFTMH/IMAN)*6 
CTIME=CTIME+(SHFTMH/IMAN)*6 
892  NSHIFT=NSHIFT+1 

shftmh~hrpsh*ncrew 

iman=ncrew 


•  u 
.0 


•+bl  IF  (SHFTMH-BINSMH)  452*453,453 
452  SHFTMH=bHFTMH+NCRE<V*HKPSH 


ri  o  r;  r#  c,  o 
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NSH I F T =»MSH  I  FT + 1 
453  ShFTMHrSHFTMH-BlNSMH 

TblNSPs (BINSMH/INSPM  )*60.0 
TBIT5=TdlT5+TdINSP 
CT IME=CTIM£+TBlN5P 
TM0LE2=TM0LE2+TRINSP 
TMIR02=THIPD2+TPINSP 
MR£ST=(NCREW-INSPM)  +i*IREST 

RESTMH=KFSTMh+(NCREW-1NSPM)*TuINSP/60.0 

IF  (LOGPRT)  463,482*483 

482  PRINT  162 

PRINT  464»CTIME»  BINSMH*  TUNNEL 

484  FORMAT  (1H  »  F10 • 3  *  5X»*  BIT  INSPECTION  **  3X»F10.3,3a,  Flo. 3) 
COUNT  NUMBER  OF  BITS  TO  BE  REPLACED*  IF  ANY 

483  KOUNT=0 

IF  (  ID  .EQ,  1)  GO  TO  668 
TOLSTRrSTROKE+TOLIT 
GO  TO  6o9 

o66  TOLSTR=STROKE-TEMSTR+TFT+TOLII 
o69  DO  415  I=1»NBITS 

IF  (MH)  416,417,41b 
417  XX(I)=XX(I)*XX(NSCF) 

416  IF  (  XX(I)  .GT.  TOLSTK  )  GO  TO  415 
KOUNT=KOUNT+l 
K=KOUl\IT 
IbIT(K)=I 
415  CONTINUE 

REPLACE  BITS,  IF  ANY 

If  {  KOUNT  ,EQ.  0  )  GO  TO  430 
DO  420  I=l»KOUNT 
ITEM=IBIT(I) 

IF  (  IMAN-MAN(ITEM»1) )  421,422*423 

423  IF  (  IMAN-MAN( ITEM»2) )  422*422*424 
429  IMAN=1MAN+MANWK 

421  MREST=IMAN+MREST 

RESTMHrRESTMH+SHFTMH 

IF  (  IMAN  .EQ.  0)  GO  TO  893 

WAITM=WAITM+(SHFTMH/IMAN)*60.u 
CTIME=CTIME+ ( SHFTMH/IMAN) *60 . 0 
693  NSHIFTzNSHlFT+1 

shftmh=hrpsh*ncrew 

IMAN=NCREW 

424  MANWK=MAN(ITEM»2) 


:  o  o  r  o  o 
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iman=iman-manwk 
GO  TO  425 
422  MANWK=IMAN 
IMAN=0 

425  ITEMH=lTEM+NblTS 
BMHrXX(lTEMH) 

429  IK  (SHFTMH  .LT.  BMH)  GO  TO  426 
ACTM= ( BMH/MANWK ) *60 . 0 
TUIT4=TBI T4+ACTM 
ctime=ctime+actm 
TM0LE2=TM0L£2+ACTM 
THIRD2=THIRD2+ACTM 
SHFTMH=bHFTMH-BMM 

IF  (  LOtoPRT)  420 >465, 420 
485  PRINT  162 

PRINT  466 »  CTIME*  BMH>  TUNNEL>ITEM 
466  FORMAT  ( 1H  »  F10.3>5X>’  BIT  REPLACING  ’»  3X>  F10.3,  3x»  F10.3, 
*  3X,’BIT  NO,»,  13,  2X# ’REPLACED* > 

42?  CONTINUE 

SUBTRACT  STROKE  FROM  BIT  LIFE 

430  CONTINUE 

IF  (  IP  .NE.  i  )  GO  TO  670 
AA=(STROKE-TEMSTR)+TFT 
GO  TO  671 

670  AA=STR0KE 

671  DO  431  1=1#NB1TS 

IF  (MM  .NE.  0)  GO  TO  432 
AA=AA/XX(NSCF) 

GO  TO  431 
432  AA-AA 

431  XX(I)=XX(I)-AA 

REPLACE  BIT  LIFE  OF  BIT  REMOVED  WITH  A  NEW  LIFE 

IF  (KOUNT  .EQ.  0)  GO  TO  438 
DO  435  1  =  1  *  KOUNT 
ir£M=IBlT(I) 

1TEMH=ITEM4NBITS 

CALL  CAuCUM ( ITEMH , X ) 

XX( ITLWH)=X 

CALL  CALCUM(IIEM,X> 

435  XX ( ITEM ) =X 

436  IF  (  ID  .NE,  1)  GO  TO  577 
Iu=ID+l 
TEMSTR=STROKE 

GO  TO  998 
677  ICFzNSCF 


r>  n  an 
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C 

CALL  CALCUM  (  lcF#X) 

C 

XX(NSCF)=X 

AuRT=XX(NSCF) 

XX(N5CF)=XX(NSCF)/60.U 

IF  (  ICYCLE  ,£Q.  0)  GU  TO  038 

c 

CALL  MUCK 
C 

030  TTM=TFT/XX(NSCF) 
TMH=TTM*NCREW/60 , 0 
IF  (  TNL-TUNNEL)  950# 950 » 951 

951  IF  (  MAXSHT-NSHIFT)  952»952»953 

953  IF  {  TIMAX-CTIME)  954»954#987 
950  PRINT  5U7 

GO  TO  520 

952  PRINT  512 
GO  TO  520 

954  PRINT  513 
GO  TO  520 


COMPLETED  INSPECTION  AND  REPLACING  OF  BITS 
ADVANCE  BY  TFT  (CASE  0 F  TEMSTR  .GE,  TFT) 


350  CONTINUE 
TFTA=TFT 

C 

IF  (  TFT  )  771 » 361 #771 
771  WTLDG=CROSEC*TFT*GAMMA/2000, 0 
SCCTM=0.0 

IF  (  ICYCLE  ,NE.  0)  GO  TO  710 
FTAD=TFT 

C 

CALL  TRANS 
C 

GO  TO  369 

710  IF  (  LS(MSW)  , EQ.  0  ,OR#  LS(NSW) 
LWTID=1 
GO  TO  702 
701  LWTID=0 
C 


CALL  TRANS 
C 


CTIME=CTIME+WTIM 

TM0LE2=TM0LE2+WTIM 

TBIT2=TBIT2+WTIM 

THIRD2=THIRD2+WTIM 

702  IF  (  WTLDG-TEMPWT)  703,7C4#705 

703  TEMPWT=TEMPWT-WTLDG 
DO  706  I=l»NLOCO 

706  IF  (  LL( I )  .E«.  NSW)  GO  TO  707 


•EQ# 1 )  GO  TO  701 


o  o  o 


PRINT  7U8 

70\F22?2!.!  1H°»,L0ADING  WAS  ATTtMPTED  WITHOUT  EMPTY  TRAIN  AT  LOADING 

•P  rUlNI  »  J 

STOP 
707  JJ=I 

CTLOC(JJ)=CTlOC(JJ)+TnUCK 
TLOC1 ( JJ)=TLOCl ( JJ) +TMUCK 
GO  TO  369 

704  TEMPWTsO.O 
ILWTID=ILWTI0+1 
IOLOADzO 

C 

CALL  TRANS 
C 

IF  (  IDLOS  ,EQ.  1)  GO  TO  520 
C 

CALL  MUCK 
C 

TEMPWTsWTMUCK 
GO  TO  369 

705  W  TLDG-WTLDG-TEMPWT 
TEMPWT=0.0 
ILWTIU=ILWTID+1 
IDLOADrO 

C 

CALL  TRANS 
C 

IF  (  IDEOS  »E0»  1)  GO  TO  520 
C 

CALL  MUCK 
C 

TEMPWTsWTMUCK 
GO  TO  710 

CHECK  WHETHER  THE  PROJECTED  BORING  CAN  BE  DONE  IN  THE  SHIFT 

369  IF  (SHFTMH-TMH)  358»359»359 

358  SHFTMH=SHFTMH+NCREW*HRPSH 
NSHIFT=tMSHIFT+l 
GO  TO  369 

359  SHFTMH=SHFTMH-TMH 
TUNNEL=TUNNEL  +TFT 
TaiTlsT3lTl+TTM 
TM0LE1=TM0LE1+TTM 
THIRD1=THIRD1+TTM 

ctime=ctime+ttm 

C 

CALL  SUPPRT(TFTA»REQMH) 

C 

AVAMH=TTM*2,0  /60.0 
IF  (  REOMH-AVAMH)  865»865,866 
866  WTTIM=( (REQMH-AVAMH)/2,0)*60.0 


o  o  on 


9] 


C 


CTIME=CTIME+WTTIM 
SWTTIM=SWTTIM+WTTIM 
TBlT2=TdIT2+WTTlM 
TM0LE2=  TM0LE2+WTTIM 
THIRD2=THIRD2+WTTIM 
MREST=MREST  +  (NCREW-2) 
RESTMH=RESTMH+WTTIM* ( NCREW-2 } /60 
WTTMH=wrTlM*NCREW/60,0 

S94  IF  (  SHFTMH-WTTMH)  860 #862, 062 

8t>0  SHFTMH=SHFTMH+NCREW*HRPSH 
NSHIFT=NSHIFT+I 
GO  TO  894 


862  SHFTMH=iHFTMH-WTTMH 
tsuppt=tsuppt+ttm+wttim 

GO  TO  895 

865  TSUPPT=TSUPPT+TTM 
095  CONTINUE 


•  0 


CALL  CALCUM(NSCF#X) 

XX(NSCF)=X 

ADRTSXX(NSCF) 

XX ( NSCF  )  =XX  ( NSCF ) /60 , 0 

IF  (  ICYCLE  ,EQ,  0)  GO  TO  837 


C 


C 

C 

C 


CALL  MUCK 


037  IF  (LOGPRT)  361,360.361 
360  PRINT  162 


_  PRiNT  362 *  CTIMEi  TMH#  TUNNEL 
362  FORMAT  ( 1H  »F10, 3,  5X»  *  TUNNEL  B0RIN( 


DETERMINE  REPAIR  TIME 


*»3X,  F10.3»3X,F10,3) 


361  DO  370  K=1»NEVENT 
M=NBITS*2+I 
N=NSCF-1 
JCOUNTsO 
ITEM=IX(K) 

DO  365  J=M»N 
JCOUNT=UCOUNT+l 

365  IF (ITEM  ,EQ,  J)  GO  To  366 
PRINT  367 »  ITEM 

367  FORMAT  (»  EVENT  NUMBER  =  », 
*AY»  SO  RUN  STOPPED*) 

STOP 

366  IR=NSCF+JCOUNT 

370  IXR(K)=IR 


13*  f COULD  NOT  BE  FOUND  IN  CF(I,J)  ARR 


START  REPAIRS 

MAKE  AN  ARRAY  OF  IXR(I)  in  DECENDlNG  ORDER 


o  n  o  o 
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IF  (  NEVENT  ,EQ»  1)  60  TO  371 
00  372  X=1>N£VENT 
IR=IXR(I) 

372  XIXR<I)-XX(IR) 

NA=NEVENT-1 

DO  373  J=1»NA 

M=J 

MA=J+1 

DO  374  K=MA,NEVENT 

374  IF  (XIXR(K)  ,GT.  XIXR(M))  M=K 
TEMP=XIXR(J) 

XIXR(J)^XIXR(M) 

373  XIXR(M)=TEMP 
MKLrlXR(l) 

DO  840  KL=2»NEVENT 
IR-IXR (KL) 

IF  (  XX  (MKL)*»XX  ( IR) )  840 » 371 ,  a40 
840  CONTINUE 

DO  375  X=1 jNEVENT 
IR=IXR(X) 

DO  37b  J=1»NEVENT 
IF  (XIXR(J)-XX(IR) )  376#377»376 

376  CONTINUE 

377  IXRR(J)-IR 

375  IXX( J)=IR-(NSCF-NBITS*2) 

DO  645  I=1»NEVENT 
ix(i)sixxd) 

645  IXR(I)=IXRR(I) 

371  ISUMlsO 
ISUM2=0 

TO  DISTRIBUTE  MEN  TO  CREWS,  THE  BIG6ER  JOBS  ARE  ASSIGNED 
MANPOWER  FIRST 

DO  380  I=1»N£VENT 
ITEM=IX(I) 

ISUMlsISUMl+MANUTEMfl) 

380  ISUM2=ISUM2+MAN<ITEM,2) 

DO  500  KK=1»NEVENT 
IR=IXR(KK) 

ITEM=IX(KK) 

I=KK 

IF  ( IMAN-ISUM1 )  381 » 362, 383 

383  IF  (  IMAN-ISUM2)  381,384,384 
382  MANAW(I)=MAN(ITEM,1) 

GO  TO  395 

384  MANAW(I)=MAN(ITEM,2) 

GO  TO  395 

381  IF  ( IMAN-MANC ITEM, 1) )  390,391,392 
392  IF  ( I MAN-MAN ( I TEM , 2 ) )  391,391,393 
390  MREST=MKEST+IMAN 


n  o  or> 
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restmh=kestmh+shftmh 

IF  (  SMAN  «EQ.  0)  60  TO  846 
WAITM=WAITM+ (SHFTMH/IMAN) *60 . 0 
CTIME=CTIME+(SHFTMH/IMAN)*60.0 
846  NSHIFT=NSHIFT+1 

shftmh=ncrew*hrpsh 

IMAN=NCREW 

393  MANAW(I)=MAN(ITEM*2) 

60  TO  395 
391  MANAW(I)=IMAN 
60  TO  395 

395  IMAN=IMAN-MANAWCI) 

C 

C  IF  A  JOB  CANNOT  BE  HANDLED  IN  SHIFT*  IT  IS  EXTENDED  INTO 
C  THE  NEXT  SHIFT 

C 

IF  (SHFTMH  .GEt  XX ( IR) )  60  TO  396 
397  SHFTMH=SHFTMH+NCREW*HKPSH 
NSHIFT=NSHIFT+1 
IMAN=NCREW-MANAW(I) 

IF  (SHFTMH  ,LT.  XX(IR)>  60  TO  397 

396  SHFTMH=SHFTMH-XX(IR) 

WORK  fit -YK ( TR 1 

500  AC TIM l I 7= (WORK ( I ) /MANAW ( I ) ) *60*0 
C 

C  ARRANGE  ACTIM(I)  IN  ASCENDING  ORDER 
C 

IF  (NEVENT  *EQ»  1)  60  TO  550 
DO  557  IrlrNEVENT 

557  ACT(I)=ACTIM(I) 

na=nevent-i 

DO  551  J=1*NA 
M=J 
MA=J+1 

DO  552  I=MA* NEVENT 

552  IF  (ACT ( I )  .LT.  ACT(M) )  M=I 
TEMP=ACT(J) 

ACT(J)=ACT(M) 

551  ACT(M)=TEMP 

DO  558  I<=1*  NEVENT 
DO  559  1=1 » NEVENT 
IF  (ACTIM(I)-ACT(KU  559*558,559 
559  CONTINUE 

558  LOG(K)=I 
GO  TO  553 

550  ACT(1)=ACTIM(1) 

LOG(l)=l 

553  IF  (  NEVENT  ,EQ.  1)  GO  TO  630 

REDUCE  THE  VALUES  OF  ACTIM(I)  BY  REDISTRIBUTE  THE  MANPOWER 
AFTER  ACTIM(l)  IS  ACHEIVED 


noon 
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J=1 

DO  600  L=NEVENT»2#-1 
I=L 

620  M=L0G(I) 

IR=IXR(M) 

ITEM=IR-(NSCF-NBITS*2) 

IF  (MAN(ITEM»2)  .EQ.  MANAW(M))  GO  TO  471 
MADsMAN ( I TEM , 2 ) -M ANAW ( M ) 

610  IF  (MAD-MANAW(J) )  60l#602,602 
602  MTB=MANAW(J) 

IF  (  J  #NE#  1)  GO  TO  603 
DNsWORK (M) -ACT ( J) ♦MANAW (M) /60 . 0 
OD=MANAW(M)+MTB 
GO  TO  604 

t>03  DN=DN-(ACT(J)-ACT(J-15  )*DD/60.0 

do=dd+mtb 

604  ACT(I)=ACT«J)+(DN/DD)^60.0 
J=J+1 

IF  (I-J)  605#606#607 

605  PRINT  60S 

608  FORMAT  (IX# ’RUN  STOPPED  BY  AN  ERROR#  SEE  STATEMENT  605  IN  MAIN’1 
STOP 

607  mad=mad-mtb 

IF  (  I  ,NE.  NEVENT)  GO  TO  609 
IF  (  ACni)-ACT(I-l) )  600#610»610 
o09  IL=I+1 

DO  611  K=IL, NEVENT 

611  IF  (  ACT(K)  ,GT.  ACT(l))  GO  TO  612 
IF  (  ACTCI)-ACT(I-I) )  600#610#610 

612  IsK 

GO  TO  620 

o06  IF  (ACTlI)-ACT(I-l) )  6l3»471#47l 

613  PRINT  614 

614  FORMAT  (IX# ’RUN  STOPPED  BY  AN  ERROR#  SEE  STATEMENT  606  IN  MAIN’) 
STOP 

THE  CASE  OF  MAD  *GE.  MANAW (J)  IS  COMPLETED 
THE  CASE  OF  MAD  .LT.  MANAW (J)  IS  BEGUN 

601  MTB=MAD 

MANAW  ( J)=MANAW  ( J)  -MTB 
IF  (J  ,NE.  1)  GO  TO  621 
DN=WORK ( M ) -AC  T ( J ) ♦MANAW ( M ) /60 . 0 
DD=MANAW(M)+MTB 
GO  TO  622 

621  DN=DN-(ACT(J)-ACT(J-m*DD/60.0 
DD=OD+MTB 

622  ACT(I)=ACT(J)+(DN/DD)*60.0 

IF  (  I  .NE.  NEVENT)  GO  TO  623 
IF  (  ACT ( I ) -ACT ( 1-1 ) )  600»471»471 
e23  IL=I+1 

DO  624  K=IL#NEVENT 

624  IF  (  ACT(K)  #GT,  ACT(D)  GO  TO  625 


o  o  o  n 
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IF  (  I-(J+1) )  626,627,628 
t>26  PRINT  629 

o29  FORMAT  <1X»»RUN  STOPPED  BY  AN  ERROR,  SEE  STATEMENT  626  IN  MAIN»1 
STOP 

b27  GO  TO  471 

028  IF  (  ACr(I)-ACT(I-l) )  600,471,471 
o25  I=K 

GO  TO  620 
OOO  CONTINUE 

COMPLETE  THE  REDUCTION  OF  THE  ACTIM(I)  VALUES 
REARRANGE  ACT(I)  IN  ASCENDING  ORDER 

471  CONTINUE 

DO  631  L=1»NEVENT 
031  JON ( L ) =0 
IM=1 

DO  632  K=l» NEVENT 
DO  633  IK=1, NEVENT 
MM=JON(IK) 

033  IF  (IM  .EG.  MM)  IM=IM+1 
M=IM 

DO  634  1  =  1 » NEVENT 
DO  636  N=l» NEVENT 
JJ-JON (N) 

635  IF  (  JJ  .EG.  1)  GO  TO  634 

IF  (ACT(I)-ACT(M) )  636,634,634 

636  M=I 

o34  CONTINUE 

CTM(K)=ACT(M) 

JON(K)=M 
032  CONTINUE 
GO  TO  47? 

030  CTM(1)=ACT(1) 

J0N(1)=L0G(1) 

472  IF  (  LOGPRT)  473,470,473 
470  PRINT  162 

473  DO  455  K=l» NEVENT 
IK=JON (K) 

M=LOG(IK) 

IRSIXR(M) 

J=IR-NSCF 

TCT=CTIME+CTM(K) 

RJOB=WORK(M) 

GO  TO  (456,457,858) »J 

456  TBIT3=Ti3IT3+CTM(K) 

IF  (  LOGPRT)  455 » 275, 455 
275  PRINT  459,  TCI »  RJQB,  TUNNEL 

459  FORMAT  <1H  ,  F10.3,  5X,»  BIT  REPAIR  »,3X,  F10,3,3X,F10.3) 

GO  TO  455 

457  TM0LE3=TM0LE3*CTM<K) 


conn 


9  (> 


IF  (  LOGPRT)  455*276,455 
276  PRINT  460*  TcT»  RJOD#  TUNNEL 
460  FORMAT  (1H  »  FlO»3»5X»f  MOLt.  REPAIR 
GO  TO  455 


tfbH  THlRD3=rHlRD3+CTM(K) 

IF  (  LOGPRT)  455*277,455 
277  PRINT  461*  TCT*  RJOB,  TUNNEL 
461  FORMAT  (1H  »F10.3,5X,»  THIRD  REPAIR 


*♦55  CONTINUE 


•  »3X,F10.3»3X»H0.3) 


♦,3X*F10.3»3X»F10.3) 


DETERMINE  CLOCK  TIME  AND  WAITING  TIME 

CTIME=CTIME+CTM(NEVENT  > 

IF  (NEVENT  .EQ»  D  60  TO  404 

NA=NEVENT-1 

BIG=CTM(NEVENT> 

DO  400  I=1»NA 
WTM=BIG-CTM(I) 

IK=JON(I) 

M=LOG(IK) 

MREST=MREST+MAN AW ( M ) 
RESTMH=RESTMH+(WTM*MANAW(M) ) /b0,0 
IR-IXR (M) 

J=IR-NSCF 

GO  TO  (401*40^*403) »  J 
401  TM0LE2rTM0LE2+WTM 
THIRD2=THIRD2+WTM 

GO  TO  400 

*+02  TbIT2=TDlT2+WTM 

THIRD2=THIRD2+WTM 
GO  TO  400 

403  TUIT2=TBIT24W1M 

TM0LE2= TMOLE2+WTM 


40C  CONTINUt 
GO  TO  405 

404  MKEST=MREST+ ( NCREW-MANAW ( 1 ) > 

RESTMHrHESTMH ♦ ( ( NCREW-MANAW ( 1 ) ) *CTM ( 1 > ) /60 . 0 


IR-IXR ( 1 ) 

J=IR-NSCF 

GO  TO  (406*407*408) »J 
406  TM0LE2=rM0LE2+CTM(l) 
THIR02=THIRD2+CTM(1) 


GO  TO  405 

407  TBIT2“T6IT2+CTM( 1 ) 

THIRD2=THIRD2+CTM ( 1 ) 


GO  TO  405 

406  TBIT2-T6IT2+CTM( 1 ) 
TM0LE2=TM0LE24CTM ( 1 ) 


SUBTRACT  XX (ITEM)  FROM  XX(I) 
END  OF  REPAIRS 


o  o  r>  r  o  o  r,  r  nn  ooo 


*♦05  ITEM=iX(l) 

M=NBITS*2+1 
N=NSCF-1 
DO  490  I=M#N 
490  XX(1)=XX(I)-XX(ITEM) 

REPLACE  XX ( ITEM)  WITH  NEW  XX(lTEM) 

DO  495  J=1»N£VENT 
IT£M= IX ( J) 

CALL  CALCUM(ITEM#X) 

XX ( ITEM) sX 
IR=IX«(J) 

CALL  CALCUM(IR,X) 

495  XX(IR)=X 

PERFORM  A  BIT  INSPECTION  IF  THE  MOLE  IS  DOWN  FOR  OTHER  REPAIRS 

DO  555  I=1#NEVENT 
IR=IXR(I) 

IDR=IR-NSCF 

555  IF  (  IDR  .EG.  2)  GO  TO  666 
GO  TO  6o7 
ob6  ID=1 

IF  (  ID  .EQ.  1)  GO  TO  481 
667  TEMSTR=rEMSTR-TFT 
998  CONTINUE 

COMPLETED  THE  CASE  OF  TEMSTR  .GT.  TFT 
CHECK  TERMINATION  VARIABLES 

IF  (  TNL-TUNNEL)505»505»506 

506  IF(MAXSHT-NSHIFT)  510»510»515 
515  IF  (  TIMAX-CTIME)  5il»511»999 
505  PRINT  507 

507  FORMAT ( 1H1 # ’SIMULATION  WAS  TERMINATED  BY  THE  MAXIMUM  ADVANCE  OF  TH 
■SE  TUNNEL*//) 

GO  TO  520 
olO  PRINT  512 

512  FORMAT  ( 1H1 »» SIMULATION  WAS  TeRMINATEDbY  THE  MAXIMUM  NUMBER  OF  SHI 
*FTS*//) 

GO  TO  520 
511  PRINT  513 

513  FORMAT  (  lHlf 'SIMULATION  WAS  TERMINATED  BY  THE  MAXIMUM  CLOCK  TIME 
*»//) 

o20  PRINT  521 

o21  FORMAT (5X» * ♦♦♦♦♦♦♦♦♦♦SIMULATION  SUMMARY  DATA***********//) 


•■jaw 


o  c.  o 


98 


PRINT  THE  SUMMARY  OF  I  HE  SIMULATION 

CTIME=CTIME/bO.O 

WAITM=WAITM/60.0 

TUIT5=Tt3lT5/60.0 

TdITl=TdlTl/60.0 

TdIT2=TbIT2/60.0 

TciIT3=TUIT3/60.0 

TbIT4=TbIT4/60.0 

TMOLE1=TMOLE1/60.0 

TMOLE2=TMOLE2/60.0 

TMQLE3=rMOLE3/60.0 

ThIRDl=THlRDl/60.0 

THIR02=THIRD2/60.0 

ThlR03=rHlRD3/60.0 

TSUPPT=TSUPPT/60.0 

SwTTIM=SWTTIM/60.0 

PRINT  525 

325  FORMAT  (1X» ’CLOCK  TIME* »3X» 'NO.  SHIFT' , 3X» ‘TUNNEL  LENGTH' »3X» 

*  »MH  ON  REST* »3X» 'IDLt  TIME »» 2X » 'BIT  INSP  HR*  //) 

PRINT  526’  CTIME»NSHIFT»TUNNEi_»RESTMH»WAITM, TBITS 

326  F0RMAT(lX»Fl0.3,6X»l5»4X»Fl0.O’6X»Fl0.3»2X»F9.2’2X»Fl0.3  ///) 
PRINT  527 

S27  FORMAT  (l3Xr»hR  WORKED ' »3X» 'HR  WAITED' »3X» 'DOWN  T1ME'»3X» 

*  'REPLACING  HR'  //) 

PRINT  528»  TBIT1»THIT2»T0IT3» T6IT4 
326  FORMAT  (  10H  BIT  »  3F12.3’  3X,  F12.3) 

PRINT  530’  TMOLE I » TM0LE2 » TM0LE3 
S3C  FORMAT  (  10H  MOLE  ,  3F12.3) 

PRINT  531 » THIRD 1 » THlRu2»  THIRD3 
331  FORMAT ( 10H  THIRD  »  3F12.3) 

PRINT  097 »  TSUPPT,SWTTIM 

697  FORMAT  ('  TOTAL  SUPPORTING  TI.v|E=*,  F12.3/'  »»♦  TOTAL  TIME  DELAY! 
*BY  THE  SLOWER  SUPPORTING  SYSTEM*,  F12.3) 

I0E0S=1 

CALL  TRANS 

STOP 

END 


COMPILATION 


NO  DIAGNOSTICS. 


o  c<  o  c. 


99 


SUBROUTINE  CAECUM < ICF. X) 

common  wtmuck»wtim»nloco»tmuck»kmax.ncars»lclas,ndc,accmax, 
't>VELMAX»UECEE»NSC  »D1 » AFT .HAUL .NSECS. NSW, WTCAR. FCAR. TPM. ADRT , 
*>CROSEC»  GAMMA,  CTIME.LWT  10  »TTM,TNL»TlMAX,NRBG» ILS, 
*CFR(13).FTA(l3).NCLAS(120).LUl0),CTLOC(10),TLOCl(t0).«TLOC(10). 
»FLOCOllo) »G<100) .0(100) »DS(lO) >WTLOAD(lO) .SPEED (10) » LOAD (10) . 
1>DISTRllO)  .TIME (10)  .TSfOPdO).  1  PASS (10)  »WTTRN(lO).CFE(13)  ,CFD(13)  . 
%C1  (13) ,WTL(l3) »OSTOP(iO) >ST(13) .CFS<13) .TLOAD(IO) »TDUMP(lO) , 

*  CF(120»13) »TV(l2C>13) »T(10.3G) .S(lO»30) »LS(20) .1EWTID 
*»  1DL0aD»SCCTM.T0IT2»  TM0LE2»THIRD2»dELTh.  ideos.tnlt 
t.FTAD.TMOLEl ,  TBIT1 .TH1RD1 .NCRtW. SHFTMH.HRPSH. NSHIFT » I  MAN.  ICYCLE 


SUBROUTINE  CALCUM  DETERMINES  AN  ABSCISSA  VALUE  ON  A  CUMULATIVE 
PROBABILITY  CURVE  CORRESPONDING  TO  A  GIVEN  RANDOM  NUMBER 
BETWEEN  0.0  AND  1.0 

I  =  ICF 
YsttAND(N) 

JN=NCLAS(I) 

DO  100  J=1»JN 

IF  (CF(I.J)  .LE.  Y  .AND.  CF(I.U+1)  ,GE,  Y)  GO  TO  102 
lOO  CONTINUE 

PRINT  10  3.  I.Y.JN 

1C3  FORMAT  (IX. *1  Y  JN  AT  103  IN  CALCUM’.  15.  F10.8.  110) 

PRINT  111 

111  FORMAT  (IX. 'RUN  TERMINATED  BY  ERROR  IN  SEARCHING  THROUuH  THE  CF(I, 
S>J)  * ) 

IF  (  JM  , GE,  100  )  GO  TO  113 
PRINT  112  (  CF(I.J).  o=l » JN) 

112  FORMAT  (  3F18.8) 

113  STOP 

102  X=TV(I,J)+(TV(I.J+1>-TV(I»U))*((Y-CF(I.J) )/(CF (I. J+1)-CF (I. J) ) ) 

return 

END 


compilation: 


no  DIAGNOSTICS, 


o  r>  o  r. 


100 


SUBROUTINE:  SUPPRT  (  TFTA»REQMH  ) 

COMMON  wtmuck»wtim»nluco»tmuck*kmax»ncars»lclas»ndc»accmax, 
t>Vt  LMAX  »l)FICEL»NSC  »D1  »AFT »HAUL*nSECS»NSW»WTCAR»FCAR»TPM» ADRT , 
bCROSfc.C»OAMMArCTlMh  »LW i ID»TTM» TNL*T lMAX,NRRG» ILS» 
bCFimj)  »FTA<13)  »NCLASU20)  »LL(iO)  »CTL0C(10)  » TL0C1 ( 10 ) r *TlOC ( 10 ) # 
tFuOCO(lu) »G(100) »D(100) »DS(10) »WTL0AD(10) *  SPEED (10) »LOA()(10) » 
*DISTR(10)  »TlML  ( 10 ) » TS TOP ( 10 )  »  TPASS  ( 10) » WTTRN (10)  »CFL  (13)  »CFD  (13)  » 
*CT ( 13) ,«TL<13) .OSTOP(IQ) #ST( 13) » CFS ( 13) , TLOAD ( 10) » TDUMP ( 10 ) , 

*  CF(12C»1J) »TV(12C»13) »T(10»3u) »S(l0»30) rLS(2Q) ,ILWTI0 
* » I OLOAP * SCCTM » T  J I  T  2 » TMOLE2 » THiRD2 » DELTH » I DEOS » TN|_T 
i»FTAD»TMOl.El»  TBIT1 » THIRD  1  »NCRtW»SHFTMH»HRPSH»NSHIFT » IMAN»  I  CYCLE 


SUBROUTINE  SUPPRT  UETlRMINES  THE  MANHOURS  OF  TIME  REQUIRED  TO 
COMPLETE  THE  SUPPORT  WORK  FOR  ONE  UNIT  OF  ADVANCE  OF  THE  TUNNEL 

SPMH=0 . u 
REQMHsO • 0 
NFTR=INT(TFTA+0.4) 

IF  (  NFTR  .EQ.  0)  00  TO  100 
DO  10  1=1 » NFTR 
Y=RANU ( N) 

DO  20  J=1»NRUG 

20  IK  (  CFK(J)  ,LE.  Y  .AND.  CFR(J+1)  ,GE,  Y  )  GO  TO  30 
PRINT  21 

21  FORMAT (1X»  *CFK( I)-FTA(I)  UlAGnAM  HaS  INCORRECT  INPUT  DATA’) 

STOP 

3r  5PMH=FTA( J)+(FTA(U+1)-FTA(U) )*( (Y-CFR(J) ) / ( CFR ( J+l )-CFR ( J) ) ) 

in  REQMH=REQMH+SPMH 

IOC  RETURN 
END 

imQ 


COMPILATION 


DIAGNOSTICS. 


n  r>  r>  r;  r> 


I  0  I 


SUBROUTINE  MOI ION(NL#NS) 

common  wtmuck#wtim#nluco#tmuck»kmax»ncars»lclas»noc»accmax, 

■&VELMAX  #DFCE'L#NSC  »Ul » AFT  #  HAUL » NSEC S » NSW » WTC AR  »FCAK#  TPM#  ADR T » 
*CROSEC#OAHMA#CTlMt;#LWlID»TTM,lNL»TlMAX,NRBG»lLS» 

1>CFR ( 13)  #FTA(13)  #NCLAS(120)  »LLU0)  »CTLOC(10)  »TlOC1(10)  »WTLOC(10)  » 
■fcFl.OCO ( 10)  #G(100)  #0(100) #DS(10) *WTLOAD(10)  # SPEED ( 10)  #LOAD(  10 )  » 
*OISTR(lu) » TIME (10) #TSTGP(10) »  fPASS ( 10 ) » WTTRN ( 10 ) »CFL ( 10) , CFD ( 13) . 
‘frCT (13) ,WTL(13) #0STOP(10) »ST(1j) »CFS(13)»TLOAD(10) »TDUMR(lO) , 

*  CF(120#13) # TV (120 #13) » T ( 1 0 » 3u ) »S ( 10 » 3q ) ,LS (20) # ILWTID 

* » I ULO AD » SC  CTM » TB I T 2 » TwOLEg #  TH i RD2 #  DELTh  » IDEOS » TNLT 

4>»F  TAD»TM0LE1,T[)IT1»TH1RD1»NCRL^»SHFTMH,HRPSH#NSHIFT»IMA|M,ICYCLE 

SUBROUTINE  MOTION  MOVES  A  TRAIN  FROM  ONE  SWITCH  TO  THE  NEXT 

DIMENSION  SGL(IO) #MSS(1C) #GREV(100) »DREV(100) »DSREV(1U) 

REVERSE  PERTINENT  VARIABLES  FuR  LOADED  TRAINS 

IF  (  LOAD(NL)  .EQ.  0)  GO  TO  luO 
NNN=C 

IF  (  D(NSECS)-DELTH)  til5»816»Bl6 
61b  NSECSsNSECS-1 
NNNsl 

616  DO  101  I=1»NSECS 
IREV=NSECS-(I-1) 

GREV(IREV)=-G(I) 

101  DREV(IREV)=D(1) 

DO  99  I=1»NSECS 
G(I)=GREV(I) 

99  D(I)=UREV(I) 

DO  102  1=1#NSW 
J=NSW-(I-1) 

102  DSREV(I)=HAUL-DS(J) 

DC  98  I=1»NSW 

98  DS(I)=DSREV(I) 

ICC  SSTzO.O 

S6L ( 1 ) =D ( 1 ) 

MSS ( 1 ) =1 

SSL (NSW) =0.0 

MSS(NSW)=0 

N9=NSW-1 

DO  103  I=2»N9 

00  104  J=l»NStCS 

SST=SST+0(J) 

If  (SST-OS(I))  1 34 » lQo# 107 

106  SOL( I ) zD ( J+l ) 

MSS(I)=J+1 

GO  TO  lo3 

107  SGL(I)=SST-DS(I) 

MSS (I )=J 


102 


GO  TO  lu3 
104  CONTINUE 
103  CONTINUE 

T1ME(NL>=0.0 

TSTOP(ML)=0.0 

TPASS(Ni.)=O.0 

l)SrOPCNL)=0.0 

START  LOCO  IN  MOTION 

IF  (  LOAD(NL)  .EQ.  0)  GO  TO  80 
JK=NSw-(NS-l) 

NS=JR 

OF  SLEFT=nb(NS+l)-OS(NS) 

GLEFT=SGL(NS) 

MM-MSS (NS ) 

FKFC=(  FLOCO(NL)*WTLOC(NL)+(  wTLOAD(NL)+WTCAR*NCARS)*FCAR) 

IK  (  LOAO(NL)  .EQ.  1)  GO  TO  7/7 
WTTRN(NL)=WTLOC(NL)+  ImCARS+WToAR 
m  GFC=(  2U.0*G(MM) )*WTTRN(NL) 

RE0TF=FUFC+5FL 

IF  (  SPEED(NL)  )  110,110*111 
110  AVATF=(T(NL,1)+T(NL#2) )/2,0 
120  ACCFC=AVATF-R£QTF 

ACCR=ACCFC*32.2/ ( WTTRw (NL) *2000 , 0 ) 

IF  (  ACCR  )  115, lib* 117 
US  TL)=0.0 
00=0.0 

IK  (  SPEED(NLJ)  130*130*131 
13C  IF  (LOAU(NL)  ,EO.  1)  GO  TO  132 

PRINT  133*  NL»MM»G(MM)»ACCR»WTLOC(NL> *WTLOAD(NL) ,WTCAR*NCARS»FCAR, 
SWTTRN(NL) »FLOCO(NL) » FKFC »GFC *kEGTF , AVAtF* ACCFC » T (NL* 1 > » T (NL * 2) 

133  FORMAT  (  1HO,*LOCO  N0.=»»  13,  ’WITHOUT  LOAD  CANNOT  NEGOTIATE  SEC f I 
$011  NUMRcR  '*14  /*  ',3A,*GRADE=’*F10.4»3X,,ACCR=’,F12.4  /»  », 

%  3F1Q.2*  15,  3F10.2  /»  *,7F12.3) 

STOP 

132  PRINT  134,  NL , MM , G ( MM ) , ACCR , WTLOC ( NL ) * WTLOAD ( NL ) , WTCAR * NCARS ,FC AR , 
•t.  WTTRN(NL)  »FL0C0(NL)  ,krfc,gfc*.<eotf*aVatf*accfc,t<nl»i>  ,T(NL»2) 

134  FORMAT  (  1HO,'LOCO=’»  13* ’WITH  LOAD  CANNOT  NEGOTIATE  SECTION  NUMBE 
*R  • , 14  /»  *,3X*' REVERSED  GRADt=’ *F10.4,3X* ' ACCR=* *F12.4  /»  *, 

*  3F10.2*  15,  3F10.2  /’  ',  7F12.3) 

STOP 

131  IF  (  GLEFT-SLEFT)  135*135,130 

135  IF  (  SPEED ( NL ) “S ( NL » 2 ) )  136,136,137 

136  ADEC=-ACCR 
TTD=SPEED(NL)/ADEC 

TUD=SPEED ( NL ) *SPEEU ( NL ) / ( 2 . 0*mOE  C ) 

TD=TDfTTD  * 

DO=TOD+OD 

IF  (  DD-GLEFT)  130,130*138 
138  DTR=DD-GLEFT 
TD=TD-TTD 


I ! 


0 


103 


L)J=no-TUD 

sp=sPhe:u(NL) 

ASArSPEED (NL)  *SPEED  (Nl.)  -2 ,  Q*Al)£C*DTT< 

IF  <  ASA  )  261  *  261 , 26*: 

«iol  PRINT  2o3#ASA 

s?opA'  llX>,«RI«L£  ASA  IN  MuTIOH  HAS  NEGATIVE  SIGN* »  fio.j, 
dt>2  CONTINUE 

SPZED  ( ML )  =SQRT  ( SPEED  ( .ML )  ♦SPEED  ( NL )  -2 •0*ADEC*DTR ) 
TT1=(SP-SPEED(NL))/ADlC  AUtCKJTK) 

TDlrOIR 
D1=D0+TL)1 
T1=TD+TT1 
Gu  TO  161 
107  AuEC=-ACCR 
SP=SPEFU(NL) 

N7-KMAX-1 


b't 

139 

i4i 


14-: 


142 


DO  I3y  l=l»N7 

IP  (  bPcED(NL)-S(NL»I))  139,1j9,60 
If  (  SPEED ( NL ) -S ( NL » I  + 1 ) )  140.l4G.i39 

CONTINUE 

PKINT  1<+1,  SRtEO(NL),,ML 

FORMAT  (  1H0, •SPEEDS*, F10. 3, *0F  L0C0  NUMBER*. 13 
»D  IN  ITi  CHARACTERISTIC  CURVE*) 

STOP 

KK-I 

T T  (  SPEED ( NL ) -S ( NL , KK ) ) /ADEC 

TU=TO+TT 0° ( NL '  *SP!:E0  <  (>“- )  -S  ( NL ,  KK )  ♦S  ( NL ,  KK  >  )  /  ( 2 . 0* ADEC ) 
DU-OD+TUO 

IP  (DD-GLEFT)  142.138.13B 
SP-SPLFj(NL) 


'COULD  NOT  bE  FOUR 


144 
09  r- 
901 


143 


M  S  (NL^KK+1  J'*S(NlIkk  )  J  **  * SPEEO (NL) -S (NL* KK ) )/ 

SPF.Ed  ( NL )  =S  (Nl  » KK  ) 

AVATF=(AVA1+T(NL.KK) )/2.0 

ACCFC=AVATF-RlQTF 

ACCR=aCCFC*32.2/(WTTRn(NL)»20uJ.0) 

IF  (  ACCR)  136.143. 144 

IP  {  ACCMAX-ACCR)  990.991,991 

ACCRsACCPAX 

SPEED ( NL ) =S ( NL , KK ) +ACCR* ( S ( NL , KK  +  1 ) -S ( NL  ,  KK ) ) / ( AOfcT+ACCu i 
T J lz ( SPEED (NL)-S(NLrKK) )/ADEC  K) )/(ADtC+ACCK) 

TOATUt?alU<NL,*SPKl;D(l',L,‘S<N1-*KKI*S<NL*KI<"/(£*0*ADEC) 

dd=dd+ddi 

TD2= { GLEFT-DD ) /SPEED ( nL ) 

DU2= (6LEFT-DD) 

Ol=DD+nu2 
T1=TD+TD2 
GO  TO  lal 


o  r  r:  r  n  r 


I  O  ', 


CALCULATE  THl-  CASE  OF  2ER0  ACCELERATION 

lit'  IF  (  SPEED (NL)  )!3C»l3o*14b 
145  J1=GLEFT 

Tl=GLLFr/SP£EU(NL) 

SP=SPEED(NL) 

GO  TO  tbl 


CALCULATE  the  case  of  positive  acceleration 


117  IF  (  aCCR-ACCMAX)  146,146,147 

147  ACCR=ACCMAX 
146  N7=KMAX-1 

DO  146  1=1 »N7 

IF  (  SPEED(NL)-S(NL»I) )  143 » 5i » 51 
bl  IF  (  SPeFU(NL)-S(NL,I+1) )  150»14Q,148 

148  CONTINUE 

PRINT  l49»NL,aPEEfJ(NL) 

149  FORMAT  (  1HO,»LOCO  NUMBER * , I3» ’  HAD  BAp  INPUT  DATA  OF  SPEED  -»,FlO 
1,2) 

STOP 
IbC  KK=I 

T 1= ( S ( ML » KK+ 1 ) -SPEED ( HE ) )/ACCK 
D1  =  T1+  (SPEED (NL)  +  (Tl*  (ACCR/2,*,)  )  ) 

D2=GLEFT 

V V=SPEE J ( NL ) ♦SPEED ( NL) +2 . 0*ACtR*D2 
T2= ( SoR T ( VV ) -SPEED (NL) )/ACCR 
IF  (  ul-D2)  bbl » 152 , is2 
lb2  D1=D2 
T1=T2 

SP=SPEED(NL) 

SPEED ( NL ) =SP+  r 1 *ACCR 
GO  TO  IdI 
obl  SP=SPEED(NL) 

SPEED (NL)=S(NL,KK+1) 

1 b 1  DST OP ( NL ) =SPEEO ( NL ) *SPCED ( NL ) / ( 2 . 0+DECEL ) 
lb3  DISTR(NE)=DISTR(NL)+D1 
T1ME(NL)=TIME(NL)+T1 

gleft=gleft-oi 

SLEFT=SLEF T-oi 

IF  (  SLLFT-O.b  )  lb4, 154,607 
b07  IF  (  GLEFT-0.5  )  156,156»157 
lb 6  MM=MM+1 

GLEFT=D(MN|) 

GFC=(20.0*G(MM) ) *WTTRN(NL) 

rlotf=fkfc+gfc 
ill  continue 

lb7  IF  (  SLtFT-GLEFT)  50l»5Cl,500 
PCI  IF  (  US TOP (NL) -SLEPT)  502,503»154 
POP  IF  (  SPEED (NL)-VELMAX)  506 , 50b » 505 
P05  D1=SLEFT-DST0P(ML) 

SPEED <NL)=VELMAX 


105 


T1-D1/SPEED(NL> 
nME(NI.)=TlME(NL)  +  rl 
DISTR(ML)=Dlsi  R  (ND+Dl 
Tb  T  OR ( r  jl ) =SPCLD ( NL ) /OLCEL 
TRASS ( ML ) =OSTOP ( NL ) /SPEED ( NL ) 

GO  TO  4u5 
ace*  gllft=sleft 

Go  TO  SUO 
a  00  CONTINUE 

IF  (  SPlED(NL)-VELMAX)  15H,21o»216 
^16  SPEED  (NL)=VELMX 
GO  TO  U6 

ibP  no  159  K=1 r KMAX 
lo  Jf  <  SPEED  (NL)  -S  (NL»K  i )  159,52*52 
0t:  IF  (  5PEFD(NL)-S(NL,K+1) )  160,159,159 
ISO  CONTINUE 

PHIMT  lol ,  NL?  SPEED ( NL ) » ( S ( NL  *  K ) ,  Kzi.KMAX) 

lblZT^'S?0°P&°y  *T  SPE“  W-fW  BA,  INPUT  0 

STOP 
Xoi'  KK*"K 

M  u?nl"kk?1  I  -S !  mJ.KKL> ;  ^K+1  >  -T  <  '>1-  KK  "  *  <  SPEED  ( ML )  -b ,  NL ,  KK ,)  / 

AVATF=(AVAl+T(fJL,KK  +  l)  )/2,0 
GO  TO  120 
ISA  SLLFT^pi+SLEFT 

OTRO=(Di+DSTOP(NL) )-SuEFT 
IF  (  ACCR  )  400,401*400 
401  Di=ni-orRD 

DI5TR ( Nu ) =D I S TR ( NL ) -0  T  RD 
TIME  (Nu)=TIME(NL)-T1 
Tl-Dl/SPEED (NL) 

TIME(NL)=TIME(NL)+T1 
TSTOP ( ML ) =SPEL0 ( ML ) /DECEL 
Tf 'ASS  ( ML )  -DSTOP  ( ML )  /SPEED  ( NL ) 

GO  TO  t|  u5 

1+01  V2=(2.0*SLEFT+SP*SP/alCR)*(AClR*DECEL/(aCcR+DECEl)  ) 

rM1/11  t,A3r  V2, SP»  SPEED (NL )  » ACCR  *  SLEFT  *D1  ,USTOP (Nl) 

*+^*  Ii^TFIO^  ***r,V,c;  MOTION  HAS  NEGATIVE  VALUE  V.  *,2X» 

STOP 

412  SPEFD(ML)=SQRr (V2) 

DISTR(ML)=DISTR(ML)-D1 
TIME  (Ml)=TIME(ML)-T1 

DST0P ( ML ) =SPelD ( NL ) *SPEEO ( NL ) / ( 2 . 0*DECEL ) 

Ul= ( SPEED ( NL ) *SPEE0 ( NL ) -SP*SP ) / ( 2 . 0*ACcR ) 

DISTR(NL)=DlsrR(NU+Di 

Tl= ( SPEED ( NL ) -SP ) /ACCR 
TIME(NL)=TIME(ML)+T1 
1S5  TSTOP (NL)=SPEED(NL) /DECEL 

TP ASS ( Mu ) =DST0P ( NL ) /SPEED ( NL ) 

405  IF  (  LOAD i NL )  ,EQ.  0)  GO  TO  4t>5 


O  Cl 
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NS=NSW-NS+1 
00  861  l  =  l*NSt-CS 
iuev=Nc:£CS-(i-i) 
GKEV{IRLV)=-G(I> 

661  DKeV(iP£V)=0(I) 

L)0  062  l  =  l»NStCS 
G ( I ) =GRLV ( 1 ) 

662  n(i)=untv(i) 

DO  065  l=l»N5h 
J=U5lN-(l-l) 

ooj  D6KEV(I)=HAUL-0S(J) 

DO  86h  1=1 » NS a 
664  US(I)=DSREV(I) 

IF  (  nNN  .EO.  0)  GO  TO  455 
NSECS=MbECS+l 
*456  Rt-TURN 
END 

COMPILATION:  no  diagnostics. 


COMMON1  i  TMUCkC  WT I M » NLOC0 .  TMUCK  •  KMAX » NC, ARSp  .R^C  i  ^Spjf  ADRT ' 
SVLLMAX .  DECEL » NSC  r Di  t  AFT » HAUL » uSF.CS » NSrt. WTC AR » FC  AN »  TPM » ADRT » 

t2  c  1ST!  "SoTTELI  2o»Tc?Kc7i2»  Ixtoei .  10  > . -TL.oc « 10 » . 

..  n  ah  .  SCCTM » TB I T2 » TM0LE2 » TH1RD2 1  DELTh  » I DEOS » TNL I 

•i.!rrAD^TM0LEl,Tnm»THiR01»NCRLW»SHFTMH,HRPSH»NSHlFTf  IMAN»ICYCLE 

SUIJROUT1NE  LOADNS  SIMULATES  THE  LOADING  OF  THE  TRAIN 


TLOAD(ML)=O,0 

WTLOAj{NL)=0.0 

Y=RAWU(N) 

NC'-NSC- 1 

60  IF  ^°CPS(K)6rLT.  Y  .AND.  CFS(K*1>  .GE.  Y  >  GO  TO  61 

62  FORMAT6?  2X,. STOPPED  oY  AN  ERROR  IN  SEARCHING  THROUGH  CFS(I)M 

61  TS^=ST(n)+(ST(K+1)-ST(K))*{(Y-CFS(K))/(CFS(K+1)-CFS(K))) 

TLOAD (ML) =TMUCK+TSfl' 

WTLOAu(NL)=WTMUCK 

WTTRN{ML)=WTLOAD(NL)+NCARS*WrcAR+WTLOC(NL) 

LOAD(NL)=l 

RETURN 

END 


compilation: 


no  diagnostics. 
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SUBROUTINE  DUMP  SIMULATES  THE  DUMPING  OF  THE  TRAIN 
Y-RANU(N) 

N5=NCJI>1 

100  IF  (°cFoli)  ?LT,  Y  .AND.  CFD(J+1)  .GE.  Y)  GO  TO  110 

120  FORMATlflHo!»lNCnuMPlNG  CYCLE  RANDs’ >FlO.€>*  2X. ’CANNOT  bE  FOUND’ 
b/*  *»  3F10.6) 

110  TDUMP(NL)=CTCJ)  +  CCT(J+l)-CT(J)»*((Y-CFDCJn/CCFD(J+l)-CFD(J))) 

L0AD(NL)=0 

RETURN 

END 


compilation: 


diagnostics, 


II 
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SUBRGUT INE  MUCK 

COMMON  ^T,MUCK»WTlM,NLUCO,TMUCK»KHAX»NCAHS»LCLAbrlMUC,ACCMAX, 

' UFLEL  '  W!*C » D1 » AH’ t HAUL , rtSECS,  NSW ,  WTC AR , FC AK » TPM » ADKl . 
iCKOSEC»oAMMA»CTTMt'»LWI  ID#TTM,  TNL»TlMAX,NRUG»  IlS» 

*CH<(13)  »FTA(1J)  »NCLAS(12Q)  »LL  ( 10 )  » CTLOC  ( 10 )  »TlOC1(10)  rwTLOCf 10)  * 
^FLOCOdo)  »G(UO)»0(10U)»DS(10)  »WTL0AD(10)  > SPEED (1U)  » LOmu ( 10 ) 
oil!  iTR  ( 1  o )  » T I  ML  CIO)  » TS  TOP  ( 10 ) »  I  PASS  ( 1 0  )  ,  WTTRN  ( 10 )  »CFL  (1 3 »  ,CFnil  *  i 
,WTL<13) .DSTOPUQ) ,ST< 1J) .CFSClsJ “iLOAotlOl »TDui?(Io)!f (13> ' 

*  ^1(,1w2,JfrTTVi1?°'l3)  'T(10»  ■J^)  *  S  ( 1 0  r  3o  )  »LS(20)  rILWTID 

*  * 1 GLOAC » see  Tm  t  TP  I T2 1 TM0LE2  *  Ti  l  i  RD2  ?  qELTh  » I OFOS , TNL  T 

5*»F  FAD, 7MOLE1,  TBIT1  »TH1RD1#MCRl.W»SHfTMH,HRPSH»NSH1FT»  iCYCLE 


ScarsJTINE  MUCK  ncrER,'IINES  Tl1^  ^ight  of  the  muck  loaded  into  the 

AFTzO.O 

TMUCK=0.0 

WTMUCKzu.O 

N3=LCLAS-1 

DO  10  1=1 »NCAKS 

Y=RANu(n) 

DO  20  U=1»N3 

t;1  ;LT>  V  •AND*  cFUJ+i)  .ge.  r  ,  so  to  3u 

‘‘■‘‘HUwV*?;:,1'1  T"E  L0ADING  CYCLE  RANDI''  F10.6.2X, -CANNOT  BE  FOB 

PRINT  2b  (  CFu(K),  K=1»LCLAS) 

26  FORMAT  (  F12.&) 

STOP 

OcJ  WTMM=WTL ( J)  +  (  WTL  ( J+l ) -WTL  ( J)  )  *  (  ( Y-CFL ( J)  ) / (CFL  f J+1  »-rPi  , 

1°  WTMUChrWTMUCK+WTMM  L  1  lChL(J+1>  CFL(J))) 

TP1 1= ADR T*CROSEC ♦GAMMA/ (60,0*2oQ0#G» 
tmuck=wtmuck/tpm 

AFT=AuRT*TMUCK/63,0 

return 

End 


compilation: 


NO  DIAGNOSTICS, 


c  p  c,  Ci  fi  o  o  n  r  n  r.  r 
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SUBROUT INE  TRANS 

DIMENSION  CFMb(l5)  »TMS(15)  »CF1“)U(15)  »TMD(15) 
DIMENSION  LLW(IO)  *LC(1G) 

DIMENSION  TLOC2(10)  #TLOC3(10)  #TLOC4(.!0)  » 


4>L  Vw  (10)  » I A ( 10 )  » IPE ( 10 )  » IDL ( 10 )  iAV 

COMMON  i*TMUCK»WTlM»NLOCO»TMUCi\»KMAX»NCARS»l.CLAS»Nl>C» ACCMAX, 
j, VELMAX , OECEL , NSC » Dl » Ah  T» HAUL » uSECS , NSW , WTC  AR , FC  AK » TPM , ADR  T  , 

*CROSEC»GAMMA,CTlME»LWHD»TTMf  TNL»TlMAX,NRHG»  ILS» 

bCFR(131 »FTA(i3) .NCLAS1120) »LL(10) »CTLOC(10) »TLOC1(10) t WTLOC (10) » 
bFLOCOI  10)  » G ( 100 )  #0(100)  »DS ( 1 0 )  » WTLOAD  ( 10 )  ►  SPEED  (10)  #LOAU ( 10 )  # 
bl)ISTR(lu)  »TIMc(10)  »TSTOP(10)  #  IPASS  ( 1 0 ) »  WTTPN  ( 10 )  »CFL  ( 13)  ,CFO(13) , 
*C I ( 13) » wTL ( 13) »DSTOP ( 10) »ST ( 13) *CFS ( 13) #TLOAD (10) # TDUMF ( 10) » 
b  CF ( 120 » 13) » TV ( 1^0  » 13) #T ( 10# 3u ) #  S ( lO #  3q ) # LS ( 20 ) » ILWTID 
i»IDL0Ar>»SCCTM»TRlT2»TM0LE2»THlRD2»DELTH»  IDEOS#TNLT 
*»FrAD»TrtOLEl»TBlTl»THiHDl»NCRtW»SHFTMH»HRPSH»NSHlFT»IMAN» 1CYCLE 


SUBROUTINE  TRANS  CONTROLS  THE  SIMULATION  OF  THE  MATERIALS  HANDLING 
SUBSYSTEM 


READ  IN  CONTROL  CARDS 

IF  (  ICYCLE  ,NE.  0)  GO  TO  650 
IF  (  IDuOS  ,EG.  1)  GO  TO  699 
IF  (  11. 5  .NE,  0)  GO  To  652 

READ  IN  INPUT  DATA  FOR  CONTINUOUS  MATERIAL  HANDLING  SYSTEM 


READ  653 »  EFFCSAf  FLUVc.L»GMUCK 

653  FORMAT  (  3F10.3) 

READ  b54»  NCPTS 

654  FORMAT  (15)  .  k  , 

READ  b55  (  CFMS(I ) »TMS(I)  »  1=1 # NCPTS) 

o55  FORMAT  (  6F10.3) 

READ  654 »  NCST 

READ  655  (  CFMD ( I ) » TMu ( I ) »  I=i»NCST> 


initialize  variables  to  zero 

SPWT=G.d 
TUELl' 1=0.0 
TBELT  2=u  •  0 
TDELT3=U,0 

DWNrM=0.0 

SJMDLYrO.O 

CALCULATF  PWT  and  spwt 

Y=RANu(N) 

DO  656  1=1 » NCPTS 
656  IF  (  CFMS(I)  ,LT.  Y 


.AND,  CFMSd  +  1)  .GE.  Y)  GO  TO  657 


nor* 


110 


0 


obB 

ot>7 


o5< 


PRINT  GbB#  y 

Sr0pAT  F12. 6f  3X#  ‘CANNOT  b£  FOUND  IN  CFMb-lMb  CURVL’I 

Z  r=PwT*6o!oT'V,S  ( 1 * 1 * "  fMS  ( 1 * '  *  1  ‘ Y'CFMS  ( 1 }  }  '  1 CFMS  ( 1  +  1 )  -ChMb(I) )  ) 

SPWT=SPrtT+FWr 

IF  (  I Lb  .EQ,  0)  GO  Tu  701 

PRINT  26 

FORMAT  l  *  ILS  WAS  NUT  ZERO*) 

STOP 


CAPMH=FFFCSA*FLUVEL*GMUCK/200y .  0 
TPI4=AURT*UROSlC*GAMMA/  (60 . 0*2u00 . 


MOTM=FTAD*60.0/ADRT 

IF  (  CAPMH-TPM)  659*6o0»660 


0 


) 


COUNT  TIME  DELAYED  AND  MH  WAIlED  DUE  To  HANDLING  SYSTEM 

t>69  AuJTPMsCAPMH 

AMOTM=MOTM*TPM/ADJTPM 
D IF TM=AMOTM-MOTM 
SUMDLYsSUMDLY+DIFTM 
CT I ME  sCTIME+uIFTM 
TDELTls  TBELTl+OlFTKi 

tditi=tditi+diftm 

TM0LE1=TM0LE1+DIFTM 
ruIRDlrTHIRDlfDIFTM 
DIFMH-D1FTM*NCREW/o0.U 
°77  IF  (  SHFTMH-D1FMH)  67a, 676# 6 7o 
o7'j  S(IFTM»i=bHFTMH+NCREW*HKPSH 
NSHIFT=NSHIFT+1 
GO  TO  677 

°?b  SHFTMHsSHFTMH-DIFHH 
TUELT2=rBELT2f0lFTM 

MOTMSAMOTM 

obO  IF  (  SPWT-MOTM)  6bl»6o2,662 
«jd1  MOTM=M0TM-SPWT 

TbELTlrTOELTl+SoWT 
SPulTsO  •  u 


MH  FOR  REPAIR 

Y=RANd(|>|) 

DO  663  j=l»NCbT 

°bJ  PUINTCfo5»JY  *LT*  Y  *AN°*  CFMU(J+D  •  OE •  Y)  GO  TO  664 
oob  FORMAT  (  lX,tRAND=*»  FI2 , 5 # 3X# » NEVER  Be  FOUND  IN  CFMD-IMD’ ) 

ot"  D“NHI '=TMD <  J )  +  ( TMB ( J+l ) -TMD < J ) ) * ( ( Y-CFMd ( J ) ) / ( cFMu ( J+l > -CFM0 <d >  > ) 

c  calculate  downtime 


Ill 


oc.l  IF  (  iHFTMH-DiuNMl I)  66o»6b7»66/ 
ot5(S  DwNMHsCwNMH-SHFTMH 

IK  (  I*' AN  .EQ.  0)  GO  10  856 
DWt  JTM«SHFTMH*bO ,  9/IMAN+OWNTM 
8bo  SHFTMh=NCR£W*HRPSH 
N5MIFTrNSHIFT+l 
l'MAN=i\»CKEto 
GU  TO  6b8 
v. 

C  DETERMINE  TBElTI  AND  WAITING  TIMES 
C 

ob7  SHFTMH=SHFTMH-DWNMH 

D«NTM=pwNTM+DWNMH->oO.O/IMAN 
TDELT  jrTFlELT  3+D',VNTM 
C  T I ME=C  T I  ME+p WNTM 
TuIT2=TUIT2+DwNTM 
rM0LE2r  TMOLE2+DWNTM 
THIRD2Z  miRD2+nWNTM 

OWNTMzn.O 

C 

c  calculate  pwt  and  spwt 

c 

Y=RANU(N) 

DO  6b y  I=1»NCPTS 

oc9  IF  (  CFMR(I)  ,LT.  Y  , AND.  CFMb(I+l)  ,Ge.  Y  )  GO  TO  670 
PRINT  668 »  Y 
STOP 

o7  :  PdTrTKS { I )  +  ( TMS { I  + 1 ) - T MS ( I ) ) * ( ( Y-CFMS ( I ) ) / ( CFMS (1  + 1 ) -CFMS ( I ) ) ) 
PrtTzPwT+60.0 
SPWT=SPrtT+PWT 
6u  TO  6o0 

ob2  SPWT=SPwT-MOTM 

TiJELTisTBELTi+MOTM 

MOTMrO.v 

GO  TO  701 

C 

-  PRINT  SUMMARY  OF  CONTINUOUS  MATERIAL  HANDLING  SYSTEM 
o9?  PRINT  672 

o72  FORMAT  (  1H0>*  SUMMARY  OF  CONTINUOUS  SYSTEM*) 

TuELT lzTBELTl/60.0 
TdELT2=TBELT2/60.0 
TdcLT6=TPELT3/60.0 
SUMOL r rSUMDLY/60 . 0 
PRINT  673 

o73  FORMAT  (  lHOf’  WORK  TIME  TIME  DELAYED  DOWN  TIME’) 

PRINT  674»  TBELTI »TBELT2»TRELT3 
o74  FORMAT  <  3F12.3) 

PRINT  665 »  SUMDLY 

0db$FSYSTLM'1XFl2°3?L  TlMt  DELAYEu  BY  THE  CONTINUOUS  MATERIAL  HANDLING 
GO  TO  7ul 
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C 

c  CYCLIC  MATERIAL  HANDLING  SYSTtM 


C 

c 

c 


If-  (  IDEOS  ,EO.  1)  GO  TO  2*19 

II)LC=0 

ILC=0 


DO  ‘io  I=l»NLOCO 
BO  LC(I)=0 

IE  (  IL.S  .EQ,  0)  GO  To  99 
IF  (  ILwTJD  .EQ,  1)  GO  TO  90 
98  IF  (  LWTIO  «NE»  0)  GO  TO  301 
WTIM=0,u 


99 

IO'1 

101 


GO  TO  301 

forma  } « M  ^°CU ' KMAX  ’  Ni«S .  NSW  -  NC  ARS ,  LCLAS ,  NDC  .NSC 
PRINT  101»NL0C0,NSECS»i\ICARS 

FORMAT  <1H1»3X#» NUMBER  OF  LOCuMOTivES  -*,I«5  /.ft*  3*. 

-numblr^of  grade  sections  =*  *  is  ;fXcars  per  train 


read  in  all  characteristic  curves  of  locomotives 


ir-  f(*J(I'J,'S(I»'J)»  J=1»KMAX)#  I=l#NLOCO) 

ii<_  format  (  gfio.3) 

c 

c  RhAD  in  THE  cumulative,  frequency  curves  of  loading  and  dumping 


c 

c 

c 

c 

c 

c 

c 

e 

c 

c 

c 

c 

c 

c 

c 

c 


RLAD  ICO  (  CFL(I) »wTL(X) » 
READ  103  (  CFS(I)»  ST ( I ) # 
READ  103  (  CFD(I) »  C T 1 1 ) # 
103  FORMA T(  6F1Q.3) 


1=1 #LCLAS) 
1=1 #NSC ) 
1=1 #NDC ) 


READ  in  THE  WEIGHT  ANU  FRICTIuN  COEFFICIENT  of  the  LOCOMOTIVES 
read  103  (WTLOC(I)fFLOCO(I),  i=l#NLOCO) 

READ  IN  THE  PROFILE  OF  THE  TUNNEL 
READ  103  (  G(1)#D(I),  1=1 f NSEcS) 

READ  IN  THE  WEIGHT  ANU  FRICTION  COEFFICIENT  OF  THE  CARS 
READ  103  »  WJCARfFCAR 


READ  IN  SPEED  LIMITATIONS 


READ  103#  ACCMAX»VELMAX»DECEL 


READ  IN  THE  DISTANCE  oETWEEN  SWITCHES#  DELTH#  AND  THE  NUMHPP  ne 
SWITCHES  BETWEEN  THE  DUMPING  STATION  AND  THE  PORTAL 


r.  r  r 


]  ]  3 


KLmD  104»  DISw »DF.L FH 
104  FORMAT  {  2F10.3,  15) 

INITIALIZE  VAKl ARLES 

DO  105  l=l#NLOC0 

niSTRmzo.o 

LOAD ( I ) =0 
WTLOAU (1)50*0 
WT  TRN ( I ) =0*0 
TLOC1 { I )=C.O 
TLOC2(I)=0.0 
TLOCJ(I)=G.O 
TLOC4(I)=0.0 
TlMEtD-O.O 
CTLOC(Ii=0.0 
SPEED (1)50.0 
TFAS5(I)=0.0 
TSTOP ( I ) =0 . 0 
□STOP ( I ) zO .0 
5  CONTINUE 
LIL=0 
(4lDE=0 
NLUL=0 
HaUL=0.u 
TSEC=0 , 0 
NS  A  50 
SAFTsO.U 
WTDzO.O 
WT  GzG  •  n 

LOCATE  and  COUNT  THE  SWITCH  POINTS 

DO  10b  J=1*NSL'P 

106  HAUL=HAUL+D(J) 

TNLT= INl+HAUL 

NSW=2 

DHSHAUL 

10<:>  IF  {  OH-DISW)  107» 1 08 » 1QS 
10M  NSW=NSW+1 
DH5DH-DISW 
GO  TO  109 

107  Nw=NSW-l 
DMSsUH 

IF  (  NN  ,EQ.  1)  GO  TO  110 
DO  111  I=2»NN 
111  DS{I)=DISW*(I-1) 

110  DS(1)=0,0 
OS(NSW)=HAUL 


r.  r.  cs  o  r.  o 


114 


COUNT  THE  NUMBER  OF  SUCTIONS  *T  THE  STftRT  OF  THE  SIMULATION  ANU 
DETEkMINE  THL  DISTANCE  TO  THt  LAST  SECTION 

Nb£CS=NSOP+l 

DU=0(NStCS) 

O(HSECS)=0.0 

LOCATE  fHt  LOCOMOTIVES  AT  THE  STARTING  POINT 

DO  199  l=l»NLOCO 
199  IDE ( I ) =U 
NLDESU 
NLHLsO 

IF  (  Nl.OCO-NSw)  311  * 112» 1 13 

311  NImLDsU 
N2SNSW-WL0C0  +1 

DO  114  1=  NSW»N2»-1 

jiNSW-I+l 

LS(I)=2 

114  LL(I)=J 
N1=N2-1 

DO  115  I  =  1»N1 

115  LS(I)=0 
GO  TO  116 

112  NWLD=1 
IUE(NLOCO)=l 
NLDE=i 

DO  117  1=1 > NSW 
JzNSW-I+l 
LS(I)=2 
117  LL(I)=J 
GO  TO  116 

113  NwLD=iMLOCO-NSW  +  l 

lmn=o 

DO  312  i=NSW,NLOCO 
LMN=LMN+1 

312  IUE(I)=LMN 
NLDE=NWcD 

Do  lib  I=1»NSW 
jsNSW-1+1 
LS  ( I )  =?. 

lie  ll(I)=j 
NEXsNSW+1 
DO  97  I=NEX,|MLOCO 
97  Lu(I)=l 

116  CONTINUE 

DO  12o  I=l»NLOCO 
12b  IUL(I)=0 

IF  (  ILS  .EQ,  0)  GO  TO  701 


r  r-  r.  o  n  n  r. 
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START  THE  SIMULATION  rtlTH  LOCoMOTlvE  NUMDtTR  1  AT  THE  LOADING  POINT 


90  NL=1 

CALL  LOADNG(NL) 

IuLOAU=i 

CCTMrwTLOAD (NL ) *2000 , 0*60 , 0/ ( CROSEC*GAMMA*ADRT ) 

SCCTM=SCCTM+CCTM 

TLOC 1 ( NL ) =TLOC 1 ( NL ) *  TuOAD ( NL ) 

C  TLOC  (  Nl  )  =CTL0C  ( NL )  +Ti-OAD  ( NL ) 

SAFT=SAFT+AFT 

WTGrWTG+WTLOAUtNL) 

N  j -LL (NL) 

LS ( NS )ri 


DETERMINE  THE  LOADING  TIME  FOk  THE  LOCOMOTIVE 

DO  119  J=2»NL0C0 
TL0C4 ( J)=TLOAD(NL) 

119  CTLOC(J)=TLOAU(NL) 
c 

CALL  MOT  ION (Nu# NS) 

C 

CT  LOC ( Nu ) =CTLOC ( Ml ) +T 1 ME ( NL ) A>0 . 0 

TLOC 3 (NL) “TLOC 3 ( NL ) +T1ME (NL ) /oO , 0 

IF  (LOAD (NL)  .EQ.  1)  GO  TO  12o 

LS(NS)=lS(NS)-2 

NS=I'JS+1 

LL(NL)sNS 

LS(NS)=LS(NS)+2 

NSlsNS-i 

NS2=Ni> 

GO  TO  121 

120  LS(NS)=LS(NS)-1 
NS=NS-1 

LL (NL) =NS 
LS(NS)=LS(NS)+1 
NSlzNS+1 
NS2=NS 

121  CONTINUE 

IF  (  LOoPRT  .NS.  0  )  GO  TO  301 
PRINT  411 

411  FORMAT  (  1H0»  ’LOCO  NO.  CLOCK  TIME  MOVING  Sv¥  NO.  LOAD'/’  *r25X 

’FROM  TO*  /) 

PRINT  4i2»  NL»CTLOC(Nt-)  »NS1»NS2»L0AD(NL) 

412  FORMAT  (  4X,i3»4X,Fl0.2»5X»I2»4X»I2»7XrIl  //) 


o  r~.  n  n  r- 
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CMOOSL  THE  LOCOMOTIVE  HAVING  THE  SMALLEST  VALUE  OF  CTLOC(I) 

•301  IL-0 
M-l 

00  122  J=2»NL0C0 

122  If  (CTi_OC  (M)  ,oT,  CTLOC(J))  Mru 
NL=M 

empty  tkain 

Ip  (  LIL  .EQ,  0)  GO  Tu  200 
DO  221  I=i»LlL 

221  IF  (  NL  .EQ,  LLW(I))  GO  TO  222 
GO  TO  200 

222  Ms I 

IF  (  LLW(M)  ,£Q.  LIL)  GO  TO  224 
LLW (M) zu 
M1=M+1 

00  223  JsMl'LIL 
KsJ-1 

223  LLW(K)=lLW(J) 

22'*  LILsLlL-i 

200  It-  (  LOAD(NL)  ,EQ.  1)  GO  TO  12-3 
NSsLL(NL) 

IF  (  NS  *EG.  NSW)  GO  fO  124 
IF  (  NS  .EG.  1)  GO  TO  125 

GJ  T0L361S+1)  *E°*  2  *°R*  LS(|',S+1)  ,EQ»  3)  60  T0  3&0 
30'  IF  (  1LC  .EQ,  0)  GO  TO  129 
NSS=NS+1 
DO  362  Isl.NLOCO 

362  IF  (  LL(I)  .EQ.  NSS  )  GO  TO  3tj3 
PRINT  364.  NSS 

304  FORMAT  (  »  NSS  AT  364  IN  TRANS’.  15) 

STOP 

363  II=I 

IF  (  wOAD(II)  ,EQ,  0)  GO  TO  3ot> 

Jl=II+l 

DO  363  KrJl.NLOCO 

36b  If  (  LL (K )  .EQ.  NSS)  GO  TO  367 
PRINT  3o4#  NSS 
STOP 
3o7  II=K 

3o6  DO  366  irl.ILC 
368  IF  (  LC(I)  .EQ.  II)  GO  TO  369 
GO  TO  129 
3o9  ILC=ILC+1 
LC(ILC)=NL 
IDLC=1 
~*»G~  TO  129 
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3o  1  CTLOC  ( ML )  =CTLOC  ( NL )  +TPASS (NL )  /GO . 0 
D 1 5TU ( ML ) =0 ISTR ( NL ) +Db 1  OP ( ML ) 

TL0C3 (ML ) =TL0C3 ( ML ) +TPASS ( NL )  /6U . 0 

CALL  MOTION  (NL>NS) 

IF  (  uOAO(NL)  .EQ.  1)  GO  TO  1^7 

LS(NS)=LS(NS)-2 

NS=NS+1 

LL(NL)=.\|S 

LS (NS)=LS (NS ) +2 

NSlrNS-i 

NS2=MS 

GO  TO  120 

127  LS(NS)=LS(NS)-1 
NS=NS-1 
LL(NL)=MS 
LS(N$)=LS(NS)+1 
NS1=NS+1 
NS2=Nb 
12F  CONTI NHL 

CTLOC ( NL ) =CTLOC ( NL ) +T ilC  ( NL ) /oO . 0 
TLOC3(ML)=TLOC3(NL)+T*fC(NL)/o0.u 
IF  (  LOGPNT  , ME ,  0)  GO  TO  413 
PRINT  411 

PRINT  412 »  NL»CTLOC  (Nl)  »NS1  »Ml»2»L0AD(NL) 
413  IF  (IL  .EG.  0)  GO  TO  OU0 
IF  (L1L  .EG.  u  )  GO  TO  225 
150  192  1  =  1  »IL 
DO  220  J=i»LlL 

220  IF  (  LW ( I )  .tQ.  LLU (0) )  GO  TO  229 
ATPS=0.0 
GO  TO  230 
229  ATPS=f-TPS 
200  JJ=LW ( I ) 

C  T  LOC ( j  J ) =CTLOC ( J J ) + ATPS+T I ML ( NL ) /60 , 0 
192  TL0C4 ( Jo ) =TL0C4 ( JJ ) +AT  PS+TIME ( NL ) /60 , 0 
DO  231  I=1»IL 
DO  232  0=1 >LIL 

2.32  IF  (  LW ( I )  .EG.  LLW(J))  GO  TO  231 
L1L=LIL+1 
LLW(LiL)=LW(l) 

2.31  CONTINUE 

FTPS=TST0P(Nl)/60.o 
DO  35  1=1 »IL 
35  LU ( I ) =0 
IL=0 

GO  TO  3oC 
225  DO  22b  1=1 »IL 
JJ=LW ( I ) 

C  TLOC  ( JJ )  =CTLOC  ( JJ )  +T1ME  ( ML )  /*>0 . 0 
22b  TL0C4(JJ)=TL0C4( JJ) +T1ME (NL) /60.0 


n  r. 
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LIL=IL 

FlPS=15IOP<NL)/f>0.  0 
DO  227  1=1#LIL 
227  LLW ( 1 ) =LW ( I ) 

Ou  3b  1=1, IL 
3G  LW(I)=0 
IL=0 
Lto ( I ) =0 
GO  TO  300 

IF  (  LS(NS)  ,LQ.  3  .Ok.  LS(N5)  .Eg,  1)  GO  TO  338 
GO  TO  339 

33*5  If-  (  ilc  ,EQ,  0)  00  TO  129 
DO  331  I=1»NL0C0 

3j.I  IF  (  LL(I)  .EG.  fJS  .AwD.  LOAD(I)  ,E0.  1)  GO  TO  332 
PI  <  I  NT  541  ° 

b41  FORMAT  (  IX, ’STOPPED  aT  541  In  TRANS* ) 

S  TOP 

332  INLC=I 

00  333  1=1 » ILC 

333  IF  (  LC(I)  .EG.  INLC)  GO  TO  334 
GO  TO  129 

334  ILC=IlC+1 
LC(ILC)=NL 
IUl.C=l 

GO  TO  129 

339  CTLOC l NL ) =CTLOC ( ML ) +TSTOP ( ML ) /60 . 0 
D I S TR  ( ML )  =1)  I STR  ( ML )  +DbTOP  ( NL ) 

TLOc3(NL)=TLOC3(NL)+TbTOP(ML)/60.0+TIME(ML)/60.0 

SPEED (ML) =0.0 
TPASS(ML)=0.0 
TSTOP(ML)=0.0 
DSTOP (Ml) =0,0 
TIME(UL)=0.0 

CLAST  s2lTCHNNEL  LENe,H  AND  HaULagE  DISTANCE  AND  RELOCATE  THE 

IF  (  SAFT  .GE,  DELTH  )  GO  TO  id 0 
GO  TO  189 

10 U  SAFT=0AFT-DEL1H 

haul=haul+uelih 

TSEC=rSEC+DELTH 

dms=dms+pelth 

IF  (  TSEC-DD)  181 , 182,182 

102  tsec=tsec-dd 

D (MSECS )=DD 
NbECS=N5ECS+l 
Du=0 (MSECS) 

1B1  D(NSECS)=TSEC 

IF  (  UMS-DISW  )  183,184,184 
194  DM3=DhS-niSW 
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OS  ( NS W  )  =HAU|_-DMS 
NSW=NSW+1 
1U3  DS(NSiV)=HAUL 

IF  (  HAUI.-TNLI)  189» lub»  185 
1Gb  PKlK'l  J86 

XUf>  GOKfoT,99H°’  COKPLETEO  AND  SIMULATION  TERMINATED'. 

109  COIJTlf'iUt. 

IF  (  LWTID  ,|\|L.  0)  GO  TO  700 
WTIM=CTLOC  (nd-ctime 
II  (  IL  .EG,  0)  GO  TO  701 
00  70*2  1  =  1  »IL 

J  J=LW ( l ) 

Lrt(I)=G 

CTC=C  (IOC  (ML) -CTLOC  (JJ) 

CTLOC( JJ)=CTLOC(NL) 

/Ov?  T|.0C4(jj)=TL0C4(  JJ)+CIC 

1L=0 

GO  TO  7ol 

70 IF  (  IOLOAO  .tO.  1)  GO  TO  600 
IOLOAL)=i 


ChLL  LOAONG(NL) 

0 


191 


oOl 


9Jl 


err  <NL  1  *20O°  * u*60  •  0/{CHOSEC*GAMMA*ADRT  .* 

SCCTM=SCCTM+CCTM 

vJTu=WTG+WTLOAlj(NL) 

CTLOC ( Hu ) =CTLOC ( ML ) +TLOAD ( ML ) 

TLOC 1 ( NL ) =TLOC 1  ( NL ) +TuOAD ( ML ) 

LS(NS)=lS(NS)-1 

saft=saft+aft 

IF  (  1L  .EQ.  0)  GO  TO  300 

00  191  1=1 # IL 

J0=LW(n 

LU ( I ) =o 

C  TLOC ( JJ) =CTLoC ( JJ)  +fL0AD(ML) 

TlOC4( JJ)=TLGC4(JJ)  +  TLOAD(Nu) 

IL=0 

GO  TO  3u0 


WVJTM=  (CTIME+S0CTM)-CTi-0C  ( NL ) 
IF  (  kvWTM  )  6C1 f 601 »6u2 
IF  (  IL  .EG.  u)  GO  TO  701 
DO  931  1  =  1 #  II 
JJ=LW ( J ) 

Lw(I)=C 


TLOC4  ( JJ)  =TL0C4  ( Jj)+{ oTLOC  { NL) “CTLOC  ( Jj) ) 
CTLOC ( JJ) =CTLOC ( JJ) + (CTLOC (NL) "CTLOC { Jj) ) 


GO  TO  701 

t>C2  CTLOC  ( NL)  =CTLOC  (ML)  +WwTM 
TL0C4(NL)=TL0C4  (NL)+Wv»TM 
IF  (  IL  .EG,  (j)  GO  TO  300 


•  n  a  r  r» 
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DO  39  I  —  1  r  I L 
JJ=L'.v  { T  ) 

Lm ( I ) =0 

CTLOC ( JJ)=CTLoC< JJ)+WwTM 
TL0C4 ( ml ) =TL0C4 ( NL ) +WwTM 

I  L-0 

GO  TO  3u0 

*29  CTLOC  ( ML )  =CTLOC  ( PJL )  +T3T0P  ( NL )  /  60 . 0 
DISTK(Ml)=OISTR(NL)+DSTOP(NL) 
rLOC3{NiL)=TLOL3(ML)+TGTOP(NL)/6C.O 
SPEED (Nu> =0,0 
TPASSlflL)=0,0 
TSTOP(ML)=0.0 
OSTOP(KL)=0.0 
T1ME(NL)=C.0 

If-  (  iOLC  .NE.  1)  00  10  150 
IDLC=0 

SSCC=CTTM-CTLOC(NL) 

IF  (  i,SCC)  335>335»33o 
o36  CTLOC  (fJL)=CTLOC(  ML  )+SbCC 
TL0C4(fiL)=TL0C4(flL)+SbCC 
035  IF  (  II.  ,EO.  0)  00  TO  340 
DO  337  1=1 »IL 
JJ=Lh ( I ) 

Lto ( I )=0 

TL0C4 ( J  J ) =TL0C4 <  J  J )  + (CTLOC ( NL ) “CTLOC ( J  J  >  > 

CTLOC (JJ)=CTL0C (JJ) + (CTLOC (NL) “CTLOC (Jj) ) 

ILC=ILC+1 
037  LC(ILC)=JvJ 
IL=0 

GO  TO  340 

15;  IL=IL+1 
LW(  IL)=ImL 

IF (IL  ,EQ.  1)  GO  TO  130 
ARRANGE  LfcdLi  IN  ASCtNDlNG  ORDER 

DO  131  1=1 » IL 

131  XA(X)=Lw(X> 

MAsiL-l 

DO  132  0=1 »NA 
M=J 

MA”J+i 

DO  133  1=MA»IL 

133  IF  (  IA ( I >  .LT.  IA (M) >  M=I 
I  IL’1P=IA(0) 

I A ( J)=XA(M) 

132  IA(M)=ITEMP 

SEARCH  FOR  THL  SHORTEST  VALUE  OF  CTLOC (I)  EXCEPT  FOR  CTLOC (1L) 

ML=1 

Nl=2 
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DO  134  1  =  1  * IL 
NJ=IA(  I) 

IF  (  NJ  ,LQ.  ML)  GO  To  56C 
Ni\=I 

60  TO  5ol 
bt i'j  Ml=ML  +  1 


ni=ml+i 


IF  (  I  .ME.  IL)  GO  TO  134 
IF  (  Ml  .EG,  Nl.OCO)  GO  TO  139 
GO  TO  501 


104  CONTINUE 
boi  Du  562  j=MK,il 
NJ=IA(J) 

IF  (  N1  .EG.  NJ)  GO 
135  N2=NJ-1 

UO  137  K=N1,N2 
107  IF  (  CTLOC(ML)  .GT. 
bb3  N1=NJ+1 
bb2  CONTINUE 

IF  (  NJ  ,EQ.  NLOCO) 
bbl  DO  140  jsNl» NLOCO 
1*»G  IF  (  CTLOC(ML)  ,6T. 
139  NL=ML 


TU  563 


CTLOC {K ) )  ML=K 


GO  TO  139 
CILOC(J))  ML=J 


GO  TO  J?uO 
130  NJ=LW(IL) 

IF  (  NJ  .EG.  1)  GO  TO  141 
ML=1 


Ml=2 

IF  (  NJ  .EG,  2)  GO  TO  142 
N2SNJ-1 

DO  143  J=N1 »N2 

l‘+3  IF  (  CTcOC(ML)  .GT.  CiLOC(J))  ML=J 
142  N3=IJJ+1 
N4=NL0C0 

IF  (  NJ  ,E0.  NLOCO)  GO  TO  144 
1)0  14b  J=N3»N4 

145  IF  (  CTLOC (ML)  .GT.  CTLOC(J))  ML=J 
144  NI.£ML 

GO  TO  2o0 
141  ML=2 
Ml=3 

N2=NL0C0 
DO  246  J=Nl t ML 

<^46  IF  {  CTLOC  ( ML )  ,GT.  C1"L0C(J))  ML=J 

NlsML 

GO  TO 

*25  IF  (  NLOE  .EQ.  1)  GO  10  146 

IF  (  IDL(NL)  .FQ.  1)  oO  TO  14o 
IF  (  11.C  .EQ,  0)  GO  TO  150 
N11=I0E(NL)-1 
00  341  I  =  1 » N 1 i 
DO  342  J=l» NLOCO 


122 


342  IF  (  infc(J)  .to.  1)  eu  TO  343 
PRINT  344 ,  I 

344  FORMAT  (  *  ID£(NL)=I  AT  344  I„  TRANS  COULD  NOT  BE.  FOUNU*, 

STOP 

343  ONLzJ 

00  34b  K=1 » ILL 

34b  IF  (  JML  .EQ.  LC(K))  bO  TO  34b 
GO  TO  341 


346 


372 


ILC=lLC+l 
LC ( ILC ) -ML 
SbCC=CTTM-CTLoC (NL) 

C  TLOC  ( NL)  =CTLOC  ( N'L )  *SbCC 
TL0C4 (NL)=TL0C4 (NL) +SbCC 
IF  (  11.  .£0.  0)  GO  TO  371 
00  3/2  t\=l>lL 
JJ=LW(«) 

Lh'  (K)=n 

ilc=ilcfi 

LC(ILC)=JJ 

TL0C4  ( JJ)=TL0C4  ( JJ)  +  (C.TTM-CTL0C  ( JJ)  ) 

C I LOC  t  J J ) =CTLOC ( J J )  +  ( c TTM-CTLOC ( J j ) j 
Il=C 


GO  TO  371 

-371  IF  *  I  *£«•  Nil)  GO  TO  425 
NIOsI+l 
GO  TO  347 

34i  CONTI NHL 
GO  TO  lb(? 

347  DO  346  L=N10,Nll 
DO  34 y  0=1»NL0C0 

J49  <  int(J)  .to.  L)  GO  TO  350 

PRINT  351 »  L 

351  FORMA |  (  IX, STOPPED  AT  351  Km  TRANS*,  15) 

3b ?  JJ=J 

TLOC 4 ( J J ) = TL0C4 ( Jj ) 4 (  CTTM-CTLOC ( Jj ) ) 

C  TLOC ( JJ)-CTTM 

ILC=IlC+1 

LC(ILC)=JJ 

34  8  CONTINUE. 

*+25  IF  (  1PL(NL)  ,E0.  NLDL)  GO  To  340 
N13=Ijf;(NL)+i 
DO  420  I=N13,nLDF 
DO  421  J=1,NL0C0 

*+21  IF  (  int(J)  ,L0.  I)  GO  TO  422 
PRINT  423 

423  FORMAT  (  IX, *  STOPPED  AT  423  IN  TRaNS*) 
STOP  M  3  1 

422  JJ=J 


XLCrlLCM 


15) 


c 

c 


LC  ( I!-C)=JJ 
SSCC=CTTM-CTLOC< JJ) 

C 1 LOC { J J ) =CTLOC ( J J ) +SiCC 
TLOC4{ JJ)=TL0l4 (JJ)+SSCC 

420  CONTINUE 
GO  TO  340 

a46  IF  {  LS(NS+1)  .FQ.  0  .ON,  LS(NS+1) 
IF  (  ILL  .EQ.  0)  GO  Tu  150 
NNSSsi-JS+1 


,  EQ , 


U  GO  To  147 


00  352  1=1 »  NLOCO 

IF  {  LL  ( I )  .Em.  fJNSS)  GO  TO  3bi 

PRINT  354 t  NNbS 

J54  FORMAT  (  IX,  t  STOPPED  AT  345  lii  TRaNS*.  IM 
STOP  ' 

553  II=I 


IF  (  LOAD (II)  ,E0,  0)  GO  TO  3b6 
J1=II+1 

DO  35b  K=Jlt NLOCO 

555  IF  (  lL(K)  .eg,  NNSS)  GO  TO  3^7 
PRINT  354,  NNSS 
STOP 
557  I I=K 


556  00  358  I=lrlLC 
550  IK  (  LC(I)  .EG.  II)  GO  TO  359 
GO  TO  lo9 

559  SSCCSCTTM-CTLOC(NL) 

TL0C4 ( NL ) =TL0C4 ( NL ) +S5CC 

CTL0C(NL)=CTL0C(NL)+S5CC 

ILC=ILC+1 

LC(ILC)=NL 

IF  (  IL  ,EQ,  0)  GO  TO  340 

DO  370  J=1»IL 

JJ=LW(J) 

Lw ( J) =0 

TL0C4 i J j ) =TL0C4 ( JJ) + (C TTM-C TLoC ( J J ) ) 
C  TLOC ( J J ) =CTLOC ( J  J )  +  <  C  TTM-C  TLoC ( J j  j ) 
ILCsILC+ft 
370  LC ( ILC ) =JJ 
IL=0 

GO  TO  340 
147  NLl)E=NLDE-l 

DO  148  I=l»NLOCO 
IF  {  IDt. ( I )  ,EQ.  0)  GO  TO  148 
IDE ( I ) =1PE ( I ) -1 
140  CONTINUE 


CALL  MOTION  (  N|_,NS> 

C ILOC { NL ) -C TLOC ( NL ) +T1ME ( NL ) /oO , 0 

ns=ns+i 

Ll(NL)=NS 

LS(NS‘=LS(NS)+2 


r  r>  n 
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IF  (  LOGPRT  .InjE .  0)  GO  TO  165 

NS1=NS-1 

NS2=Nb 

PRINT  411 

PRINT  4l2»  NL»CTLOC(No) »NSl,Nb2,L0AD(NL) 

165  IF  (  1L  .EQ.  0)  00  TO  300 

DO  151  1=1 t IL 
JJ=LW(I) 

LW(I)=0 

C TLOC  ( J 0 )  =C TLOC  { J  J )  +T 1ME  ( NL ) /oU .  0 
151  TL0C4 ( J J ) =TL0C4 { J J ) +T I  ME (NL)/uO,u 
IL=0 

GO  TO  300 

TRAIN  LOADED 

123  NS=LL(NL) 

IF  (  NS  .EQ.  1)  GO  TO  155 
IF  (  NS  .EQ.  NSW)  GO  10  156 
IF  (  NS  .EQ.  2)  GO  TO  157 

IF  (  LS(MS-l)  .EQ.  1  .OR.  Lb(NS-i)  .EQ.  3)  GO  TO  460 
GO  TO  157 

460  IF  (  ILC  .EQ.  0)  GO  TO  129 
NSSsNS-1 

Do  462  1=1»NL0C0 

462  IF  (  LL(I)  .EU.  NSS)  uO  TO  460 
PRINT  464 1  NSS 

464  FORMAT  (1X» ’STOPPED  AT  464  IN  TRANS' t  15) 

STOP 

403  11=1 

IF  (  LOAD(II)  .EQ.  1)  GO  TO  4o6 
Jlsll+l 

DO  465  K=Jl»NLOCO 

465  IF  (  LL(K)  .EQ.  NSS)  GO  TO  46/ 

PRINT  464 1  NSS 

STOP 
467  II=K 

4o6  DO  466  1=1 » ILC 

466  IF  {  LC ( I )  .EQ.  II)  GO  TO  469 
GO  TO  129 

469  ILC=ILC+1 
LC(IU)=NL 
IULC=1 
GO  TO  129 

157  CTLOC(ML)=CTLOC(NL)  +IPASS(NL)/6U.O 
0I5TR  ( NL )  =D I STR  ( NL )  +05T0P  ( NL ) 

TL0C3 ( No ) =TLOC3 { NL ) +TPASS ( NL ) / 60 . 0 
LS ( NS ) =LS ( NS ) -1 


C 


CALL  MOTION(NL’NS) 
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CTLOt  (NU)  =CTLOC  <N|_ )  +T1ME (HL) /o« .  0 

TLUC3  (ML )  =TL0C3  ( N|_ )  +T  IMF  <  NL  )/oO,0 

NS=NS-1 

LL(NL)=nS 

LS (NS) =LS (NS) 4 1 

IF  (  LCOPRT  ,wE.  0)  Gu  TO  414 

NSl=NS+i 

NS2=N5 

PRINT  4x1 

PRINT  412 1  NL»CTLOC(Nl.)  »NSl»Nbii»LOAD(NL) 
414  IF  (  MS  .EG.  1)  GO  TO  160 
IF  (  IL  .EG.  C)  GO  TO  3C0 
IF  (LIL  ,EQ,  0  )  GO  Tu  525 
DU  162  1=1 >IL 
DO  52b  J=1»LIL 

526  IF  (  LW ( I )  .Eg,  LLw(J))  GO  TO  52 9 
ATPS=0.U 
GO  TO  550 
529  ATPS=FTPS 
53.2  JJ=LMI) 

C  TLOC ( JJ ) =CTLOC ( J J ) +A \ PS+T I ME ( »C ) /60 . 0 
162  TLOC4  ( J  J )  =TL0C4  ( J  J )  +/\  1 PS+TIME  ( ML )  /60 . 0 
DO  531  1=1 »IL 
DO  532  J=1»LIL 

532  IF  (  LW(I)  .EG.  LLW(JJ)  GO  TO  531 
LIL=LIL+1 
LLW(LxL)=LW(I) 

531  CONTINUE 

FTPS=TSTOP(NL)/60.u 
DO  37  1=1 »IL 
37  LW ( I ) =G 
IL=0 

GO  TO  3oO 

525  DO  526  1  =  1 » I|_ 

JJ=LVJ  ( I ) 

CTLOC ( JJ ) =CTL0C ( JJ ) +T1ME (NL ) /bu . 0 

526  TLOC4 ( J J ) =TL0C4 ( J J ) +T 1ME ( NL ) /od . 0 
L1L=IL 

FTPS=TSr0P(NL)/60.o 
DO  5 27  1  =  1 » LIL 

527  LL.V(I)=LW(I) 

DO  36  1=1 r IL 

36  Lrt(I)=0 
IL=0 

GO  TO  3u0 

160  CTLOC ( ML ) =C TLuC ( ML ) +T5T0P ( ML ) /60 . C 
DlSTR ( ML ) =D I STP ( NL ) +05  TOP ( ML ) 

TL0C3 ( Ml ) =TL0C3 ( NL ) +T5T0P ( ML ) /60 . 0 
SPEED (ML) =0,0 

nldl=mlul+i 

IDL(ML)=NLUL 
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IF  (  IL  .E.Q,  o)  GO  TO  159 

DO  165  1=1# IL 

JJ=LW(I) 

LW ( I ) =p 

CTl.OC  ( JJ)  =CTLoC  ( JJ) +TSTOP ( NL ) /60 •  Q+TIMe  (NL)  /60  »  0 

163  TLOC4(JJ)=TLOC4(JJ)+T^TOP(NL)/60.0+TIME(NL)/60.0 
Il“0 

159  IF  (  NLDL  .E<3.  1)  GO  TO  161 
GO  TO  3oO 
C 

161  CALL  OUMP(NL) 

C 

C  TLOC ( NL ) =CTLOC ( ML ) +TUUMP ( NL ) 

T  L0C2 { ML ) =TL0C2 ( Ml ) +TLUMP ( ML ) 

MLDL=l\ILuL.-l 

I0L(IML)=IDL(NL)-1 

nlde=nldf+i 

IL)£(NL)=NLDE 
LS(NS)=LS(NS)-1 
WTTRN(NL)=WTTKN(NL)“W I  LOAD (ML) 

WTD=W  TD+WTLOAD ( ML ) 

WTLOAU ( NL ) =0 , 0 
GO  TO  3o0 

GO  loiir1’  *EQ*  1  *°R*  *S)  60  TO  500 

50 C  IF  (  iLC  .EQ,  0)  GO  TO  150 
NN5S=NS-1 
DO  5C1  1  =  1 #  NLOCO 

501  IF  (  LL  ( I )  .EO.  MNSS)  GO  TO  5U2 
PRINT  503 

503  FORMAT  (  IX,  ^TOPPED  503  lix  TRANS*) 

STOP 

J02  II=I 

IF  (  LOAD (II)  .EO,  1)  GO  TO  5W4 
J1=II+1 

DO  505  K=JlrNLOCO 

505  IF  (  LL(K)  .Eg.  MNSS)  GO  TO  5u6 
PRINT  507 

507  FORMAT  (  1X,» STOPPED  «T  507  In  TRANS*) 

STOP 

5C6  I I=K 

504  00  506  1=1# ILC 

508  IF  (  LC ( I )  .EG.  II)  GO  TO  509 
GO  TO  150 

509  SSCC=CTTM-CTLOC(ML) 

TL0C4 ( NL ) =TL0C4 ( NL ) +S5CC 
CTLOC (Ml)=CTLOC (ML) +S5CC 
ILC=ILC+1 

LC  ( ILC )  =NL 

IF  (  XL  .EQ.  0)  GO  TO  340 
DO  510  J=1 # IL 
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JJ=LW{ J) 

Lw ( J)=0 

TL0C4  ( JJ )  =TL0C4  ( JJ)  +  ( cTTM-CTLuC  ( J  J )  ) 

CTLOC  { JJ )  =CT|_0C  ( J  J )  ♦  { CTTM-CTLuC  { J  J )  ) 
ILC=1LC+1 
bir  LC(ILC)=JJ 
IL=0 

GO  TO  340 

bll  LS<NS)=LS<NS)-1 

CALL  MOTION  (  NLHJS) 

NS-NS-1 

L5(NS)=lS(NS)+1 
LL (NL)=NS 

Cl LOC ( ML ) =CTLOC ( NL ) +T 1 ME ( NL ) /o0 . 0 

TL0C3 ( NL ) =TLOC3 { NL ) +T iME <  NL ) /toO . 0 

IF  (  LOGPRT  .HE.  0)  GO  TO  165 

NS1=NS+1 

NS2=Nb 

PRINT  411 

PRINT  412»  NL»CTLOC(NU »NSl»Nb2»L0AD(NL) 

GO  TO  le>5 
165  NLOL=NLUL+l 
IJL(Nl)-NLOL 

IF  (  NLOL  .EG.  1)  GO  10  167 
IF  (IUL(NL)  ,E9.  1)  GO  TO  160 
IF  (  ILL  .EQ.  0)  GO  TO  150 
Nll=IbL(NL)-l 
DO  301  I=1»N11 
DO  302  J=1»NL0C0 

302  IF  (  10L(J)  .tQ.  I)  GO  TO  383 
PRINT  304 1  I 

J04  FORMAT  (  lX»tbT0PPE0  AT  304  IN  T«AN5*»  15) 
STOP 

363  JNL=J 

DO  385  K=1»ILC 

305  IF  (  JNL  »EQ.  LC(K))  GO  TO  30o 
GO  TO  381 
386  ILC=ILC+1 
LC ( ILC ) -NL 
SSCC-CTTM-CTLOC (NL) 

CTLOC (NL)=CTLOC (NL)+SbCC 
TL0C4 ( ML ) =TL0C4 ( ML ) +SbCC 
IF  (  1L  .EQ.  0)  GO  TO  391 
DO  392  K=1»IL 
JJ=LW(K) 

Lfl(K)=0 

ILC=ILC+1 

LC(ILC)=JJ 

TL0C4 ( JJ) =TL0C4 ( JJ) + (CTTM-CTLOC ( JJ) ) 

39?  CTLOC { JJ) =CTLOC ( JJ) + (oTTM-CTLoC { JJ) ) 
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IL=0 

-391  IF  (  I  ,EQ.  Nil)  GO  TO  42 b 
N13=I+1 
GO  TO  307 
<301  CONTINUE 
GO  TO  100 

087  00  380  LrMlOf Nit 
DO  389  JrlrNLOCO 

•309  IF  (  IDL(J)  ,EQ,  L)  GO  TO  400 
PRINT  4ul 

401  FORMAT  {  IX,  'STOPPED  AT  401  XN  TRANS*) 
STOP 
400  JJZJ 

TL0C4 ( JJ)=TL0C4 ( JJ) + (CTTM-CTLoC (OJ) ) 
CTLOC( JJ)=CTL0C(JJ)+(CTTM-CTL0C( JJ) ) 
ILC=ILC+1 
LC(ILC)=JJ 
308  CONTINUE 

426  IF  (  IDL(NL)  ,EQ.  NLDL)  GO  TO  340 
N13=I0L (ML ) +1 
DO  430  I=N13,nLDL 
DO  431  JrlrNLOCO 

431  IF  (  IDL(J)  ,EQ.  I)  GO  TO  432 
PRINT  433 

433  FORMAT  <  IX, 'STOPPED  AT  433  Iu  TRANS') 
STOP 

432  jjr J 
ILCrlLC+1 
LC ( ILC ) r JJ 
SSCC=CTTM-CTLOC ( JJ) 

C  TLOC ( J J ) rCTLOC ( J J ) +SSCC 
TLQC4 l JJ)=TL0C4 ( JJ) +S6CC 
430  CONTINUE 
GO  TO  3'40 

167  IF  (  LS (NS)  .EQ,  1)  GO  TO  169 

LS(NS)rv.S(NS)-l 
GO  TO  lo6 

lo9  LS (NS) rLS (NS) +1 
GO  TO  166 

166  IF  (  LS(NS)  ,E0.  3  )  GO  TO  l6o 
LS ( NS ) r lS ( NS ) +2 
C 

166  CALL  UUMP(NL) 

C 

CTLOC ( NL ) =C  TLOC ( NL ) +TUUMP ( NL ) 

TL0C2 ( NL ) =TL0C2 ( Nl ) +TUUMP ( NL ) 

NLDL=NL0l.-l 

NLOE=NLuE+l 

IDE(NL)rNLDE 


C'.  n  n 
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WTO=WTI)+WTLOAU(NL) 

WrrKN(NL)=WTTi<N(ML)-WlLOAD(NL  ) 

WTLOAD(nU  =  0.0 

If:  (NLOL  .EQ.  1)  GO  Tu  170 

DO  171  i=l»NLOCO 

If-  (  lOL(I)  ,hO.  G)  GO  TO  171 

IDL(I)=1DL<I)-1 

171  CONTINUE 
GO  TO  172 

170  IUL(NL)=IDL<Nl)-1 

172  IF  (  IL  .£0.  o)  GO  TO  300 
DO  173  1=1 » IL 

JJ=U'.V  ( I ) 

LW ( I ) -0 

C  T LOC ( J J ) =CTloC ( J J ) +TOUMP ( ML ) 

173  TL0C4( JJ)=TL0C4( JJ)+TuUMP(MU 

1L=0 

GO  TO  300 

300  DO  17b  1=1 » NLOCO 

175  IF  (  CTLOC(I)  ,GE,  TlwAX  )  GO  TO  176 
IK  (  IDCOAD  ,tO.  OJ  GO  TO  301 

tttm  ,L2*-r1  *MND*  NL  ,LT*  NL°CO  )  GO  TO  301 

CTTMZCTIME  ♦  SCCTm 

IF  (  ILL  .EQ,  0)  GO  To  733 

DO  734  1  =  1 » NLOCO 

DO  73b  J=1»ILC 

70b  IF  (  1  .EG.  LC(J))  GO  TO  734 

IF  (  CTLOC(I)  ,LT.  CT l M )  GO  Tu  734 
ILC=ILC+1 
U(ILC)=I 
734  CONTINUE 
GO  TO  736 

733  DO  720  i=l»NLOCO 

IF  {  CTLOC(I)  ,LT.  CTTM)  GO  To  720 
ILC=ILC+1 
LC ( ILC ) =1 
720  CONTINUE 

736  IF  (  ILC  .EQ,  0)  GO  TO  301 
34T  IF  (  ILC  .GE.  NLOCO)  «o  TO  701 

FIND  THE  SMALLEST  VALUE  OF  CTLOC(I)  EXCEPT  FOR  CTLOC(lLC) 
M=1 

727  DO  72b  1=1 >ILC 

725  IF  (  M  .EG.  LC ( I ) )  GO  TO  726 
GO  TO  728 

726  M=M+1 

GO  TO  727 

728  Ml=M+i 

IF  (  Ml  .GT,  NLOCO)  Go  TO  841 
00  729  1=M1, NLOCO 
MJ=  I 
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DU  7.50  J=1*ILL 

730  IF  (  HJ  ,EU.  LC(J))  GU  TO  729 

IK  (  CTLOC(M)  .LE.  CTU)C(MJ))  GO  TO  72q 
M=MJ 

id 9  CONTINUE 
641  NL=M 

GU  TO  2u0 

176  PKINT  177 

177  FORMAT  (  1H0, ’SIMULATION  TERMINATED  BY  MAX  CLOCK  TIME  ALLOWED  * ) 
IUEOS=l 
GO  TO  7J1 

PKINT  THE  SUMMARY  OF  SIMULATION  AT  THE  END  OF  RUN 
299  PKINT  193. HAUL .WTG.WTu 

193  FORMAT  (  1H1, ’LENGTH  OF  HAULAGE  LlNE= * , F12 , 3/ *  »,*WT  OK  MUCK  GENER 
*ATCO=»,  F12.3/*  *.  *WI  OF  MUCK  DUMPED-’ .Fi2, 3  //) 

PKINT  194 

194  FORMAT  (1H1.  ’LOCO  NO.  CLOCK  TIME  LOADING  TIME  DUMPING  TIME  KUNlJ 
'KING  TlMc.  WAITING  TIME  UlSTANCt  TRAVELED’) 

DO  195  i=l»NLOCO 
CTLOC ( I ) =CTLOC < I ) /60 . u 
TLOCl(I)=TLOCi(I)/GO.O 
TL0C2(I)=TL0C2(I)/o0.u 
TL0C3 ( I ) s  TLOC J ( I ) /GO , 0 
TLOC4 ( I ) =TL0C4 ( I ) /60 . 0 

195  PRINT  196*  I  * tTLOC ( I ) *  TLOC1 ( I ) »TL0C2 ( I ) »TLOC3 ( I ) , TL0C4 ( I ) *DlSTR ( I ) 

196  FORMAT  (2X .  13. 3X * F  12.3 »  4(3X»M0.3)»  F15.3) 

IDEOS=l 

701  RETURN 
END 

compilation:  no  diagnostics. 


